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(57) Abstract: 
This invention 
describes the use 
of peptide profiling 
to identify, 
characterize, and 
classify biological 
samples. In 
complex samples, 
many thousands of 
different peptides 
will be present 

at varying concentrations. The invention uses liquid chromatography and similar methods to separate peptides, which are then 
identified and quantified using mass spectrometry. By identification it is meant that the correct sequence of the peptide is 
established through comparisons with genome sequence databases, since the majority of peptides and proteins are unannotated and 
have no ascribed name or function. Quantification means an estimate of the absolute or relative abundance of the peptide species 
using mass spectrometry and related techniques including, but not limited to, pre- or post-experimental stable or unstable isotope 
incorporation, molecular mass tagging, differential mass tagging, and amino acid analysis. 



Dc novo 

Identification of 
amino tod sequence 
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PROTEIN EXPRESSION PROFILE DATABASE 



CROSS REFERENCE TO RELATED APPLICATION 

This application claims priority from Canadian Patent application No. 2,349,265, 
which is incorporated by reference herein. 

FIELD OF THE INVENTION 

The field of this invention relates to the fields of peptide separation and 
proteomics, bioinformatics, metabolite profiling, medicine, drug screening and 
computer databases. - 

BACKGROUND OF THE INVENTION 

Modern biochemistry and molecular medicine is entering the post-genomic era. 
While genome sequencing has generated a large amount of genetic data, the 
focus in the biological sciences is now changing to the full characterization of 
proteins. Protein posMranslational modifications, protein localization, protein- 
protein interactions, and analysis of protein structure and folding have become 
subjects of major importance. 

Proteomics is the study of patterns of protein expression by complex biological 
sy stems. I t involves, in principle, the det ermi nation of the relative abundance, 
post-translational modification, and/or stability of large numbers of cellular 
proteins at specific time-points within the life cycle of an organism. 



There is growing recognition that qualitative and quantitative analysis of protein 
expression profiles on a genome-wide scale will accelerate the development of 
powerful new diagnostic tools and therapeutics, including novel biomarkers and 
drug targets, as well as lead to a better understanding of the basic molecular 
logic that governs cell biology. This is because most, if not all, complex biological 
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alone. This is because molecular regulation of proteins, and not simply their 
corresponding genes, holds the key to the function of most, if not all, complex 
biological processes. 

In contrast to genomics, which captures DNA information that is largely stable 
throughout the lifetime of an organism, proteomics efforts seek to summarize the 
protein-expression patterns of dynamic biological systems at different times. 
While there are a finite number of genes in a given genome, a cell's proteome is 
constantly fluctuating in response to environment and cellular perturbations. 
Hence, understanding how proteins work together requires systematic data on 
the entire spectrum of protein status in a cell at any given time. 

Biology Enters the Post-Genomic Era 

By the late 1990's the DNA sequences of numerous bacterial and eukaryotic 
organisms had been published and in 2000 the nearly complete DNA sequence 
of Homo sapiens was completed. The availability of large-scale genomic 
sequencing efforts now offers investigators a unique opportunity to perform 
comparative analysis from an evolutionary perspective which can both help to 
annotate and validate completed genome sequences and also help identify 
conserved protein function, regulation, or pathways based on protein sequence 
homology. 

Today several disciplines, in particular bioinformatics, functional genomics, and 
proteomics, are converging in efforts to exploit this newly-available genome 
sequence information. The long-term objective of these efforts is to understand 
the function and interrelationships of the many thousands of genes and proteins 
present in human cells, with the implicit expectation that this understanding will 
lead to dramatic progress in the clinical sciences. 
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large sets of proteins (Flores et al, 1999). inoluding whole yeast proteome (Ito et 
al.. 2000; Uetz et a,, 2000). Second, comprehensive soreening of mutant genetic 
loc, as a means for disserting networks of interacting gene products has recently 
been adapted to automated high-throughput formats. Finally, powerful 
experimental to* for identifying the components of protein samples, including 
large complexes such as the ribosome (Link et al., 1999) and nuctear pore (Rout 
et al., 2000). and most recently whole organelles and whole cells have been 
described. 



Tandem Mass S pe>ntrr>mo*ry 



Because the amino acid sequence of a protein is encoded in DNA, and because 
the rules for determining the primary amino acid sequence of a protein are 
known, vast numbers of hypothetical proteins with no known function await 
classification and characterization. Clearly, many of these genes and proteins 
Play a role in human disease and other phenomena of biologica. or commercial 
interest. 

The emerging Held of proteomics research relies on enabling technologies tha, 
can accurately and rapidly characterize the numerous diverse proteins typically 
found inbiological samples. This requires scalable, robust, and automated _ 
methods for protein analysis. 

To reveal biochemical pathways and regulatory networks, and help define new 
targets tor structure-function analysis, proteomics studies require high-resolution 
high-sensrtivrty techniques for separation, detection, and quantitation of proteins ' 
as well as methods for .inking proteins to their corresponding cognate gene 
sequences. 
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Mass spectrometry is the study of gas phase ions as a means to characterize the 
structures, and hence identities, of molecules. Proteomics began with the 
commercialization of soft ionization techniques in the 1990s, in particular 
electrospray ionization (ESI) and matrix assisted laser desorption ionization 
(MALDI), which permitted analysis of proteins for the first time. Commercial MS 
instruments are designed as high performance instruments for structural 
characterization of ions produced by these soft ionization techniques and have 
largely replaced traditional Edman chemical sequencing for the analysis of 
proteins. MS has proven to be very successful at identifying limited numbers of 
proteins, such as single polypeptide bands cut from polyacrylamide gels, and it is 
currently possible to identify proteins at picomolar to sub picomolar levels. 

Recent advances in mass spectrometry and data analysis described below are 
providing the necessary tools for implementation of high-throughput protein 
identification and characterization. As the scope of protein analysis has shifted 
from a molecule-by-molecule approach to a genomic scale, the ability of both 
academia and industry to generate new MS data has dramatically outstripped the 
ability to validate, manage, and interrogate the data. 

For these studies, routine access to state-of-the-art mass spectrometry 
instrumentation with an adequate infrastructure is essential. Two new ionization 
techniques, MALDI and ESI, have revolutionized the analysis of proteins. The 
MALDI and ESI techniques can be coupled with various types of mass analyzers, 
such as quadrupoles (Quad, Q), time-of-flight (TOF), ion-trap, Fourier transform 
ion cyclotron resonance (ICR) and hybrid instruments with two different mass 
analyzers (Q-TOF). Each kind of instrument has advantages and disadvantages 
and, in practice, the achievement of high throughput in conjunction with reliable 
protein identification requires access to both MALDI and ESI instruments. 
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separation techniques such as affinity chromatography, HPLC or capillary 
electrophoresis. 

Tandem mass spectrometry (MS/MS) provides a means for fragmenting a mass- 
selected ion and measuring the mass-to-charge ratio (m/z) of the product ions 
that are produced during the fragmentation process. The MS/MS process used 
most often is based on collision-induced dissociation (CID), in which a mass- 
selected ion is transmitted to a high-pressure region of the instrument where it 
undergoes low energy collisions with inert gas molecules. 

As a molecular ion collides, a portion of its kinetic energy is converted into excess C 
internal energy rendering the ion unstable, and driving unimolecular 
fragmentation reactions prior to leaving the collision cell. Detailed structural 
information is generated as a result of fragmentation. The mass selectivity of 
many commercial MS systems permit the isolation of single precursor peptide 
.ons from mixtures, thereby removing the contribution of any other peptide or 
contaminant from the sequence analysis step. The product ion spectra can 
subsequently be interpreted to deduce the amino acid sequence of a protein. 

A protein to be identified by MS is first digested enzymatically with a site-specific 

protease such as trypsin (which cleaves after lysine and arginine residues) in 

order to produce peptides with structures suitable for MS. Tryptic peptides are C 

particularly amenable to MS/MS analysis since mobile protons localize to the N- 

terminal amine and the side chains of the carboxy-terminal arginine or lysine 

res,dues at which proteolysis occurs. These protons cause peptides to fragment 

in a somewhat predictable manner following activation in a tandem MS leading 

to production of two broad classes of fragment ions - the so-called amino- 

terminal b-type ions and carboxy-terminal y-type ions. Recognition of the 

members of these series is a fundamental process of MS-based protein 
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trypsin digestion of protein sample is initially determined by the mass 
spectrometer. The peptides are then isolated based on their mass/charge 
properties, fragmented using low energy collision with inert gas (or with 
resonance excitation), and the fragments are analyzed using a second round of 
mass spectrometry. 

The relative abundance of daughter product ions in peptide tandem mass spectra 
varies considerably, and some are not observed. This variation reflects subtle 
differences between favored and disfavored fragmentation sites, the nature of the 
amino acid side chains, and their position on the peptide backbone. CID of 
protonated peptides also leads to other fragmentation reaction products that can 
complicate spectral interpretation. Molecular losses of water or ammonia for 
instance, are commonly observed in the product ion scans of tryptic peptide ions. 
Spectra often also contain non-peptide noise peaks. Because of this, de novo 
interpretation of spectra is extremely difficult to automate and most MS-based 
identification techniques rely on reducing the computational scale of the problem 
by searching protein sequence databases using a relatively simple correlation 
algorithm. 

The fragmentation patterns of the peptides can be used to obtain amino acid 
sequence information by comparison with predicted patterns obtained from 
translated protein databases. In addition, advances in tandem mass spectrometry 
mean that polypeptides can now be identified at a low picomolar to femtomolar 
level in a rapid, sensitive, and versatile manner. By revealing the composition of 
biologically relevant, low abundance protein complexes, the technology can 
provide fundamental insight into the circuitry of interacting proteins. 

Tryptic peptides are particularly amenable to MS/MS analysis since mobile 
protons localize to the N-terminal amine and the side chains of the carboxy- 
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(a typical MS/MS peptide spectra showing prominent b- and y-ions is shown 
below). 

The fragmentation pattern reflects the dissociation of the peptides along the 
pept.de bond backbone, and therefore correlates with the sequence of amino 
acds for those peptides. Recognition of the members of the b- and y-ion series 
is a fundamental process of MS-based protein sequence interpretation. Since de 
novo interpretation of spectra is difficult to automate, most MS-based 
identification techniques rely on reducing the computational scale of the problem 
by searching protein sequence databases using a relatively simple correlation 
algorithm. The SEQUEST program (US Patent 5,538,897), for instance uses 
uninterpreted product ion spectra to search databases of theoretical spectra 
derived from protein and translated gene sequence databases. 

Recent developments in tandem mass spectrometry (MS/MS) now allow for the 
identification of hundreds of proteins per samp.e in a single run using available 
technology. This represents a major breakthrough compared to traditional 
methods, for example, 2D gel electrophoresis, and permits, for the first time, 
protein analysis on a truly proteomic scale. 

Accurate mass measurement of p eptide s derived from proteins provides 
information not available from DNA sequence, such as post-translational 
modifications and correction to errors in the DNA databank. Database searching 
with masses of peptides obtained from proteolytic digests is a well-established 
technique in many laboratories around the world. The searching of databases 
with partial sequence information obtained from MS/MS sequencing experiments 
is even more reliable because it imposes statistical constraints on the 
identification. 
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molecular mass to one of the protein samples. By combining the samples after 
this treatment, the relative abundance of different protein species in each sample 
can be estimated by comparing the signal intensities of the corresponding 
peptides in the mass spectrometer. 

Another quantitative approach, limited to cultivable organisms, is to label growth 
media with stable isotopes such as N15. The isotope becomes incorporated into 
the peptide or protein and the isotope-treated peptide is offset in the mass 
spectrum by multiples of 1 amu (the difference in mass between the naturally 
abundant isotope N14 and the heavy isotope derivative N15) depending on the 
number of N atoms in the peptide. These spectra can be deconvoluted to 
determine the relative abundance of the labeled and unlabeled peptide species. 
Alternatively, non-isotopic mass tags, whereby the labeled' or tagged species is 
offset by the mass of the tag, can be used. Thus methods suitable for high- 
throughput and efficient identification and quantitation of large numbers of 
proteins from complex mixtures are now available. 

HPLC 

High-resolution separation techniques are required to separate the peptide 
components of complex biological mixtures prior to mass spectrometry. A 
particularly powerful approach to identifying the components of complex protein 
mixtures is direct analysis of the protease-digested proteins using high- 
performance, high-resolution multi-dimensional liquid separation techniques 
coupled online to mass spectrometry/database searching (HPLC-MS/MS)(Link et 
al., 1999). This strategy enables the separation of very complex peptide mixtures, 
such as the whole cell extracts or nuclear extracts (Washburn, 2000). One aspect 
of the method separates complex peptide mixtures by strong cation exchange in 
the first dimension and by reverse phase in the second. However, many 
combinations of separation media and more than two dimensions could be used. 
One advantage of the strategy is that it eliminates the need to separate proteins 
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Bioinformatics 

The interpretation of peptide mass spectra for the purposes of generating protein 
identifications can be carried out manually but requires experience and skill and 
is prohibitively time-consuming. For this reason, computer algorithms have been 
developed that, while not capable of interpreting all spectra they encounter, can 
easily outperform human identifications for even minimally complex peptide 
mixtures. Any of several generally available algorithms may be used for this 
purpose. For instance, the SEQUEST program (Eng et at. 1994) uses 
uninterpreted product ion spectra to search databases of theoretical spectra 
derived from protein and translated gene sequence databases, SEQUEST first 
generates a list of theoretical peptide masses for each entry in the database that 
match the experimentally determined peptide mass, producing a list of candidate 
peptides. The program then calculates the fragment ion masses expected for 
each of the candidate peptides, generating a predicted MS/MS spectrum. Finally 
the experimentally determined MS/MS spectrum is compared with the predicted ' 
spectra using a correlation function. Each comparison receives a score, and the 
highest-scoring peptide(s) are reported. When high scoring matches are 
detected, one effectively jumps from spectral data directly to a peptide identity, 
which in turn can be linked to the entire amino acid and DNA sequence of the 
corresponding gene. Ideally, a protein is positively identified when the spectra of 
one or more peptides in a tryptic digest can be matched unambiguously. 

Mass spectral reference libraries representing stored tandem mass spectra, or 
validated, chemical signatures, are routinely used for the identification of small 
chemical compounds by MS (eg. Wiley Registry, N 1ST database). Unknown 
compounds can then be both identified by searching experimental spectra 
against a comprehensive database of these reference mass spectra, which are in 
turn derived from pure compounds, so that only hits of strong similarity or identity 
are produced. A similar reference spectral database approach would likewise 
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holds great promise for rapid genome functional analysis. It is plausible that the 
protein expression profile could serve as a universal and rich cellular phenotype: 
provided that the cellular response to disruption of different steps of a given 
biochemical process or pathway is similar, and that there are sufficiently unique 
cellular responses to the perturbation of most cellular pathways, systematic 
characterization of novel genetic mutants could be carried out with a single 
genome-wide protein expression measurement. 

To date the only studies focusing on peptides or proteins that includes a 
quantitative component has been the separation of bacterial and yeast cell 
lysates on 2-dimensional electrophoretic gels (refs). These approaches do not 
directly identify the resolved proteins, are relatively insensitive, and are unlikely to 
scale up to the study of larger proteomes (e.g. that of vertebrates). Furthermore, 
no attempt was made to use the data to identify or characterize unknown 
samples. 

SUMMARY OF THE INVENTION 

The protein profiling approach proposed has both a qualitative and a quantitative 
component such that each profile generated can be directly compared to other 
profiles present in a reference database. 

This invention describes the use of peptide profiling to identify, characterize, and 
classify biological samples. In complex samples, many thousands of different 
peptides will be present at varying concentrations. The invention uses liquid 
chromatography and similar methods to separate peptides, which are then 
identified and quantified using mass spectrometry. By identification it is meant 
that the correct sequence of the peptide is established through comparisons with 
genome sequence databases, since the majority of peptides and proteins are 
unannotated and have no ascribed name or function. Quantification means an 
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The principle experimental strategy of the present invention is centered on rapid 
h.gh-throughput protein identification using coupled tandem mass spectrometry 
(MS/MS) and sequence database searching. Quantitation is based on either 
metabolic labeling with stable isotopes or with chemical derivation. Below an 
example of a non-isotopic tag based on the lysine-specific guanidylation reagent 
O-methy.isourea is described in detail. Significant patterns of peptide expression 
are .dentified with software and data mining algorithms. Below, a method is 
described for identifying, classifying and characterizing functions of known and 
unknown gene products, peptides and proteins , for characterizing metabolic and 
other functional pathways in cells, and for identifying the proteins and pathways 
targeted by drugs and other reagents. The method is based on the comparison 
of protein profiles obtained following global proteomics or other comprehensive 
protein studies from cells, cell fractions, tissues, organisms or other defined 
sources. 



The invention further contemplates the use of high-throughput robotic screening 
of diverse chemical compound libraries to systematically identify small molecules 
that perturb cellular pathways associated with disease. The protein targets of the 
lead compounds will be isolated and identified by the tandem mass spectrometry 
profiling techniques described herein. Protein profiling acts as an optimal assay 
since the profile of a healthy cell or tissue is the goal. _ 

The invention relates to a method for identifying the constituent proteins for a cell 
type, tissue or pathological sample using a database comprising peptide profile 
libraries wherein the libraries have multiple peptide sequences, comprising- 

1 . deriving a plurality of peptides from the cell type, tissue or pathological 
sample; 

2. identifying the peptide species by liquid phase tandem mass spectroscopy 
sequencing; 

3. compiling a data set or peptide profile containing the collection of peptide 
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a) obtaining a peptide-containing extract of the cell type, tissue or pathological 
sample; 

b) digesting the extract producing peptides with an enzyme, the enzyme 
capable of localizing mobile protons to the N-terminal amine and the side 
chains of the carboxy-terminal arginine or lysine residues; 

c) separating the peptides by high pressure liquid chromatography apparatus; 

The enzyme preferably comprises one selected from the group consisting of 
trypsin and endoproteinase LysC. The step of digesting the extract producing 
peptides preferably further comprises the steps of: 

a) dividing the extract into two equal portions; 

b) derivatizing completely one of the two equal portions with a reagent, the 

reagent comprising one selected from the group consisting of o- 
methylisourea, homoarginine, canavanine, hydrazine, phenylhydrazine, and 
butyric acid derivatives. 

c) combining the two portions. 

The methods of the invention may be used in toxicology analysis. The methods 
optionally comprise administering a candidate compound to a cell. As described 
above, samples suitable for MS anaylsis are generated and a peptide profile is 
produced. Relative abundance of peptides in samples is also preferably 
determined. This candidate compound peptide profile is compared to peptide 
profiles in a database or library (for example, profiles showing the cell in a normal 
state and in varied states of toxicity). If the candidate compound sample profile is 
highly similar to (for example, greater than 90%, 95%, or 99% similarity), or 
identical to a profile in the database or library, then that similarity shows the 
amount of toxicity of the candidate compound to the cell. If the candidate 
compound sample profile is highly similar to a normal cell profile, then the 
candidate compound is less likely to be toxic than if the candidate compound 
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and relative abundance towards a normal, healthy profile and relative abundance 
wth substantial similarity (eg. Over 90%, 95%, 95% similarity), or identical to the 
healthy profile and relative abundance, the drug compound is likely to be useful 
as a therapeutic. 



Another embodiment relates to a method for identifying a peptide sequence for a 
cell type, tissue or pathological sample using a database comprising peptide 
profile libraries wherein the libraries have multiple peptide sequences, 
comprising: 

a) obtaining a peptide-containing extract of the cell type, tissue or pathological ~ 

sample; 

b) digesting the extract producing peptides with an enzyme capable of localizing 

mobile protons to the N-terminal amine and the side chains of the carboxy- 
terminal arginine or lysine residues; 

c) separating the peptides by high pressure liquid chromatography apparatus- 

d) identifying the peptide species by tandem mass spectroscopy sequencing- ' 

and 

e) compiling a data set or peptide profile containing the collection of peptide 

sequences obtained thereby. 
The enzyme is preferably selected from the group consisting of trypsin and 
endoproteinase LysC. The step of digesting the extract producing peptides 
preferably further comprises the steps of: 

a) dividing the extract into two equal portions; 

b) derivati2in9 completely one of the two equal portions with a reagent, the 
reagent comprising one selected from the group consisting of o- 
methylisourea, homoarginine, canavanine, hydrazine, phenylhydrazine, and 
butyric acid derivatives. 

c) combining the two portions. 
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a) deriving a plurality of peptides from each sample of the cell type, tissue or 
pathological sample; 

b) identifying the peptide species by tandem mass spectroscopy sequencing; 

c) compiling a data set or peptide profile containing the collection of peptide 
sequences obtained thereby; 

d) cross-tabulating with a collection of peptide sequences in the database of 
peptide sequences; and 

e) determining the relative abundance of the proteins. 

In the methods of the invention, a pathological sample may have been contacted 
with a candidate drug compound and the peptide profile and/or relative 
abundance of the peptides and/or proteins is compared to a database comprising 
peptide profile libraries of the cell in varied states of toxicity (ie. exposed to known 
toxic compounds which injure and/or kill the cell). The toxicity of the candidate 
drug compound may be determined by comparison of the profile and relative 
abundance for the cell type, tissue or pathological sample exposed to the 
candidate drug compound with the profile and relative abundance for the cell 
type, tissue or pathological sample in varied states of toxicity and a normal state. 
A similar method may be used to determine whether a compound is likely to be 
useful as a therapeutic, for example by comparison of the profile and relative 
abundance for a pathological (diseased) cell type, tissue or sample exposed to 
the candidate drug compound with the profile and relative abundance for the cell 
type, tissue or sample in a normal, healthy state. 

The invention includes a method for quantitating the relative abundance of 
proteins in two samples of a cell type, tissue or pathological sample using a 
database comprising peptide profile libraries wherein the libraries have multiple 
peptide sequences, comprising: 

a) deriving a plurality of peptides from each sample of the cell type, tissue or 
pathological sample; 
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d) determining the degree of relatedness of a collection of peptide sequences 
■n the database of peptide sequences using clustering and related statistical 
methods 

The step of deriving a plurality of peptides in two samples preferably further 
comprises the step of: 

a) obtaining a peptide-containing extract of each sample; 

b) digesting separately the extracts producing peptides with an enzyme, the 

enzyme capable of localizing mobile protons to the N-terminal amine and 
the s.de chains of the carboxy-terminal arginine or lysine residues; 

c) combining the two extracts; and - - - . - 

d) separating the peptides by high pressure liquid chromatography 

The enzyme preferably comprises one selected from the group consisting of 
trypsin and endoproteinase LysC. 

The step of digesting the extracts preferably further comprises the step of 
derivatizing completely one of the two extracts with a reagent, the reagent 
comprising one selected from the group consisting of o-methylisourea 
homoarginine, canavanine, hydrazine, phenylhydrazine, and butyric acid 
derivatives. 

The invention a.sc includes a method for identifying a peptide sequence for a cell 
type, tissue or pathological sample, comprising: 

a) obtaining a peptide-containing extract of a cell type, tissue or pathological 

sample; 

b) digesting the extract producing peptides with an enzyme capable of localizing 

mob.le protons to the N-terminal amine and the side chains of the carboxy- 
terminal arginine or lysine residues; 

c) separating the peptides by high pressure liquid chromatography apparatus- 

d) .dentifying the peptide species by tandem mass spectroscopy sequencing; ' 
and 
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The step of digesting the extract producing peptides preferably further comprises 
the steps of: 

a) dividing the extract into two equal portions; 

b) derivatizing completely one of the two equal portions with a reagent, the 
reagent comprising one selected from the group consisting of o- 
methylisourea, homoarginine, canavanine, hydrazine, phenylhydrazine, and 
butyric acid derivatives. 

c) combining the two portions. 

Another embodiment of the invention is a computer system for identifying 
quantitative peptide profiles, comprising: 

(a) a database including peptide profile libraries for a plurality of types of 
organisms wherein the libraries have multiple peptide profiles each profile 
comprising an array of at least 50 peptide species each having a unique 
identifier cross-tabulated with quantitative data indicating relative and/or 
absolute abundance of each peptide species in a sample; and 

(b) a user interface capable of receiving a selection of one or more queries to 
the database for use in determining a rank-ordered similarity of peptide 
profiles in the database. 

The invention includes a method of producing a computer database comprising a 
computer and software for storing in computer-retrievable form a collection of 
peptide profiles for cross-tabulating with data specifying the source of the 
peptide-containing sample from which each peptide profile was obtained. 
Optionally, at least one of the sources is from a sample known to be free of 
pathological disorders. Optionally, at least one of the sources is a known 
pathological specimen. 

The invention also includes a method of comparing quantitative peptide profiles 
using a database of a plurality of peptide profile libraries, the method comprising: 
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The correlation of a peptide profile against selected peptide profile libraries may 
be determined by 

Px.y = [1/n (p! to „) I (Xj - Px ) (Yj - Uy)]/ia x . 3 y ] 

where peptides common to two profiles score T and peptides not shared 
between profiles score '0\ 

The peptides profiles are preferably of cell fractions, the cell tractions comprising 
h,gh molecular weight proteins, soluble proteins, membrane proteins, modified 
proteins, phosphoproteins, peptides terminating in lysine o, arginine or the 
specific products of proteolytic enzymes or chemical derivatives of those 
products, peptides containing rare amino acids, and proteins isolated by binding 
to disease-specific affinity reagents. * ( 

The specific products of proteolytic enzymes may be comprise chemical 
derivatives of these products wherein de novo sequencing or relative abundance 
measurements of the peptides is facilitated. 

The chemical derivatives may be obtained by guanidinylation and related 

mediations. The rare amino acids may comprise tryptophan and cysteine and 

ammo acids comprising 50/0 or less of the amino acid representation 

The disease-specific affinity reagents may comprise polyclonal antibodies, toxin 

or drugs. The peptide profiles may be of peptide sequences, the peptide 

sequences comprising mammalian peptide sequences. Thee peptide profiles 

may be of peptide sequences, the peptide sequences comprising microbial 

peptide sequences. ^ 

The step of receiving a selection of two or more of the peptide profile libraries for 
comparison may include receMng-a user selection from two or more pulldown 
menus using a graphical user interface. The step of receiving a selection of two 
or more of the peptide profite libraries for comparison may comprise command 
l.ne entry using a computer. The step of receiving a selection of two or more of 
the peptide profile libraries for comparison may comprise receiving an 
electronically transmitted file containing sequence and quantitative data The 
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database. The method may further comprise the step of displaying the peptide 
profiles common to the selected peptide profile libraries. The method may further 
comprise the step of displaying the peptide profiles unique to the selected peptide 
profile libraries. 

The invention also includes a method of identifying peptide profiles common to a 
set of environments, organisms, organs, tissues, cells, cellular fractions or 
isolated molecular complexes using a database comprising peptide profile 
libraries for a plurality of types of organisms wherein the libraries have multiple 
peptide sequences, the method comprising: 

(a) displaying at least one list of peptide profile libraries; 

(b) receiving a selection of one or more peptide profile libraries from at least 
one list of peptide profile libraries; 

(c) determining peptide profiles common to the selected peptide profile 
libraries; and 

(d) displaying the results of said determination. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Embodiments of the invention will be described by way of example and 
with reference to the drawings in which: 

FIG. 1 is a diagram of the MCAT approach for peptide sequencing and relative 
protein abundance determination. 

Fig. 2 is diagram showing how MCAT enables identification and quantitation of 
complex protein mixtures. 

Figs. 3A and 3B are diagrams showing de novo sequencing of a yeast peptide 
and a human peptide using MCAT approach. 

Figs. 4A and 4B are diagrams showing relative abundance ratios of positively- 
identified peptides. 

Fia. 5 is a oeotide Drofile aenerated hv a nnR-Him^nRinnai I ckar fmm Hiv/crco 
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Fig. 7 shows .he differences between protein expression of the seven human 
hssues h,ghhghted by applying agglomerate clustering algorithms 
Frg . • is a similarity dendrogram for different human tissue consumed using 
peptide profiling. y 

Fig. 9 is a comparison of peptide profiles of different cell compartments 

22 ^s COmPariSOn °' PeP,We Pr0ffleS Umrea,ed 3nd '~ ated — 

Fig. 11 shows peptide profiling to distinguish species 

Fig. 12 is a representation of a reference database of protein profiles. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

A Quantitative Peptide Profile serves as a precise fingerprint o, peptides that can 
be successful isoiated. identified and quantified from the myriad of prote n s 
expressed in cells under any given condition. This profile, in turn, can serve as a 
un, que .denfifier of eel, state. This document describes a method to use 
quantitative pepfide profiles ,o compare biological samples, from any tissue or 
cel., among afferent types of cel. (e.g. nervous tissue cells,, or even in samples 
where little or no mRNA is made (e.g. blood platele, cells). ' 

The present invention is distinct from the established method of mRNA 
exp ression profiling in three important respects. 

First as menfioned above, the relative abundance of an mRNA is not predictive 
°«he abundance of .^corresponding p ro , e in or cognate peptides. W Is 
because many factors affec. protein expression subsequent*, me even, of 
mRNA producon, including splicing, protein terminal processing, protein 
loca„ 2 a.,on, protein degradation, protein modification, codon usage the levels of 
ava„ab,e amino acids and the subce„u,ar localization of «he prote n' mRNA 
express,on profiling is unable ,o account for or predict these events 
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includes a method for detecting and quantitatively analyzing peptides in a 
biological sample, comprising: 

a) obtaining a biological sample in a form suitable for coded abundance tagging; 

b) identifying and quantitating the peptides in the sample by mass coded 
abundance tagging. 

In one aspect, the method involves: 

obtaining an extract of the biological sample, such as a cell extract, 
digesting the sample, preferably with an enzyme, such as trypsin, to generate 
peptides with a terminal amine group, such as a terminal lysine, 
contacting the peptides with mass differential reagent, such as a guanidination 
compound (eg. Lysine guanidination compound, such as o-methylisourea, which 
modifies the epsilon-amine of the C-terminal lysine\), 

separating the peptides, preferably with liquid chromatography, such as high 
throughput capillary liquid chromatography, and 

generating mass spectra for the peptides, preferably with electrospray tandem 
mass spectrometry. 

The method is preferably carried out in both orientations, with a sample divided in 
two and either modified or unmodified. Peptides are alternatively unmodified and 
modified with o-methylisourea differ by the mass differential encoded by the mass 
differential reagent (e.g. 42 amu for O-methylisourea). The method preferably 
further involves sequencing the peptides and/or determining relative abundance 
of the peptides. Methods of sequencing and determining relative abundance are 
described below. Sequencing preferably involves comparing pair-wise sets of 
spectra (MS/MS spectra) to identify identities of y-ion peaks. One can use a 
short sequence of contiguous amino acid sequence from a peptide (e.g. 5-10 
amino acids or greater than 10 amino acids) to identify a corresponding protein. 
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The invenBon includes a method of identifying a test sample by obtaining a 
peptide profile for the tes, sample, preferably by MS. This peptide profile is then 
compared to peptide profiles in a database or library to determine if the test 
sample profile is highly similar to (for example grater than 90"/., 95% or 99% 
s,mi«arity> ,0 a profile in the database or library. Re.ative abundanoe information 
may similarly be used to identify the test sample. 

The methods of the invention may be used in toxicology analysis. The methods 
opbonally comprise administering a candidate compound to a cell. As described 

above, samp.es suitable for MS anaylsis are generated and a peptide profile is 

produced. Relative abundance of peptides in samples is also preferably ( 
defended. This candidate compound peptide profile is compared to peptide 
profiles in a database or library (for example, profiles showing the eel, in a norma, 
state and in varied states oftoxicfty). »«. candidate compound sample pro* is 
highly s,m,lar to (for example, greater than 90%. 95%, or 99% similarity) or 
identical to a profile in the database or library, then that similarity shows 'the 
amount of toxicity of the candidate compound to the cell. If the candidate 
compound sample profile is highly similar to a normal cell profile then the 
candidate compound is less likely to be toxic than if the candidate compound 
sample profile is similar to the peptide profile of the cell in state of toxicity The 
relative abundan ce of th e tes, sa mple peptides is a,so prefera bly compart to 
othrer profiles to determine the amount of toxicity of a candidate compound " In a < 
s,m,lar manner, candidate drugs compounds may be screened against cells, such 
as diseased cells. If the candidate drug shifts the profile from a disease profile 
and relative abundance towards a normal, healthy profile and relative abundance 
w,th . substantia, similarity (eg. Over 90%. 95%, 95% similarity,, or identical ,0 the 
healthy profile and relative abundance, the drug compound is very likely ,0 be 
useful as a therapeutic. 



SOOCID: •cWO_02097703A2_L> 



21 

SUBSTITUTE SHEET (RULE 26) 



WO 02/097703 PCT/CA02/00801 

Using a comprehensive database of reference peptide expression profiles, the 
pathway(s) perturbed as a consequence of an uncharacterized mutation, 
pharmaceutical treatment, or developmental or disease state would be 
ascertained by simply asking which expression patterns in the database the 
resulting profile most strongly resembles. The database or library will include one 
or more profiles and/or relative abundance determination and may be electronic 
or in a hard copy form. A sufficiently large and diverse set of profiles obtained 
from different mutants, chemical treatments, and environmental conditions would 
also result in a relatively comprehensive identification of coordinate protein 
expression sub patterns, allowing hypotheses to be drawn regarding the 
functions of gene products based on their relationship to other proteins (Eisen et 
al., 1998). 

There are several advantages to this profiling approach compared to the analysis 
of single peptides or proteins. First, there is no requirement for prior knowledge 
about the functions of the responsive peptides or parental proteins. Second, 
protein functions deduced from comparisons of profiles in a database can be 
derived from very subtle physiological responses. For instance, even though 
peptide levels may change only slightly in response to an experimental treatment, 
coordinate changes among many measured peptide abundances can be 
sufficient to characterize that phenotype. The large numbers of peptides 
measured make it unlikely that an unrelated physiological state will have an 
identical profile, even though this may not be apparent when using conventional 
experiments that measure the levels of one or a few proteins. Third, closeiy 
related profiles can be classed together, thus improving our understanding of the 
underlying biological basis of the classifications. 

The invention includes proteins, including drugs, and other compounds identified 
using methods of the invention. 
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Examples 

Example 1: Measurement pf r ^ ein mlMum ^.^ncr 



in complex mixtures 



The method relies on modification of peptides at e . amlne of |ys|ne ^ 
O-methylisourea. Peptides so modified can be readily detected by mass 
spectrometry because their mass is increased by 42Da (per lysine residue in the 
sequence). Therefore, the reiafrve abundance of a single peptide from two 
drfferent samp.es can be determined following differential modification wilh O- 
methylisourea by comparing** signal intensifies forthe pair in a mass " 
spectrometer. ( 



The steps of the MCAT procedure are as follows (Fig. 1): 

(D Two protein mixtures, obtained following different experimenta. treatments 
of a sample, are digested enzymatically with trypsin. 

(2) One digest is treated with O-methylisourea and the other with control 
buffer. 

(3) The digests are desalted using ZipTip reverse phase extraction 

(4) The two mixtures are combined and analyzed by automated e.ectrospray 
LC-MS/ MS. Using.either one-dimensional (reverse phase) or two- 
d.mensiona. (cation exchange and reverse phase) liquid chromatography 
the peptides are separated as they are introduced to the mass 
spectrometer. The instrument is run in automated multistage mode 
whereby the following cyole is implemented. First, a full MS scan (400- 
1600 m/z) is used to record the relative intensities of peptide ions 
emergmg from the column. Next. MS/MS scans of selected ions are used 
to collect spectra suitable for peptide identification. The instrument then 
reverts back to full scan mode, but is programmed to exclude MS/MS 
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(6) For identified peptides, the single ion intensity profile is reconstructed from 
the full scan data and the relative abundance of modified and unmodified 
peptides calculated by integrating the area under the curve. 

In order to correct for systemic errors, for instance preferential labeling by O- 
methylisourea of one sample, the experiment is carried out in both orientations, 
that is both samples are divided in two and either modified or unmodified. The 
fractions are then combined with the corresponding modified or unmodified 
fracton from the other sample. 

Table 1 shows some top scoring peptides from this analysis and their relative 
abundance as estimated by the area-under-curve of their respective selected ion 
tracings. For nearly all peptides, the ratio of unmodified to modified signal is 
slightly less than the expected 1:1. The variation from ideal 1:1 ratio is not the 
result of reduced ionization efficiency or MS signal of the modified peptides 
relative to their unmodified forms because the effect was consistently observed in 
subsequent experiments independently of which sample was chosen for 
modification. More likely, it results from preferential recovery of unmodified 
peptides during the Zip Tip desalting step. 

For this reason, when comparing two samples A and B using the MCAT 
procedure, four mass spectrometry analyses are routinely carried out: I) A versus 
A mod , II) A versus B mod , III) B versus B mod , and IV) B versus A mod . The ratios of 
unmodified to modified peptide signals obtained in I and III were used to 
normalize II and IV respectively, and the combination of III and IV served to 
independently confirm the quantitative observations. 



ISDOCID: <WO_02O97703A2_l_> 



24 

SUBSTITUTE SHEET (RULE 26) 



WO 02/097703 

PCT/CA02/00801 



1. Identification and quantitation of peptides from a yeast whole cell digest 



I Protein 


Peptide 


Z a 


I Score 6 


Observed 
ratio 


Expected 
ratio 


YLR044C 


~ aqyneiqgwdhlsllptf" 
gak (seq id no:1) 


2 


2.3993 


1:0.29 


1:1 


YLR044C 
YLR044C 


ttyvtqrpvylglpan lv 

DLNVPAK (SEQ. ID N02> 


2 


2.6639 


1:0.2 


1:1 


KLIDLTQFPAFVTPMGK 
(SEQ ID NO:3) 




3.3881 


1:0.67 - 


1:1 


J YHR174W 


WLTGVELADMYHSLMK 
(SEQ ID NO:4) 




I A ftcro - 


1:0.73 


1:1 ~~ 


YHR174W 


GVMNAVNNVNNVIAAAFV 
K (SEQ ID NO:5) 


O 
£. 




1 .0.48 


f 1:1 


f YBR118W 


TLLEAIDAIEQPSRPTDKP 
LRLPLQDVYK (SEQ ID 
NO:6) 


o 
O 


o oooo 
o.oooo 


1 .0.63 


1:1 


YBR118W 


VETGVIKPGMWTFAPAG 
VTTEVK (SEQ ID NO'7) 


2 


2.5458 


1:0.23 


1:1 


YEL034W 


VHLVAIDIFTGK (SEQ ID 
NO:8) 


1 


3.0798 


1:0.15 


1:1 


YKL060C 


SPII LQTSNGGAAYFAGK 
(SEQ ID NO:9) 


2 


3.6709 


1:0.73 


11 


YCR012W 


ALENPTRPFLAILGGAK 
(SEQ ID NO: 10) 


2 


2.7650 


1:0.33 


1:1 


YDR441C 
| YGR192C 


GFVPIRRVGKLPGEC* 
(SEQ ID NO:11) 


2 


1.1770 


1:1.07* 


1:1 


VINDAFGIEEGLMTTVHSL " 


2 ~ 


3.1456 


1?031 h 
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Next, mixtures derived from yeast whole cell extracts containing varying 
proportions of MCAT-treated and MCAT-untreated sample were analyzed (Fig. 
2). 

Relative abundance signal from five peptides with high SEQUEST scores 
showed linearity across two orders of magnitude (Fig. 2). Beyond this range, the 
weaker signal of the two abundances is indistinguishable from background noise. 

Table 2 shows variation in the measured relative abundance for two peptides 
from the same parent protein (and therefore are present in equimolar 
concentrations) in three replicate experiments. Experiment-to-experiment 
variation for these peptides is within 25% and variation within a single experiment 
for peptides derived from the same protein is within 20% (Table 2). 

Table 2. Identification and quantitation of two peptides derived from YLR044C in 
three replicate experiments (A, B, C). 



Protein 


Peptide 


Ratio 
A:A 


Ratio 
A:B 


Ratio 
A:C 


YLR044C 


KLIDLTQFPAFVTPM 
GK (SEQ ID NO:13) 


1.00:1.00 


1.00:0.78 


1.00:0.87 


YLR044C 


AQYNEIQGWDHLSL 
LPTFGAK (SEQ ID 
NO: 14) 


1.00:1.00 


1.00:0.79 


1.00:1.03 



Ratio of unmodified to modified peptides (normalized to A:A) 



This invention also includes computer systems including software and hardware 
to implement the above methods. Such systems include a database with the 
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Example 2: De Novo Peptide Sequencing and Quantitative Profiling of Complex 
Prote,n Mixtures Using Mass Coded Abundance Tagging 

Introduction 

There is growing recognition the, qualitative and quantitative analysis of proteins 
on a genome-wide sca.e will accelerate the development of powerful new 
d,agnos„c tools and therapeutics, and lead to a better understanding of the 
molecular logic that governs cet. behavior. This is because regulation o protein 
abundance holds the k ey to the proper funcuon of most biological processes 
(Pandey & Mann. 2000,. Proteomics studies depend on sca,ab,e, robust and 
automated methods for protein identification and quantitation that can routinely 
characters the numerous diverse proteins typically found in biologica! samptes. 

Mass spectrometry (MS) is currentty the technology of choice tor identifying 
protems present in biologica, mixtures. The primary advantages of MS are its 
h,gh sensitivity, accuracy and capacity. Tandem mass spectrometry (MS/MS, 
proves a means for fragmenting mass-seleCed precursor peptide ions and 
measunng the mass-to-charge ratio (m*, of any product daughter ions produced 
Andersen e, a,.. 1996). The process usually produces two principte classes of 
ragmen, ,ons. the so-called M-termina, b-type ions and Carmine, y-type ions 

IZZV'V^ ~ " ™° show 

MS^tS ana,ys,s s.nce mobile protons the, stimulate the fragmentation process 
readily associate with the side chains of thP r tor^i^^i 

at which protects occurred ^ reSjdUeS 

If accurate sequence information is available, computer database search 
a gonthms can rapidly and accurate.y identify proteins analyzed by MS/MS (Eno 
et al., 1994; Mann & Wilm 1994- T M i nr a i„k ,™ 9 

vv,,m ' 1994 . Taylor & Johnson, 1997, Qin et al., 1997) in 
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proteomes of higher organisms, a facile peptide sequencing method that is 
independent of sequence databases is desirable. 

Manual interpretation of peptide MS/MS spectra for the purposes of protein 
identification (a process usually referred to as de novo sequencing) is often 
prohibitively challenging. Factors such as variation in favored fragmentation sites, 
the effects of the chemical nature of the amino acid side chains and their relative 
order in a peptide backbone, and the presence of side-products such as neutral 
loss ions and non-peptide noise peaks. To address this issue, Mann and 
coworkers pioneered a post-experiment stable isotope labeling strategy whereby 
the C-termini of tryptic peptides are labeled with deuterated water in order to 
reduce spectral complexity. Comparison of the modified and unmodified peptide 
MS/MS product ion spectra allows the C-terminal y-ions to be readily 
distinguished and, hence, the peptide sequence discerned. The impact of this 
approach has been restricted, however, by the prohibitive cost of the stable 
isotope and the high mass resolution required to distinguish the labeled products. 

Functional genomics studies using DNA microarray technologies have been used 
successfully to compare the abundance of thousands of mRNA species from 
distinct cell states. In contrast, only limited analogous quantitative data has been 
obtained for protein abundance. As the scope of protein analysis has shifted from 
a molecule-by-molecule approach to a genomic scale, the ability to generate 
quantitative protein data has lagged considerably. Chait and coworkers reported 
the potential of stable N 15 isotope labeling of proteins as a means to determine 
the relative abundance of select subsets of proteins isolated from cultured yeast 
cells (Oda et al., 1999). As the isotope becomes incorporated, the mass of the 
protein becomes offset in a mass spectrum by multiples of 1 amu (the difference 
in mass between the naturally abundant N 14 isotope and the heavy N 15 isotope 
derivative) depending on the number of labeled N atoms. Although powerful, this 
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containing moiety ,o differentially la be. cysteine-containing peptides as a means 
to obtam relative abundance data for proteins found in two distinct sampies in a 
s,ng,e anaiysis. Other approaches based on differentia, stable isotope labeling 
have been devised (Munchbach et at. 2000). The .CAT method is unique in that 
,« specifically enriches for peptides containing the relatively rare amino acid 
cys erne, thereby simplifying complex protein mixtures for subsequent MS 
analys,, The relative abundance of proteins can then be determined by 
momonng the ratios of pairwise seta of selected peptide species which are oflj 
by 8 amu. While representing a major advance, the ICAT approach is based on a 
soph,s.,ca.ed proprietary chemistry that analyzes ,e,a,.ve«y rare cysteine- 
containing peptides. 

Here a complementary protein identification and quantitation strategy is 
descnbed, which is tanked Mass Coded Abundance lagging (MCAT). based on 
the drfferenual post-experiment labeling of tryptic peptides with the .ysine 
guan,da,,on agent O-methylisourea followed by high throughput capillary liquid 
chromatography electrospray tandem mass spectrometry (LC-MS/MS). MCAT 
permrts fade de novo sequencing of proteins present a, pico- to femtomole 
levels ,n complex biological mixtures and provides for ro bus, determination of me 
relate abundance of proteins in various ce,l states in a systematic, reproducible 
an d stra, 9 htforwar d mannerThe d evelopment a nda B p,ILcations.of.a.systematio 
prote.n expression profiling strategy based on the MCAT approach outlined here 
should serve as a powerful means for characterizing the physiologica. 
development or disease state of cells or organisms at, he prbteome level. 

Results 

De novo Peptide Sequencing using MCAT 

T encin9 me,hod re,ies on ,he seiec,ive and <-"«^ («. 

complete) medication of the e-amine of C-tarmina, .ysine residues of tryptic 
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shown that it can be used to sequence multiple individual peptides from complex 
mixtures in a single high-throughput electrospray LC-MS/MS analysis. 

The MCAT de novo sequencing approach is based on two principles. First, a 
short sequence of contiguous amino acid sequence from a peptide (5-10 
residues) usually contains sufficient information to identify a corresponding 
unique protein. Second, peptides alternatively unmodified and modified with O- 
methylisourea differ by the mass differential encoded by the MCAT reagent (42 
amu). This allows the identities of the informative y-ion peaks to be readily 
delineated by comparing pair-wise sets of MS/MS spectra, allowing for 
systematic sequence determination. The MCAT labeling procedure is simple, 
economic and easy to perform with complex protein mixtures. 

The steps of the MCAT peptide sequencing procedure are as follows: (1) A 
protein mixture, which can be a purified polypeptide or protein complex, a cell 
fraction, or a crude cell extract, is first digested enzymatically with trypsin; (2) Half 
of the digest is derivatized to completion following incubation with an excess O- 
methylisourea; (3) The digests are desalted by C18 solid phase extraction and 
combined; (4) The pooled peptide mixture is fractionated by reverse phase HPLC 
and analyzed by automated ESI MS/MS. The mass spectrometer is operated in 
an automated dual mode whereby successive scans alternatively record a) the 
m/z of modified/unmodified peptide pairs as they elute from the column and b) 
the MS/MS fragmentation pattern of each peptide that has undergone collision- 
induced dissociation (CID); (5) Following MS analysis, the data are processed to 
obtain the amino acid sequence identities of the components of the protein 
mixture. The process is illustrated schematically in Figure 1B. 

Inspection of pair-wise peptide spectra indicates that most ion peaks, notably the 
b-ion and y-ion series, are retained upon modification (Table 1). Since the C- 
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recorded m/z values for b-ions and chemical noise remain unchanged. Therefore 
comparison of MS/MS spectra for each unmodified/modified peptide pair allows 
ready determination of the y-ion peaks. With high quality spectra, discrimination 
of a well-defined and continuous y-ions series allows the amino acid sequence of 
a peptide to be readily deduced. This simplifies the spectral interpretation 
process, allowing for systematic sequence determination by assigning amino acid 
masses that correspond to y-ion peak distances using a reference table of 
moncsotopic amino acid masses. If required, a delta mass corresponding to a 
poss,ble post-translational modification (e.g. + 80.0 amu for phosphorylation on 
senne, threonine or tyrosine residues, or neutral loss (eg. water or ammonia) can 
be incorporated into this table. 

In a systematic series of studies using a crude yeast cell extract (Table 1) it is 
established that MCAT provides an effective method for sequencing multiple 
peptides analyzed by LC-MS/MS. Firs,, the ionization, charge and fragmentation 
properties of peptides were not greatly affected by the chemical derivation 
procedure. Peptides generally have one of three different charge states (-1 +2 
or -3), each of which results in a unique spectrum for the same peptide The 
spectra of numerous unmodified and modified peptide forms showed similar 
^formation content and could be correctly interpreted using database search 
algorit hms with similar efficiency. Seco nd, the modification of lysinfccontaining 
pephdes-ocourred-in a robust, unbiased and reproducible" manner. TNrd the 
mass tag (42 amu) added to the treated peptides was easily resolvable by MS 

pept.de modifications. Even for a charge state of +3. the delta mass is 14 units 
well within the resolution of a mass spectrometer. Fifth, the process simplified the 
spectral interpretation process so that the area of combinatorial sequence space 
to be searched was easily within the limits of modem computing technology 
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using a computer database search algorithm. The SEQUEST algorithm (and 
similar algorithms) can detect MCAT modified lysine residues unequivocally 
because modification of a C-terminal lysine following trypsin digestion alters the 
m/z of y-series ions but not b-series ions relative to the unmodified peptide. 

Although carried out manually here, the MCAT sequencing process may be 
formalized to facilitate automation. First, the mass of the tag (or a factor of it 
resulting from multiple charges) is added to each peak observed in the 
unmodified spectrum (above some threshold). The spectrum of the modified 
peptide is searched for peaks corresponding to these 'mass-tagged 1 peaks, any 
such peaks being candidate y-ions. Peaks appearing in both spectra are likely to 
represent b-ions or other ion products and are excluded from the initial analysis. 
Next, the mass differences between all candidate y-ions are calculated. Mass 
differences matching the known masses of single or double amino acids are 
noted and attempts are made to extend the sequence from this starting point in 
both directions (i.e. higher and lower m/z) using known single or double amino 
acid masses. The putative sequences can be ranked using a score incorporating 
factors such as unbroken peak series and correlation of observed peaks with 
theoretical peaks. Moreover, for each putative y-ion series, the remaining peaks 
(L. those conserved in the unmodified and modified spectra) are candidate b-ions 
and therefore can be used to impose further statistical limits on the y-ion 
designations. In other words, for any identified y-ion sequence ACDEFG, the 
corresponding sequence GFEDCA should be observed, and the extent of the 
presence or absence of the corresponding peaks can be factored into the overall 
score. 

Our results are typical of peptide MS/MS experiments in that incomplete y-ion 
series were generally observed. For high mass y-ions (yn, yn-1), this may occur 
because of charge repulsion; for low mass y-ions (y2, y3), because ion trap 
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Table 2 shows that MCAT reagent selectively modifies all lysine-terminated 
trypt.c peptides present in the mixture in a quantitative and robust manner In 
order to show that modification by the MCAT reagent is specific and that peptides 
so modified are recognizable by spectral identification algorithms, LC-MS/MS on 
a control yeast extract and a yeast lysate that had been treated with O- 
methylisourea was performed. The acquired MS/MS spectra were typically of 
high quality, with distinct b-series ion patterns the same for modified and 
unmodified spectra and the y-series offset by 42 Da, confirming that a C-terminal 
lysine had been modified (Fig. 2). Moreover, the SEQUEST scores for both 
modified and unmodified peptides were comparable and typical of high fidelity 
.dent.fications. Importantly, in no case was an unmodified peptide detected in the 
treated sample (i.e. yielding high SEQUEST scores). The corollary was also true 
w.th no peptides being significantly scored as being modified in an untreated 
sample (Table 2). 

Comprehensive LC-MS/MS analysis of an untreated and an O-methylisourea 
mod.fied yeast cell lysate yielded significant SEQUEST scores for 291 peptides 
For peptides treated with O-methylisourea, the rate of modification of non-lysine 
rescues, such as arginine or alanine, by O-methylisourea was negligible (data 
not shown), asreported by others (Ki mmel, 1967; Ha leetal., 2000; Beardsley et 
at. 2000). Greater than 95% of SEQUEST-validated peptides containing lysine 
residues were classified as modified at lysine. In contrast, less than 3% of 
untreated peptides were scored as modified by SEQUEST. the same rate of 
false-positive scoring observed for arginine-containing peptides. These false- 
pos,t,ves may result from poor quality spectra, or from acetylation or 
tnmethylation of amino acids that generate a gain in mass (monoisotopic) of 
42.0106 Da or 42.0471 Da respectively. Such false positives can be easily 
eliminated upon inspection of MS/MS spectra because the y-ions series do not 
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isobaric nature of certain amino acids (e.g. leucine and isoluecine). The MCAT 
approach is limited to peptides that terminate with a lysine residue. Tryptic 
fragments ending with arginine resdues are not modified and, therefore, cannot 
be sequenced by this approach. If necessary, endoproteinase LysC can be used 
instead of trypsin to generate peptides ending exclusively in lysine residues 
(apart from peptides derived from the C-terminus). Finally, it should be noted that 
incomplete trypsin or LysC digestion can potentially complicate the MCAT 
sequencing process by causing a mass shift in a subset of b-ions. However, the 
presence of modified internal lysine residues can be readily detected a priori by 
searching for parent ion mass shifts of multiples of 42 amu (adjusted for the 
charge on the ion). 

Relative Protein Abundance Determination Using MCAT 

The MCAT approach allows the relative abundance of proteins to be compared in 
two different samples following differential modification of peptides from one of 
the samples with O-methylisourea. By combining the peptides after treatment, the 
relative abundance of different protein species present in each sample can be 
estimated by measuring the signal intensities of the peptide pairs in a full scan 
MS analysis. The basic MCAT approach for measuring protein abundance is 
outlined in Figure 1C. 

In general, a first test sample and a second test sample may be an experimental 
sample (e.g. a sample exposed to a test compound of interest) and a control 
sample (not exposed to the test compound), respectively. Both samples are 
preferably enzymatically digested, for example in trypsin, and then one of the 
samples is treated (derivated) with a reagent to create a mass differential. This 
reagent may be called a mass differential reagent and is preferably a lysine 
guanidination compound. It may be, for example, o-methylisourea or any 
compound suitable for MCAT, that creates amino acids terminating in lysine or a 
homoarginine ending group or variant (memetic) thereof. The peptide of each 
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curve in a single ion intensity profile. Preferably, the peptide profile and relative 
abundance in the first and second sample is carried out in both orientations. 

MCAT protein quantitation is based on two principles: First, pairs of peptides 
alternatively unmodified and modified with O-methy.isourea can be discriminated 
dunng a single MS run, thereby serving as mutual internal references for 
accurate relative quantitation. In MS, the ratios between the recorded signal 
intensities of the lower and upper mass components of these ion pairs provide a 
d.rect measure of the relative abundance of the two forms of a peptide and by 
inference, the corresponding proteins in the original cell pools: Second 'the 
-dentrty of the peptides can be obtained by performing MS/MS during the same 
analysis. 



The steps of the MCAT peptide quantitation procedure are as follows- (1) Two 
protein m.xtures to be compared are obtained following different experimental 
treatment of a cell or tissue and are digested enzymatically with trypsin; (2) One 
d,gest is derivatized with O-methylisourea; (3) The peptides are desalted by C18 
sobd phase extraction, combined, and the isolated peptides are separated and 
analyzed by automated multistage LC-MS/MS. The mass spectrometer is 
operated in a dual mode where two alternative scans cycle repeatedly. First a full 
MS scan monitors the signal intensity <rf peptides eluting from the capillary 
column. Second, peptide sequence information is generated by selecting peptide 
.ons for CID fragmentation in MS/MS mode. Sequence identification can be done 
us,ng the de novo approach described above or using a protein database search 
algorithm. (4) Peptides are quantified by comparing the relative signal intensities 
of pa,rs of peptide ions with identical sequence that differ in mass due to lysine 
guamdination. In practice, an ion intensity profile is reconstructed for each 
sequenced peptide using the MS data and the relative abundance of modified 
and unmodified peptides calculated by integrating the area under the curve The 
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The MCAT approach serves as an effective method for determining relative 
abundance of proteins by LC-MS/MS since: (1) O-methylisourea derivatizes all 
lysine-containing peptides present in the mixture in a quantitative manner; (2) the 
agent adds a mass tag to the treated peptide that is easily resolvable by the 
mass spectrometer and that does not overlap with common adducts or peptide 
modifications; (3) the modification preserves the charge and ionization properties 
of peptides such that the efficiency of ionization and signal intensity are 
equivalent; and (4) the modified peptides generally co-elute during standard 
reverse phase chromatographic separation. 

To illustrate the process, the relative abundance determination of the peptide 
LPWFDGMLEADEAYFK (SEQ ID NO:15) from two replicate yeast whole cell 
extract experiments is shown in Figure 3. Base peak chromatograms show many 
peptides eluting over a 60min run, while selected ion tracings for the predicted 
doubly-charged unmodified and modified forms of the peptide show both eluting 
at 35-36min (Fig. 3A). A single full scan of an ion trap mass spectrometer 
operated in MS mode is shown in Figure 3B. Two prominent ion species are 
discernable and indicated with respective m/z values 21 m/z units apart (Fig. 3B). 
The fact that the ions co-elute, have a detected mass difference of 21 m/z units, 
and have identical sequences (data not shown) identifies them as a pair of doubly 
charged sister peptides. Over the course of the 60 minute elution gradient, more 
than 2,000 MS scans were automatically acquired. Figure 3C shows 
reconstructed ion chromatograms for each of the peptide species. The relative 
quantities were determined by integrating the curves contouring the respective 
eluting peaks. The ratio (unmodified:modified) was determined as 0.88 (Table 2). 
The peaks in the reconstructed ion chromatograms appear serrated because the 
MS system alternates between MS and MS/MS modes in order to both measure 
ion intensity as well as generate a mass spectrum of selected peptide ions for the 
purpose of protein identification. 
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select peptides throughout an entire chromatographic run typically showed 
isolated peaks with the unmodified form co-eluting. or eluting slightly earlier, than 
the modified form (Fig. 3A and C). For nearly all peptides examined, the ratio of 
unmodified to modified signal was close to the expected 1:1. The range of signal 
intensities were generally within two-fold of the unmodified form and the 
percentage error (the difference between the observed and expected 
abundances) ranged from 1 to 62% (Table 2). Some exceptions were evident and 
excluded from the analysis. These included peptides that could be positively 
identified but whose signal is very weak, and peptides containing arginines that 
were modified in addition to lysine at low frequency. Another category of ion 
found unsuitable for quantitation were singly-charged ions. It is unclear why this 
is the case but the signal from singly-charged ions is typically lower than that for 
doubly- or triply-charged ions, possibly rendering them less likely surpass the 
intensity threshold required for accurate quantitation. 

Figure 4 shows variation in the measured relative abundance for two peptides 
from the same parent protein (and therefore are present in equimolar 
concentrations) in three replicate experiments. Importantly, multiple peptides 
independently analyzed for several proteins gave similar linear responses. 
Experiment-to-experiment variation for these peptides is within 25% and variation 
within a single experiment for peptides derived from the same protein is within 
20%: The variation from ideal 1:1 ratio is nofthe result of reduceT ionization 
efficiency or MS signal of the modified peptides relative to their unmodified forms 
because the effect was consistently observed in subsequent experiments 
independently of which sample was chosen for modification. More likely, it results 
from modest variations in peptide recovery during sample workup. 

In order to correct for any possible systemic labeling errors, for instance 
preferential labeling by O-methylisourea of one sample, MCAT quantitation can 
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sample A versus modified sample A; IV) unmodified sample B versus modified 
sample B. The ratios of unmodified to modified peptide signals obtained in 
experiments III and IV can be used to systematically normalize and control for 
variations in the data obtained in experiments I and II, respectively. In practice, 
the MCAT analysis can be simplified into a two-tiered reciprocal experiment set, I 
and II, which should independently confirm any significant quantitative 
observations obtained in a sample comparison. 

To confirm the quantitative nature of the MCAT approach, mixtures of modified 
and unmodified peptides derived from a common crude yeast cell extract were 
prepared at various ratios and analyzed by a 30 minute LC-MS/MS analysis. The 
MS/MS spectra acquired were used to search a non-redundant genome 
database using the SEQUEST algorithm (Eng et al., 1994) to identify the proteins 
present in mixtures. The relative ratios of 5 peptide sister pairs was quantified as 
described above (Fig. 4B). This analysis shows the relative abundance of 
proteins can be accurately determined (i.e. exhibits a linear response) over a >30 
fold dilution series. Beyond this range, the weaker signal of the two abundances 
was indistinguishable from background noise in these experiments. 

It should be emphasized that the data were acquired for polypeptides present at 
a pico- to femtomole level in a highly complex protein mixture. The loading 
capacity of capillary reverse phase columns for complex peptide mixtures 
imposes a strict limit on the detection of low abundance proteins by LC-MS/MS. 
With a purified protein, most current MS systems generally exhibit a practical 
dynamic range of roughly three orders of magnitude based on maximal signal to 
noise ratios that can be acquired (using a purified or low complexity protein 
preparation). However, sophisticated chromatographic separation techniques can 
be coupled to fractionate complex peptide mixtures prior to MS in order to 
substantially improve the detection limits of MS protein analysis (Link et al., 1999; 
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An experimental approach for systematically sequencing and quantifying proteins 
isolated from complex biological mixtures using basic chemistry and mass 
spectrometry techniques is described and validated. De novo sequencing 
expands the range of organisms that can be analyzed and removes the reliance 
on DNA sequence databases that may be incomplete, erroneous, or that fail to 
account for complexities introduced by alternative splicing, protein modifications, 
or protein polymorphism. The quantitative capabilities of the method also 
overcome a significant limitation of current proteomics technologies, whereby the 
determination of protein abundance on a large-scale is generally low throughput, 
expensive, and tedious, for instance, radiolabelling of proteins before analysis by 
two-dimensional gel electrophoresis and quantitation following isolation of 
individual spots (that may contain one or more polypeptides). 

The ICAT method reported by Aebersold and coworkers (Gygi et al., 1999) may 
significantly improve throughput and reduce sample complexity by enriching for 
proteins containing the underrepresented amino acid cysteine. These features 
are useful for sampling a mixture whose proteome complexity could overwhelm 
the ability of current LC-MS technology to resolve it. The MCAT strategy 
described here is not limited to any particular affinity chemistry and in principle 
can be coupled to analogous affinity-based enrichment steps. For this reason, 
MCAT can potentially be used to identify and quantify all the proteins presentin a 
biological sampler In combinatioh with powerful multi-dimensional LC protein 
separation techniques, such as that described by Yates and coworkers (Link et 
al., 1999; Washburn etal, 2001), considerable depth in proteome coverage may 
be achieved. Quantitative data describing patterns of peptide or protein 
expression for many hundreds or thousands of proteins can be used to identify or 
classify protein 'profiles' in a similar manner to that routinely used for gene 
expression data. The combined MCAT approach can therefore be used for 
identifying, classifying and characterizing functions of known and unknown gene 
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First, the approach is simple and effective. It builds on established MS techniques 
and principles that are flexible and can easily be adjusted for large-scale projects, 
including efforts to generate peptide or protein profiles describing the effects of 
environment, mutation, disease or experimental interventions such as drug 
treatment. Significant patterns of expression can be identified with appropriate 
software and data mining algorithms. 

Variations of the MCAT approach can easily be devised, including strategies to 
address other quantitative aspects of protein expression, those searching for 
post-translational modifications, or those screening for mutant proteins. It is likely 
that the number of unique peptide species per organism will be multiplied 
significantly by the presence of post-translational modifications compared to 
genome predictions. Because the mass of many common important modifying 
groups are known, and because their preferences for particular amino acids are 
often known, the database can be searched for ions predicted to result from 
peptides with specific modifications. 

Finally, the addition of a dynamic component to the molecular descriptions of 
protein activities is likely to prove critical to our understanding of the biochemical 
circuitry within cells. Consequently, the development of robust analytical 
methods, such as the MCAT approach described here, that allow for efficient 
identification and quantitation of large numbers of proteins from complex mixtures 
can be expected to have a major impact. 
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Experimental protocols 

Materials. Media, standard-grade and HPLC-grade laboratory chemicals were 
obtained from Fischer Scientific (Fair Lawn, NJ). O-methylisourea (S- 
methylisothiourea hemisulfate salt) was from Sigma-Alderich (St. Louis, MO). 
Poroszyme immobilized trypsin was from Applied Biosystems (Framingham, MA). 

Preparation of protein extracts. The protease-deficient S. cerevisiae yeast 
strain BJ5460 was grown to late-log phase (OD -3) at 30"C and protein whole 
cell extracts prepared as follows: Cells were harvested, frozen, and mechanically 
lyzed by grinding in the presence of dry ice. The cells were thawed in lysis buffer 
(8M urea. 1 mM CaCfe, 100 mM Tris-HCL, pH8.5). Insoluble debris was pelleted 
by a high-speed (20 K x g) spin and the supernatant diluted to 2M urea using 
digestion buffer (100 mM Ammmonium bicarbonate, pH8.5, 1 mM CaCI2. A 
bacterial whole cell extract was similarly prepared using the E. coli DH5a strain. 
Human nuclear extracts were prepared using a commercial kit (Pierce), and 
diluted into digestion buffer. 



Tryptic Digestion and Peptide Derivatization. Porozyme immobilized trypsin 
beads were added to an aliquot of each protein extract at a 1:500 protein ratio 
and the digests incubated at 30°C for two days with tumbling. The extracts were 
aliquoted into two microtubes. Solid O-methylisourea was added to one of the 
tubes to achieve a final concentration of 1M. Base (NaOH) was added to 0 5N to 
adjust the pH to >10. The reaction was incubated at 37°C overnight. The peptide 
mixtures were extracted by solid-phase extraction using SPEC-PLUS PTC 18 
cartridges (Ansys Diagnostics, Lake Forest, CA) according to the manufacturer's 



41 

SUBSTITUTE SHEET (RULE 26) 



WO 02/097703 PCT/C AO 2/00801 

MCAT peptide sequencing. Each sample was subjected to microcapillary LC- 
MS/MS analysis with modifications to the general method described by Link and 
coworkers (1999). A quaternary Surveyor HPLC pump (ThermoFinnigan Canada) 
was directly coupled to a Finnigan LCQ-DECA ion trap mass spectrometer 
equipped with a custom microLC electrospray ionization source. A fused-silica 
microcapillary column (100 \im i.d. x 365 jam i.d.) was pulled with a Model P-2000 
laser puller (Sutter Instrument Co., Novato, CA) as described. The microcolumn 
was packed with 10 cm of 5 \vn Ci 8 reverse-phase material (Zorbax XDB-C18, 
Hewlett-Packard). Approximately 100 \ig of the unmodified fraction and 100 \ig of 
the derivatized peptide fraction were combined and loaded onto a single 
microcolumn for sequence analysis. After loading, the column was placed in-line 
with the ion source system setup as described (Link et al, 1999). A fully 
automated 30 min 100% buffer A (5% ACN, 0.1% formic acid) to 80% solvent B 
(95% ACN, 0.1% formic acid) binary gradient was run at a flow rate of -0.3 
ul/min. Eluted peptides were analyzed by automated MS/MS as described by 
Link and coworkers (1999) except that a full scan range of 400-1600 m/z was 
used. 

SEQUEST analysis. The SEQUEST algorithm (Eng et a!., 1994) was run on 
each dat set against sequence databases obtained from the National Center for 
Biotechnology Information (Bethesda, MD). Positive sequence identification was 
based on several criteria (XCorr and DCn score, and the presence of tryptic 
termini) described at http, and all identifications were confirmed manually. 

MCAT protein quantitation. Pairs of samples to be compared were subjected to 
automated uLC-MS/MS analysis with modifications to the general method 
described above. Approximately 200 [ig of the unmodified fraction and 200 |ig of 
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-C-MS-MS/MS techniques as described above. There was a consistent slight temporal 
fference in the elution of unmodified/modified peptide pairs, with the unmodified light 
•alog eluting slightly before the heavy form. Selected ion traces for each peptide pair 
sre quantified using the ADDXPRESS program by which the peak area of each eluting 
iptide was reconstructed and used in the ratio calculation. 



Me 1. De novo peptide sequencing from complex mixtures using MCAT 
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^- ( 

3 and b* refer to unmodified and modified b-ion series respectively 
i and y* refer to unmodified and modified y-ion series respectively 

A indicates a match between expected and observed m/z values (tolerance of 2.0 m/z units) 
*b, Difference between observed b and b* m/z values 
ly, Difference between observed y and y* m/z values 

*(y,y+l), Difference in observed m/z between successive y series ions, adjusted for charge state of ion 
Predicted AA, Amino acid residue predicted using A(y,y+1) 

indicates a match between M CAT- predicted and SEQUEST-predicted amino acid. 



3/>/e 2. Identification and quantitation of peptides from a yeast whole cell digest 
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BR118W 


SVEMHHEQLEQGVPGDNV 
GFNVX fSEO. ID. NO: 19) 


2550.8/ 
2592.8 


2 


1276.4/ 
1297.4 


2.2433/ 
2.5321 


gA 


X i 


X 


gA 


1.00 


0.76 


24 ±4 




TLLEAIDAIEQPSRPTDKPL 
RLPLQDVYK* (SEQ. ID. 
NO:20) 


3320.8/ 
3404.8 


3 


1107.9/ 
1135.9 


3.3888/ 
3.3370 


gA 


X | 


X j 


gA 


1.00 


0.63 


37 ±5 




VETGVIKPGMWTFAPAGV 
TTEVX# (SEQ. ID. 
NO:21) J 


2430.9/ 
2472.9 


2 


1216.4/ 
1237.4 


2.5458/ 
2.1831 


■A 


X 


X 


gA 


1.00 


0.38 


62 ± 12 


CR012W 


ALEN FTRPFLAILGGAK 
(SEO. ID. NO:22) 


1768.1/ 
1810.1 


2 


885.0/ 
906.0 


1.7773/ 
1.4083 


»A 


X 


X 


gA 


1.00 


0.57 


43 ±5 


DR155C 


HWFGEWDGYDIVK 
(SEO. ID. NO:23) 


1675.9/ 
1717.9 


2 


838.9/ 
859.9 


3.7988/ 
3.6211 


*A 


X 


X 


gA 


1.00 


0.71 


29 ± 5 


DR487C 


HGIPUSIEELAQYLK 
(SEO. ID. NO:24) 


1824.2/ 
1866.2 


2 


913.1/ 
934.1 


2.1238/ 
1.6387 


■A 


X 


X 


gA 


1.00 


0.86 


14 ± 1 


GR063C 


LPAEWELLPHYKPR 
(SEO. ID. NO:25) 


1761.1/ 
1803.1 


2 


881.5/ 
902.5 


2.0444/ 
1.9739 


gA 


X 


X 


gA 


1.00 


0.66 


34 ±6 


GR192C 


INDAFGIEEGLMTTVHSLT 
ATQK 

(SEO. ID. NO:26) 


2476.8/ 
2518.8 


2 


1239.4/ 
1260.4 


2.9164/ 
4.1100 


gA 


X 


X 


gA 


1.00 


0.52 


48 ±28 




VIN D AFG IEEGLMTTVH 5L 
TATQK 

(SEO. ID. NO:27) 


2575.9/ 
2617.9 


2 


1288.9/ 
1309.9 


3.1456/ 
3.3717 


gA 


X 


X 


gA 


1.00 


0.44 


56 ± 17 




VPTVDVSWDLTVK 
(SEO- ID. NO:28) 


1512.7/ 
1554.7 


2 


757.3/ 
778.3 


3.2279/ 
3.1548 


gA 


X 


X 


gA 


1.00 


1.29 


29± 11 


GR214W 


NVQVHQEPYVFNARPDGV 
HVINVGK (SEQ. ID. NO: 
29) 


2817.2/ 
2859.2 


3 


940.0/ 
954.0 


1.8494/ 
2.2204 


.A 


X 


X 


gA 


1.00 


0.61 


39 ±10 


GR254W 


AQYNBQGWDHLSLLPTF 
GAK (SEQ. ID. NO:30) 


2388.7/ 
2430.7 


2 


1195.3/ 
1216.3 


2.4748/ 
3.0844 


*A 


X 


X 


gA 


1.00 


0.81 


19 ±2 1 




YPIVSIEDPFAEDDWEAW 
SHFFK 

(SEO. ID. NO:31) 


2829.1/ 
2871.1 


3 


944.0/ 
958.0 


3.1108/ 
3.2183 


gA 


X 


X 


gA 


1.00 


0.61 


39 ±9 


HR174W 


WLTGVELADMYHSLMK 
(SEO. ID. NO:32 


1894.2/ 
1936.2 


2 


948.1/ 
969.1 


4.0552/ 
3.8246 


gA 


X 


X 


gA 


1.00 


0.77 


23 ± 3 


JR105C 


TVTFTHGVEPTVWSSK 
(SEQ. ID. NO:33) 


1800.1/ 
1842.1 


2 


901.0/ 
922.0 


1.5600/ 
1.8810 


«* 


X 


X 


gA 


1.00 


0.75 


25 ±4 


KL060C 


SPIILQTSNGGAAYFAGK 
(SEO. ID. NO:34) 


1795.0/ 
1837.0 


2 


898.5/ 
919.5 


3.6709/ 
4.2032 


gA 


X 


X 


gA 


1.00 


0.73 


27 ±5 




TGVTVGEDVHNLFTYAK 
(SEQ. ID. NO:35) 


1863.1/ 
1905.1 


2 


932.5/ 
953,5 


3.2735/ 
2.6813 


A 


X 


X 


gA 


1.00 


0.75 


25 ±4 


LR044C 


KUDLTQFPAFVTPMGK# 
(SEO. ID. NO:36) 


1906.3/ 
1948.3 


2 


954.1/ 
975.1 


3.5845/ 
3.9361 


-A 


X 


X 


gA 


1:00 


0.83 


17 ±2 


LR058C 


EVLYDLEN PIN FSVFPQHQ 
GGPHNHTIAALATALK 
(SEQ. ID. NO:37) 


3772.2/ 
3814.2 


3 


1258.4/ 
1272.4 


1.8356/ 
2.5693 


gA 


X 


X 


gA 


1.00 


0.73 


27 ±6 



a. Molecular mass of unmodified/modified peptides ions. 



b. Mass-to-charae ratio of unmodified/modified oeotides. 
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urther discussion of the figures related to MCAT 



J The MCAT approach for peptide sequencing and relative protein 
bundance determination. 

ee Figure 1 . (A) The guanidination reaction is specific for the side chains of lysine, 
hich is selectively converted to homoarginine. (B) For sequencing using MCAT, 
-otein mixtures are first digested with trypsin, which generates peptides suitable for MS 
nalysis that terminate with lysine or arginine residues. Half of the sample is treated 
ith the MCAT reagent O-methylisourea. Peptides ending in lysine are modified, which 
Jds 42 amu to the mass of the peptide but does not alter the properties of the peptide 
Jring LC-MS analysis. The peptides mixtures are combined at a 1:1 ratio, separated 
/ reverse phase LC and introduced online into a MS instrument using electrospray 
nization. Following tandem MS analysis, peptide sequence is determined by 
sparing MS/MS spectra of unmodified and modified peptides. The fragmentation 
attern of both sister peptide pairs are similar except for the shifted y-ion series, which 
an be deconvoluted to reveal the amino acid sequence of the peptide. (C) For relative 
Dundance measurements, samples representing different cell states are alternatively 
.odifiedor unmodified with MCAT. Full MS spectra are recorded for sister peptide 
Decies and their relative abundance determined by measuring the respective trace 
tensities on reconstructed single ion chromatograms. 
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(2) MCAT enables identification and quantitation of complex protein 
mixtures. 

See Figure 2. (A) Ion chromatograms recorded for the base peak (top), an 
unmodified peptide ion [LPWFDGMLEADEAYFK+2H]* 2 (middle) and its 
corresponding 0-methylisourea(MCAT)-modified form (bottom). When mixtures 
of untreated and MCAT-treated protein digests are resolved by reverse phase 
LC, the modified peptides elute with a minor delay compared to the respective 
unmodified forms (35.9 vs. 35.7 min respectively in this example). (B) Depending 
on charge and the number lysine residues, the m/z signals observed for pairs of 
unmodified or modified peptide ions during MS are offset by 42, 21 or 14 m/z 
units (for plus 1, 2 or 3 ions respectively). In this example, the peak signals 
recorded for the unmodified (967.07 m/z) and modified (988.08 m/z) forms of the 
peptide are offset by 21 m/z units, indicating a +2 charge. The peptide ions are 
then independently selected and automatically fragmented by MS/MS. 
Comparison of the y-ion series allows the amino acid sequence to be determined. 
(C) The relative abundance of individual peptides can be determined by 
reconstructing the chromatograms for the unmodified and modified forms of the 
peptide ions and calculating the ratio of signal intensities using area under curve 
integration. 

(3) De novo sequencing of a yeast peptide and a human peptide using 
MCAT approach. 

See Figures 3A and 3B. (A) The peptide WDLVEHVAK (SEQ ID 
NO:38)analyzed by MCAT LC-MS/MS in a digest of yeast whole cell extract. A 
representative MS/MS spectrum of the unmodified peptide (top) and the 
corresponding spectrum for the modified form (below) are shown. Because the 
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for b- and y-series ions of the unmodified and modified peptides are given (right, 
w„h those observed in the experiment underlined. The amino acid order is 
resolved by measuring the mass difference between successive y-ion peaks (B, 
The pepbde VAPEEHPVLLTEAPLNPK (SE Q ,D NO :3 9 ,was identtfied in a digei 
ETUUVP^Eft' Ml * *' S P * pl " te a stre4ch of ten amino acids (A- 

h^tT T ' °" te iden,Bied by mappi " 9 y " ions to — ">-* ^ed 
by 42 mft unrts ,n the modified spectrum (bottom, relative to the unmodified 

spectrum (top,. The dominant peak at 892.9 in the unmodified specttom is 
approximates 21 m/z units from an dominant unassigned peak a, 914,4 in ,he 
modified spectrum. These peaks probably represent doubty-charged y16 ions " 
that term,nate in with proline, an amino acid commonly observed to fbnn 
dominant peaks during CD. The other major peak in bo* spectra (1292.6 and 
1334.5 ,n the upper and lower panels respecfively, is a singly-charged y,2 ion 
ha, a.so terminates wfih proline. Therefore, an addttiona. advantage of the MCAT 
techmque ,s the resolution of such ambiguous peaks through chame 
determination. In the case of both yeas, and human peptides, me identical 
molecular masses o, leucine and isoleucine preven, their resolution by MS 



(4) 



The MCAT method is reproducible and quantitative. 



See F,gures 4A_and 4B. (A, A yeast whole cell was digested wHh frypsin in Ihree 
rep,,cateexpe- fl mehta(A,B;^ 

one o, wh,ch was tteated w«h O-methylisourea. Each pair of mixte res was then 

* f 1:1 raU ° a " d protei " «—» determined by the MCAT LC- 
MS/MS. The retafive abundance ratios (expressed a. the ratio o, modified to 
unmodrfied peptide signal, o, a subset of posittvely-idenfified peptides is given for 
each ana,ys,s. ,B, Untreated and MCAT-tabeled yeas, protein ttypflc digLs 
were combined in varying proportions ranging from from 16:1 (modified to 
unmodffied, to 1:16 effecfive concentrations. The measured relative abundance 
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Peptide Profiling 

Below examples are shown of the utility of peptide profiling as a means to 
characterize and classify diverse human tissues, to characterize subcellular 
fractions of individual tissues, and to illustrate how a database of such peptide 
profiles can serve as a depository of protein expression information that can be 
mined rapidly and accurately for knowledge about the status of an unknown 
sample. This process is robust, sensitive and reproducible. Although the method 
is generally applicable, the following serve to illustrate select uses of the 
approach. 

Example 3: Use of peptide profiles to characterize human tissue 
The invention includes methods of characterizing human tissue. The method 
comprises generating samples suitable for MS analysis and producing a peptide 
profile. The relative abundance of peptides in samples is also preferably 
determined. The peptide profile that is generated is compared to peptide profiles 
in a database or library using common algorithms in order to identify cognate 
proteins, preferably those that are considered important therapeutic targets, as 
well as metabolic enzymes and structural proteins. 

Table 1 shows 40 peptides sequenced and quantified from a human lung tissue 
lysate sample in a single LC-MS analysis that are then used to construct a unique 
peptide profile. The peptides in turn allowed for the identification of cognate 
corresponding proteins present in the sample (a total of 867 proteins were 
unambiguously identified in this analysis). Note that the peptides sequences 
obtained by a generic database search algorithm were both preceded by, and 
terminated with, a K or R residue as a result of cleavage of the input proteins by 
trypsin. The sequence of a total of 1896 peptides were determined in this one 
analysis with high accuracy and sensitivity, demonstrating the ability of the 
approach to generate a detailed profile or fingerprint of protein expression of a 
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K^HGPGLEGGLVGKPAEFTIDTK.G 

K.AAPLSLCALTAVDQSVLLLKPEAK.L ^ ^ "° m 
(SEQ. ID N0 42) 

•CAAQAHEDIIHGSGKT (SEQ. ID NO 43) K^HVSFianvTnrlo'i ^ ^ ^ 

K.AASLGSSQPSRPHVGEAA7^ATK:V K.AHVSFKPTVAQQR.I (SEQ. ID. NO:62) 
(SEQ. ID NO*44) ^ 

K.AASWLTHQGSFHGAFR.S (SEQ ,D N045) KA.'^AGC^r SK* 1 * N ° 63) - 

K.AAVFNHFISDGVKK.T (SEQ ID NO"46) KAKF^C^ ^^ 0 - ' D ' N ° :64) 

K.AAVLWELHKPFTIEDIEVAPPICA K.AKFENLCKL (SEQ. ID. NO:65) 
(SEQ. ID NO*47) 

K.AAVSGLWGKV (SEQ ID NO -4m K-AKPWSFIAGITAPPGR.R (SEQ. ID. NO:66) 

KADIIYPGHGPVIHNAEAK.I (SEQ. ID. NO:50) ^ Q L °° , ^ LVN NAALVIMQPFLEVTK.E 

^(.SKS™, KALHASGAKWAVTR.T (SEQ. ID. NO:71) 

K^Q N ^|EKT, SE Q, D . Na54) EJSSSE^^^ 

K.AFISWCTAlAATNPAVR.T(SEO ID NO'5S> ^T L W ™f^ LSQNQ(3GKLA 

^SSST^^ ' ' S^WLKPNMVrPOHACTQKP 

K.AGHPFMWNEHLGYVLTCPSNLGTGLR G K-ALTGGIAHLFKQ (SEQ. ID. NO:76) 
<SEQ.,D.NO:58> ICALVKPQA 1 KPK M( SEQ.,D.MO : 77, 

A further embodiment of the invention includes using profiles suoh as this to 
compare different tissues or experimental samples. For instance, a comparison 
of the peptide profiles for human pancreatic and heart tissues can be made with a 
stmple 2-dimensional plot that can be extended to •„• different planes as required 
(for W types of tissue, samples, or patients). Comparison of the peptide profiles 
of these samples can be done using standard computational methods (e g 
aggiomeraave Custering). In the case of human pancreatic tissue, the analysis 
showed that although several proteins are shared between the tissues, many are 
not Therefore, a further embodiment of the invention is the use of peptide 
profiles to characterize tissues and thereby categorize samples 



( 
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profiling, small molecule metabolite profiling; these methods preferably involve 
tagging the compounds of interest and performing LC-MCAT to generate a lipid 
profile, phosphoprotein profile, small molecule metabolite profile. The methods 
can provide identity and relative abundance information by readily adapting the 
methods described herein with peptides.). 

Table 2 shows some of the corresponding proteins (of the 867 unique proteins 
identified in this analysis) identified by searching the SwissProt Protein database 
using the identified peptide sequences (http://www.expasy.ch/sprot/). 

Table 2. Proteins identified using peptides isolated from human lung tissue. 

P47915 60s ribosoma! protein I29. 5/2000 [MASS=17456] 

P48025 tyrosine-protein kinase syk (ec 2.7.1.112) (spleen tyrosine kinase). 11/1997 
P48147 prolyl endopeptidase (ec 3.4.21.26) (post-proline cleaving enzyme) (pe). 10/1 
P48444 coatomer delta subunrt (delta-coat protein) (delta-cop) (archain). 11/1997 [M 
P48634 large proline-rich protein bat2 (hla-b-associated transcript 2). 2/1996 [MASS 
P48735 isocitrate dehydrogenase (nadp], mitochondrial precursor (ec 1.1.1.42) (oxalo 
P49023 paxillin. 7/1998 [MASS=60937] 

P49137 map kinase-activated protein kinase 2 (ec 2.7.1.-) (mapk-activated protein ki 
P49182 heparin cofactor ii precursor (hc-ii) (protease inhibitor leuserpin 2). 11/19 
P49321 nuclear autoantigenic sperm protein (nasp). 7/1998 [MASS=85191] 
P49327 fatty acid synthase (ec 2.3.1.85) [includes: ec 2.3.1.38; ec 2.3.1.39; ec 2.3 
P49407 beta-arrestin 1. 7/1999 [MASS=46969] 

P4941 1 elongation factor tu, mitochondrial precursor (p43). 12/1998 [MASS=49542] 

P49773 hint protein (protein kinase c inhibitor 1) (pkcM). 7/1998 [MASS=13671] 

P50096 inosine-S'-monophosphate dehydrogenase 1 (ec 1.1.1.205) (imp dehydrogenase 1) 

P50552 vasodilator-stimulated phosphoprotein (vasp). 11/1997 [MASS=39830] 

P50748 hypothetical protein kiaa0166. 11/1997 [MASS=250749] 

P50851 cdc4-like protein (fragment). 7/1998 [MASS=213599] 

P51174 acyi-coa dehydrogenase, long-chain specific precursor (ec 1.3.99.13) (lead). 

P51660 estradiol 17 beta-dehydrogenase 4 (ec 1.1.1.62) (17-beta-hsd 4) (17-beta-hydr 

P51790 chloride channel protein 3 (clc-3). 7/1998 [MASS=84793] 

P51812 ribosomal protein s6 kinase fi alpha 3 (ec 2.7.1.-) (s6kii-alpha 3) (p90-rsk 

P51885 lumican precursor (lum) (keratan sulfate proteoglycan). 7/1998 [MASS=38351] 

P51991 heterogeneous nuclear ribonucleoprotein a3 (hnrnp a3) (fbmp) (d10s102). 7/19 

P52272 heterogeneous nuclear ribonucleoprotein m (hnrnp m). 10/1996 [MASS=77469] 

P52480 pyruvate kinase, m2 isozyme (ec 2.7.1.40). 7/1999 [MASS=57756] 

Cursory examination of this list shows that many interesting and therapeutically 
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A common criticism of current proteomics technologies based on two- 
dimensional polyacrylamide gels is that they are insensitive and only identify high 
abundance metabolic proteins, ie. proteins that are not normally critical 
determents of disease (although these can be important effectors of disease) 
especially since drug development strategies nearly always target low abundance 
proteins important for counteracting a disease phenotype. 

It is clear from the above table that peptide profiling can successfully describe 
many proteins that are considered important therapeutic targets, and not just 
metabolic enzymes and structural proteins. - 



Table 3 shows how proteins from various therapeutically important categories 
were readily identified and quantified in a single analysis. This list was made 
using keywords present in the sequence annotation databases and therefore 
represents the minimum representation of such classes - the vast majority of 
sequenced mammalian proteins await functional annotation. 

By contrast, a recently published study (Proteomics 1,1303-19 A database of 
protein expression in lung cancer.Oh JM, Brichory F, Puravs E, Kuick R, Wood C, 
Rouillard JM, Tra J, Kardia S, Beer D, Hanash S. 2001) where over 1300 2D gels 
were analyz ed from a vari ety of different lung cel l lines and tumors, identified less 
than 200 proteins, the majority of which were metabolic and structural proteins of 
high abundance, and provided no quantitative information. 



Table 3. Peptide profiling identifies therapeutically important proteins. 



Kinases j 

Phosphatases 

Integrins j 

Channel proteins 
Apoptosis proteins 



•» — * _ • . . 



Peptide profiling 



46 
12 



12 



Conventional approach (Oh et al) 
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Example 4. Peptide profiling to characterize diverse human tissues 

One-dimensional LCMS was used to obtain peptide profiles from diverse human 
tissues (Fig. 5). The one-dimensional approach has 2- to 10-fold lower resolution 
compared to two-dimensional approaches but was used in this case to example a 
large number of samples to illustrate the principle. Table 4 shows the number of 
peptides and proteins identified for different human tissues. 

Table 4. The peptide profiling approach can be applied to diverse tissues. 





Proteins 


Peptides 


Brain 


359 


734 


Heart 


114 


231 


Testes 


78 


136 


Liver 


56 


83 


Muscle 


72 


66 


Plasma 


288 


846 


Pancreas 


202 


283 









It is assumed that diverse tissues may express many similar proteins (for 
instance ribosome associated proteins), yet express a subset of unique proteins 
that functionally distinguishes one tissue from another. Similarly, the proteome of 
diseased tissue may be different to healthy tissue. Although this may seem self- 
evident, very few studies have addressed these issues by directly comparing the 
proteomes from different samples. This is largely because of the technical 
impediments mentioned above - conventional techniques generally characterize 
only the most abundant proteins and peptides, and these peptides are least likely 
to differ from tissue to tissue. Figure 6 shows how many proteins were identified 
using MCAT based peptide profiling for a preliminary study of seven human 
tissues. Notably, the peptide and protein profiles of each tissue is distinct. Even 
with this preliminary low resolution analysis, each tissue evokes a different 
signature when subjected to peptide profiling. 
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proteins are highlighted in the large rectangular box, while heart- and brain- 
specific proteins are highlighted in the smaller rectangular boxes. 
Furthermore, the degree of relationship between these tissues can be 
established by comparing such peptide profiles (Fig. 8). Although the principle 
was illustrated here using different human tissues, such analysis can be used to 
detect other proteomic changes, for instance human heart tissue following 
exercise or myocardial infarction, or following administration of drugs. 

^sue' e 4 ' Peptlde pmm,na to char»nt^ n ze subcellular fractions nfa *inrj i~ 

In another embodiment of the invention, peptide profiling can be used to analyze 
the subfractions of a cell, preferably into nuclear, cytoplasmic and membrane 
fractions. This discriminatory power of peptide profiling is illustrated here, where 
the method is used to examine the subfractions of a single clonal cell line. 
Cultured human myoblast cells were processed into nuclear, cytoplasmic and 
membrane fractions and analyzed using the peptide profiling technique (Fig. 9). 
Significantly, over 400 membrane-localized proteins were identified. This class is 
normally very difficult to analyze using conventional proteomics methods yet is of 
particular pharmacologic/therapeutic interest, being the site of receptors and 
channels with critical signaling and transport functions. 

Tables-7=and S^hWhow peptide profiling can be applied to different cellular 
subfractions and used to identify compartment-specific proteins. 

Table 7. Peptide pfoWihg applied to different cell compartments. 



! Peptides \ Proteins 


Cytoplasmic 


2220 


994 


Nuclear 




804 


428 


Membrane 




727 


403 



c 
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Table 8. Peptide profiling identifies compartment-specific proteins 





CytoptasmicMembrane 


Nuclear 


Unique 


805 


249 


262 


Total 


994 


428 


403 


Percent unique 


80 


58 


65 



Example 5: Use of peptide profiles to characterize human cell lines 
In another embodiment of the invention, this invention includes methods of 
characterizing human cell lines. The method comprises generating samples 
suitable for MS analysis and producing a peptide profile. The relative abundance 
of peptides in samples is also preferably determined. The peptide profile that is 
generated is compared to peptide profiles in a database or library using common 
algorithms in order to identify cognate proteins, preferably those that are 
considered important therapeutic targets, as well as metabolic enzymes and 
structural proteins. In a further embodiment, these profiles can comprise a small 
prototype database or library, against which novel samples may be screened. 

A number of peptides from four human cell lines of distinct cellular origin are 
identified by mass spectrometry and linked to their parent proteins. This profile is 
one-dimensional because no addition information about the peptides (e.g. 
quantitative information) is included. Table 6 shows the number of peptides and 
proteins identified for the different human cell lines. 

Table 6. Peptide profiling of different cultured cells 





Proteins 


Peptides 


Myoblasts 


576 


1373 


HeLa 


974 


2067 


NYP17 


192 


290 


Raji-Jurka 


233 


376 



Here, an independent extract of one of the four cell lines is screened and 
demonstrates how this extract can be conclusively shown to be highly similar or 
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using an ion trap mass spectrometer (Deca, Thermoquest, USA) following 
separation of digested peptides using online HPLC. The mass spectrometer was 
programmed to collect primary MS spectra from parent ions, as well as tandem 
mass spectra of daughter ions generated from the first, second and third most 
abundant ions observed in the program window. These spectra were then used 
to search nonredundant genome databases using the SEQUEST algorithm 
(Yates et al., 1995) to identify the peptides and proteins present in the samples. 

The following table shows the top-scoring peptides identified in the analysis of 
one of these cell lines, Juricat: experiment, After statistical filtering, 74, 91 96 
123 peptides were used to identify 55, 62, 49, 59 different proteins in the ' 
respective cell. The peptides for all four cell lines were deposited into a database 
in this case a Microsoft Access file. The protein profiles are graphically 
represented below (5922, 4091, 5644 and 4166 tryptic peptides were observed 
from MCF7, TPA, Jurkat and K566 cells respectively, visual representation)- 
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If these profiles are considered as a small index or database, novel profiles can 
be searched against them using any common correlation test. For instance here 
the correlation is calculated by: 
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Correlation scores, P x , y , for one-dimensional peptide profiles obtained from four 
human cell lines: 





MCF7 


TPA 


Jurkat 


K556 


? I 


MCF7 


T 


0.0105 


0.33596 


0.09 


0.07 


TPA 


0.0105 


1 


0.33596 


0.31714 


0.26733 


Jurkat 


0.33596 


0.33595 ~ 


1 


~O09 


.8644 


K556 


0.09 


0.31714 


0.09 


1 


.0.09 



_ This preliminary analysis suggests that-the peptide profiles obtained from Jurkat 
and MCF7, and Jurkat and TPA nuclear extracts are more similar than those 
obtained for other combinations. More importantly, when the peptide profile 
obtained from an independent preparation of Jurkat nuclear extract (labeled '?' in 
the above Table), it received a high score and could be identified as being most 
closely related to the Jurkat cells. 

Applications of Protei n Expression Datasftts 
Relevance to Disease 

As an example of the approach, its potential use in the diagnosis and study of 
human disease is descri bed, for exam ple i n infectious disease or a genetic 
disease such as cancer. The invention may be used to systematically identify, 
compare, classify, and characterize and investigate biological or clinical samples 
from normal and virus- or bacterially-infected cells and tissues, similar cells 
obtained over a course of infection, or similar cells obtained over the course of a 
therapeutic treatment. Similarly, the invention may be used to systematically 
identify, compare, classify, and characterize and investigate biological or clinical 
samples from normal and cancerous cells and tissues, cancerous cells and 
tissues obtained from a variety of related or unrelated liquid or solid tumors, cells 



( 



( 
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The resulting datasets or profiles may therefore (i) identify robust signatures of 
disease states that can be used to facilitate diagnostic and prognostic medical 
procedures, (ii) refine current models of disease and highlight productive areas 
for focusing further basic and applied investigative approaches. 

Uses in Toxicology Studies 

As another example of the use of the invention, quantitative peptide profiles may 
be used for investigation of toxic effects in human or other tissues or cells, for 
instance the side-effects of candidate drug compounds. This is because the 
toxicity may be represented by changes in the expression patterns of peptides 
and proteins in the cells. Currently, such toxic effects are investigated using 
general marker enzymes such as cytochrome oxidase. In many ways, this is a 
'blunt tool\ failing to differentiate between different types of toxicity, and/or the 
severity of the toxic effect. Quantitative peptide profiles are likely to be discrete 
for individual compounds while profiles generated in response to related 
compounds would be expected to be also related to each other. 

A database of profiles can be assembled that describes the protein complements 
of tissues treated with known toxic agents. Large numbers of drug candidates 
can then be screened and their profiles compared to those in the reference 
database. Accordingly, the invention includes methods of determining the toxicity 
of a candidate drug compound. The method comprises administering the 
candidate compound to a cell. As described above, samples suitable for MS 
anaylsis are generated and a peptide profile is produced. Relative abundance of 
peptides in samples is also preferably determined. This candidate compound 
peptide profile is compared to peptide profiles in a database or library (for 
example, profiles showing the cell in a normal state and in varied states of 
toxicity). If the candidate compound sample profile is highly similar to (for 
example, greater than 90%, 95%, or 99% similarity), or identical to a profile in the 
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peptides is also preferably compared io other profiles to determine the amount of 
toxicity of a candidate compound. 

Profiles obtained from drug candidates that are similar to those obtained from 
damaged tissue alert the investigators to potential toxicity problems associated 
with that compound. Because each single profile comprises a large dataset 
(many individual proteins and their relative abundances), comparison of the 
profiles is statistically powerful. This reduces dependence on animal toxicity 
trials, where large numbers of animals may be necessary to obtain statistically 
relevant data. 

Healthy cells, and cells treated with toxic agents, will be analyzed by liquid 
chromatography-tandem mass spectrometry (LC-MS/MS) using a novel semi- 
quantitative approach, resulting in a protein profile for each treatment that serves 
as a signature of the cell state. The profile comprises data relating tens to 
hundreds of individual proteins and therefore represents a highly specific and 
sensitive description of the protein complement of the cell or tissue in that 
particular state. 



c 



Even without knowledge of protein function, the profiles from cells treated with 
novel compounds can be compared to those from healthy cells or cells treated 
withToxic compounds, the "method may therefore be predictive of toxic effects at 
an early stage of drug development. Further, where the test profile matches the 
profile produced by treatment with a characterized compound or family of 
compounds, the mechanism of toxicity may be similar to that produced by the 
reference class. This application of the invention can be applied to any primary or 
transformed cell line, or to tissues obtained from animal models, preferably 
mammalian and more preferably human, or to experimental or clinical samples. 
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leptin. For both treated and untreated samples, over 400 proteins and 900 
peptides were identified. Of these, 170 were uniquely observed in one or other 
sample. In Figure 10, a screenshot of this analysis shows peptides present in 
one or other sample (green or red) and peptides unique to either sample (blue). 
This experiment demonstrates that the invention can be used to examine the 
effect of drugs and other treatments on proteome mixtures. 

Example 7. Peptide profiling to characterize tissue from different organisms 

As further proof of principle, the peptide profiling approach was applied to 
different organisms - two microbes (Escherchia coli and Saccharomyces 
cerevisiae) and two mammals (Homo sapiens - humans and Mus musculis - 
common lab mouse). A standard MCAT LC-MS peptide profiling analysis was 
used to follow expression of hundreds of proteins for each species (Tables 9, 10). 

Table 9. Peptide profiling of microbial species. 





Proteins 


Peptides 


Yeast 


233 


519 


Bacteria 


542 


1647 



When the peptide profiles of the highly divergent microbial species were 
compared, 516 of the 519 yeast proteins were unique. In contrast, when a similar 
analysis was done for peptide profiles of the two mammalian species, 44 of 197 
mouse peptides were similarly observed in the human profile (representing 
homologous protein/peptide species). Thus, these preliminary analyses indicate 
that peptide profiling can both distinguish species, and that the peptide profile 
may reflect the degree of relatedness of organisms (Figure 11). 

Table 10. Peptide profiling of mammalian species. 





Proteins j Peptides | 


Mouse 


142 | 197 
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Because peptide profiling relies on the use of many data points to assess the 
degree of relatedness of many different samples, it is critical that the method be 
reproducible. This is confirmed on the samples described here. One such 
example, involving the peptide profile of yeast whole cell lysate, is shown here 
(Table 11, Table 12). 

Table 1 1. Peptides observed for two repeat samples. 



I Total 


Shared 


Samplel 


776 


686 


Samp!e2 


723 


686 



Table 12. Proteins observed for two repeat samples. .-_ 

'. j Total jSharedM 

Sam plel I 304 j 259; 
Sarnple2 j 288: 2591 

This analysis establishes the reproducibility of the process. 

Figure 12 is a representation of a reference database of protein profiles, 
incorporating both the identity, relative quantities, and overlap of peptides or 
proteins in various samples. 

It will be appreci ated that th e description abo ve re lates to the preferred 
embodiments by way of example only. Many variations on the computer system 
and methods for delivering the invention will be obvious to those knowledgeable 
in the field, and such obvious variations are within the scope of the invention as 
described and claimed, whether or not expressly described. 

All references, including journal articles, patents and patent applications, in this 
application are incorporated by reference herein in their entirety. 
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We claim: 

1 . A method for identifying the constituent proteins for a cell type, tissue or 
pathological sample using a database comprising peptide profile libraries 
wherein the libraries have multiple peptide sequences, comprising: 

a) deriving a plurality of peptides from the cell type, tissue or pathological 
sample; 

b) identifying the peptide species by liquid phase tandem mass spectroscopy 
sequencing; 

c) compiling a data set or peptide profile containing the collection of peptide 
sequences obtained thereby; and 

d) cross-tabulating with a collection of peptide sequences in the database. 

2. The method of claim 1 , wherein the step of deriving a plurality of peptides 
from the cell type, tissue or pathological sample further comprises the step 
of: 

a) obtaining a peptide-containing extract of the cell type, tissue or 
pathological sample; 

b) digesting the extract producing peptides with an enzyme, the enzyme 
capable of localizing mobile protons to the N-terminal amine and the side 
chains of the carboxy-terminal arginine or lysine residues; 

c) separ ating t he peptides b y hig h pressur e liquid chromatography 
apparatus; 

3. The method of claim 2, wherein the enzyme comprises one selected from 
the group consisting of trypsin and endoproteinase LysC. 

4. The method of any of claims 2 to 3, wherein the step of digesting the 
extract producing peptides further comprises the steps of: 

a) dividing the extract into two equal portions; 
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5. A method for identifying a peptide sequence for a cell type, tissue or 
pathological sample using a database comprising peptide profile libraries 
wherein the libraries have multiple peptide sequences, comprising: 

a) obtaining a peptide-containing extract of the cell type, tissue or 
pathological sample; 

b) digesting the extract producing peptides with an enzyme capable of 
localizing mobile protons to the N-terminal amine and the side chains of the 
carboxy-terminal arginine or lysine residues; 

c) separating the peptides by high pressure liquid chromatography 
apparatus; 

d) identifying the peptide species by tandem mass spectroscopy sequencing; 
and 

e) compiling a data set or peptide profile containing the collection of peptide 
sequences obtained thereby. 

6. The method of claim 5, wherein the enzyme comprises one selected from 
the group consisting of trypsin and endoproteinase LysC. 

7. The method of any of claims 5 to 6, wherein the step of digesting the 
extract producing peptides further comprises the steps of: 

a) dividing the extract into two equal portions; 

b) derivatizing completely one of the two equal portions with a reagent, the 
reagent comprising one selected from the group consisting of o- 
methylisourea, homoarginine, canavanine, hydrazine, phenylhydrazine, and 
butyric acid derivatives. 

c) combining the two portions. 

8. A method for quantitating the relative abundance of proteins in two 
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b) identifying the peptide species by tandem mass spectroscopy sequencing; 
compiling a data set or peptide profile containing the collection of peptide 
sequences obtained thereby; 

c) cross-tabulating with a collection of peptide sequences in the database of 
peptide sequences; and 

d) determining the relative abundance of the peptides and/or proteins. 

9. A method for quantitating the relative abundance of proteins in two samples 
of a cell type, tissue or pathological sample using a database comprising 

_. peptide profile libraries wherein the libraries have multiple peptide 
sequences, comprising: 

a) deriving a plurality of peptides from each sample of the cell type, tissue or 
pathological sample; 

b) identifying the peptide species by tandem mass spectroscopy sequencing; 

c) compiling a data set or peptide profile containing the collection of peptide 
sequences obtained thereby; 

d) determining the degree of relatedness of a collection of peptide sequences 
in the database of peptide sequences using clustering and related statistical 
methods 

10. The method of any of claims 8 to 9, wherein the step of deriving a plurality of 
peptides in two samples further comprises the step of: 

a) obtaining a peptide-containing extract of each sample; 

b) digesting, separately the extracts producing peptides with an enzyme, the 
enzyme capable of localizing mobile protons to the N-terminal amine and the side 
chains of the carboxy-terminal arginine or lysine residues; 

c) combining the two extracts; and 

d) separating the peptides by high pressure liquid chromatography. 
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extracts with a reagent, the reagent comprising one selected from the group 
consisting of o-methylisourea, homoarginine, canavanine, hydrazine, 
phenylhydrazine, and butyric acid derivatives. 

13. A method for identifying a peptide sequence for a cell type, tissue or 
pathological sample, comprising: 

a) obtaining a peptide-containing extract of a cell type, tissue or pathological 
sample; 

b) digesting the extract producing peptides with an enzyme capable of 
localizing mobile protons to the N-terminal amine and the side chains of the 
carboxy-terminal arginine or lysine residues; 

c) separating the peptides by high pressure liquid chromatography 
apparatus; 

identifying the peptide species by tandem mass spectroscopy sequencing; and 

d) compiling a data set or peptide profile containing the collection of peptide 
sequences obtained thereby. 

14. The method of claim 13, wherein the enzyme comprises one selected from 
the group consisting of trypsin and endoproteinase LysC. 

15. The method of any of claims 12 to 14, wherein the step of digesting the 
extract producing peptides further comprises the steps of: 

a) dividing the extract into two equal portions; 

b) derivatizing completely one of the two equal portions with a reagent, the 
reagent comprising one selected from the group consisting of o- 
methylisourea, homoarginine, canavanine, hydrazine, phenylhydrazine, and 
butyric acid derivatives. 

c) combining the two portions. 
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cross-tabulated with quantitative data indicating relative and/or absolute 
abundance of each peptide species in a sample; and 

b) a user interface capable of receiving a selection of one or more queries to 
the database for use in determining a rank-ordered similarity of peptide profiles in 
the database. 

17. A method of producing a computer database comprising a computer and 
software for storing in computer-retrievable form a collection of peptide profiles 
for cross-tabulating with data specifying the source of the peptide-containing 
sample from which each. peptide profile was obtained. . . ... 

1 8 - Tne method of claim 1 7, wherein at least one of the sources is from a 
sample known to be free of pathological disorders. 

1 9. The method of claim 1 8, wherein at least one of the sources is a known 
pathological specimen. 

20. A method of comparing quantitative peptide profiles using a database of a 
plurality of peptide profile libraries, the method comprising: 

a) receiving a selection of two or more of the peptide profile libraries; 

b) determining the peptide profiles common to the selected peptide profile 
libraries and identifying profiles unique to each of selected peptide profile library; 
and 

c) displaying the results of the determination. 

21 . The method of claim 20, wherein the correlation of a peptide profile 
against selected peptide profile libraries is determined by 

Px.y = [1/n 0=, to n) I (Xj - Ux ) (Yj - Uy )]/[3 X . d y ] 

where peptides common to two profiles score T and peptides not shared 
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lysine or arginine or the specific products of proteolytic enzymes or chemical 
derivatives of those products, peptides containing rare amino acids, and proteins 
isolated by binding to disease-specific affinity reagents. 

23. The method of claim 22, wherein the specific products of proteolytic 
enzymes comprise chemical derivatives of these products wherein de novo 
sequencing or relative abundance measurements of the peptides is facilitated. 

24. The method of claim 23, wherein the chemical derivatives are obtained by 
guanidinylation and related modifications. 

25. The method of any of claims 21 to 24, wherein the rare amino acids 
comprise tryptophan and cysteine and amino acids comprising 5% or less of the 
amino acid representation. 

26. The method of any of claims 21 to 25, wherein the disease-specific affinity 
reagents comprise polyclonal antibodies, toxin or drugs. 

27. The method of any of claims 21 to 26, wherein the peptide profiles are of 
peptide sequences, the peptide sequences comprising mammalian peptide 
sequences. 

28. The method of any of claims 21 to 26, wherein the peptide profiles are of 
peptide sequences, the peptide sequences comprising microbial peptide 
sequences. 

29. The method of any of claims 21 to 28, wherein the step of receiving a 
selection of two or more of the peptide profile libraries for comparison includes 
receiving a user selection from two or more pull-down menus using a graphical 
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31 . The method of any of claims 21 to 28, wherein the step of receiving a 
selection of two or more of the peptide profile libraries for comparison comprises 
receiving an electronically transmitted file containing sequence and quantitative 
data. 



32. The method of any of claims 21 to 31, wherein the results of the 
determination comprise a unique identifier for related peptide profiles. 

33 - "^e method of any of claims 21 to 31, wherein the results of the 

determination comprise annotated information relating to the related peptide 
profiles obtained from a public database. 

34. The method of any of claims 21 to 31, wherein the results of the 
determination comprise quantitative or relative abundance information relating to 
the related peptide profiles obtained from a public database. 

35. The method of any of claims 21 to 34, further comprising the step of 
displaying the peptide profiles common to the selected peptide profile libraries. 

36. The method of any of cl aims 21 to 34, further comprising Jhe step of 
displaying the peptide profiles unique to the selected peptide profile libraries. 

37 ■ A m ethod of identifying peptide profiles common to a set of environments 
organisms, organs, tissues, cells, cellular fractions or isolated molecular 
complexes using a database comprising peptide profile libraries for a plurality of 
types of organisms wherein the libraries have multiple peptide sequences, the 
method comprising: 

a) displaying at least one list of peptide profile libraries; 



( 
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Figure 4A 
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The peptide or protein profile of different tissues are unique. 
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Bioinformatics techniques reveal general and tissue-specific proteins. 
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Figure 7 
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Similarity dendrogram for different human tissue constructed 
using peptide profiling. 




Figure 8 
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Note: Fract refers to a chromatographic subtraction. 




Figure 9 
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Comparison of peptide profiles for untreated- and leptin-treated human 
muscle ceils. 




Figure 10 
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METHOD TO IDENTIFY CONSTITUENT PROTEINS BY PEPTIDES PROFILING 

CROSS REFERENCE TO RELATED APPLICATION 

This application claims priority from Canadian Patent application No. 2,349,265, 
5 which is incorporated by reference herein. 

FIELD OF THE INVENTION 

The field of this invention relates to the fields of peptide separation and 
proteomics, bioinformatics, metabolite profiling, medicine, drug screening and 
10 computer databases. 

BACKGROUND OF THE INVENTION 

Modern biochemistry and molecular medicine is entering the post-genomic era. 
While genome sequencing has generated a large amount of genetic data, the 
is focus in the biological sciences is now changing to the full characterization of 
proteins. Protein post-translational modifications, protein localization, protein- 
protein interactions, and analysis of protein structure and folding have become 
subjects of major importance. 

20 Proteomics is the study of patterns of protein expression by complex biological 
systems. It involves, in principle, the determination of the relative abundance, 
post-translational Wdification, and/or stability of large numbers of cellular 
proteins at specific time-points within the life cycle of an organism. 



25 There is growing recognition that qualitative and quantitative analysis of protein 
expression profiles on a genome-wide scale will accelerate the development of 
powerful new diagnostic tools and therapeutics, including novel biomarkers and 
drug targets, as well as lead to a better understanding of the basic molecular 
logic that governs cell biology. This is because most, if not all, complex biological 

so processes are ultimately regulated by means of protein turnover and not simply 
through the control of gene expression. 

The study of protein expression will bring researchers closer to the actual 
biological function of genes than studies of gene sequence or gene expression 
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alone. This is because molecular regulation of proteins, and not simply their 
corresponding genes, holds the key to the function of most, if not all, complex 
biological processes. 

5 In contrast to genomics, which captures DNA information that is largely stable 
throughout the lifetime of an organism, proteomics efforts seek to summarize the 
protein-expression patterns of dynamic biological systems at different times! 
While there are a finite number of genes in a given genome, a cell's proteome is 
constantly fluctuating in response to environment and cellular perturbations. 

10 Hence, understanding how proteins work together requires systematic data on 
the entire spectrum of protein status in a cell at any given time. 

Biology Enters the Post-Genomic Era 

is By the late 1990's the DNA sequences of numerous bacterial and eukaryotic 
organisms had been published and in 2000 the nearly complete DNA sequence 
of Homo sapiens was completed. The availability of large-scale genomic 
sequencing efforts now offers investigators a unique opportunity to perform 
comparative analysis from an evolutionary perspective which can both help to 

20 annotate and validate completed genome sequences and also help identify 

conserved protein function, regulation, or pathways based on protein sequence 
homology. 

Today several disciplines, in particular bioinformatics, functional genomics, and 
25 proteomics, are converging in efforts to exploit this newly-available genome 

sequence information. The long-term objective of these efforts is to understand 
the function and interrelationships of the many thousands of genes and proteins 
present in human cells, with the implicit expectation that this understanding will 
lead to dramatic progress in the clinical sciences. 

30 

In the last few years, laboratories have begun to investigate the functions of the 
protein products of genes and their respective regulatory pathways in a 
systematic global manner. Several approaches are now commonly used. First, 
systematic two-hybrid experiments can be used to define interactions among 

2 
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large sets of proteins (Flores et al, 1999), including whole yeast proteome (Ito et 
al., 2000; Uetz et al, 2000). Second, comprehensive screening of mutant genetic 
loci as a means for dissecting networks of interacting gene products has recently 
been adapted to automated high-throughput formats. Finally, powerful 
experimental tools for identifying the components of protein samples, including 
large complexes such as the ribosome (Link et al., 1999) and nuclear pore (Rout 
et al., 2000), and most recently whole organelles and whole cells have been 
described. 



Tandem Mass Spectrometry 

Because the amino acid sequence of a protein is encoded in DNA, and because 
the rules for determining the primary amino acid sequence of a protein are 
known, vast numbers of hypothetical proteins with no known function await 
classification and characterization. Clearly, many of these genes and proteins 
play a role in human disease and other phenomena of biological or commercial 
interest. 



The emerging field of proteomics research relies on enabling technologies that 
can accurately and rapidly characterize the numerous diverse proteins typically 
found in biological samples. This requires scalable, robust, and automated 
methods for protein analysis. — 

To reveal biochemical pathways and regulatory networks, and help define new 
targets for structure-function analysis, proteomics studies require high-resolution, 
high-sensitivity techniques for separation, detection, and quantitation of proteins 
as well as methods for linking proteins to their corresponding cognate gene 
sequences. 

Mass spectrometry (MS) is currently the method of choice for identifying proteins 
present in biological mixtures. The primary advantages of MS are its high- 
sensitivity, accuracy and capacity. 



3 
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Mass spectrometry is the study of gas phase ions as a means to characterize the 
structures, and hence identities, of molecules. Proteomics began with the 
commercialization of soft ionization techniques in the 1990s, in particular 
electrospray ionization (ESI) and matrix assisted laser desorption ionization 

5 (MALDI), which permitted analysis of proteins for the first time. Commercial MS 
instruments are designed as high performance instruments for structural 
characterization of ions produced by these soft ionization techniques and have 
largely replaced traditional Edman chemical sequencing for the analysis of 
proteins. MS has proven to be very successful at identifying limited numbers of 

10 proteins, such as single polypeptide bands cut from polyacrylamide gels, and it is 
currently possible to identify proteins at picomolar to sub picomolar levels. 

Recent advances in mass spectrometry and data analysis described below are 
providing the necessary tools for implementation of high-throughput protein 
is identification and characterization. As the scope of protein analysis has shifted 
from a molecule-by-molecule approach to a genomic scale, the ability of both 
academia and industry to generate new MS data has dramatically outstripped the 
ability to validate, manage, and interrogate the data. 

20 For these studies, routine access to state-of-the-art mass spectrometry 

instrumentation with an adequate infrastructure is essential. Two new ionization 
techniques, MALDI and ESI, have revolutionized the analysis of proteins. The 
MALDI and ESI techniques can be coupled with various types of mass analyzers, 
such as quadrupoles (Quad, Q), time-of-flight (TOF), ion-trap, Fourier transform 

25 ion cyclotron resonance (ICR) and hybrid instruments with two different mass 
analyzers (Q-TOF). Each kind of instrument has advantages and disadvantages 
and, in practice, the achievement of high throughput in conjunction with reliable 
protein identification requires access to both MALDI and ESI instruments. 

30 Mass spectrometry is the most powerful physical technique in its ability to resolve 
and identify rapidly the thousands of proteins expressed by a genome. Mass 
spectrometric techniques are particularly effective when coupled with classical 
biochemical techniques such as proteolytic digestion, immunoprecipitation and 

4 
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separation techniques such as affinity chromatography, HPLC or capillary 
electrophoresis. 

Tandem mass spectrometry (MS/MS) provides a means for fragmenting a mass- 
selected ion and measuring the mass-to-charge ratio (m/z) of the product ions 
that are produced during the fragmentation process. The MS/MS process used 
most often is based on collision-induced dissociation (CID), in which a mass- 
selected ion is transmitted to a high-pressure region of the instrument where it 
undergoes low energy collisions with inert gas molecules. 

As a molecular ion collides, a portion of its kinetic energy is converted into excess 
internal energy rendering the ion unstable, and driving unimolecular 
fragmentation reactions prior to leaving the collision cell. Detailed structural 
information is generated as a result of fragmentation. The mass selectivity of 
many commercial MS systems permit the isolation of single precursor peptide 
.ons from mixtures, thereby removing the contribution of any other peptide or 
contaminant from the sequence analysis step. The product ion spectra can 
subsequently be interpreted to deduce the amino acid sequence of a protein. 

20 A protein to be identified by MS is first digested enzymatically with a site-specific 
protease such as trypsin (which cleaves after lysine and arginine residues) in 
order to produce peptides with structures suitable for MS. Tryptic peptides are 
particularly amenable to MS/MS analysis since mobile protons localize to the N- ( 
term.nal amine and the side chains of the carboxy-terminal arginine or lysine 

25 residues at which proteolysis occurs. These protons cause peptides to fragment 
m a somewhat predictable manner following activation in a tandem MS, leading 
to production of two broad classes of fragment ions - the so-called amino- 
terminal b-type ions and carboxy-terminal y -typ e ions. Recognition of the 
members of these series is a fundamental process of MS-based protein 

30 sequence interpretation. 

Tandem mass spectrometry is a uniquely powerful technology for identifying the 
components of low abundance protein complexes (Andersen et al., 1996) Using 
th,s technique, the molecular weight of individual ionized peptides resulting from 
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trypsin digestion of protein sample is initially determined by the mass 
spectrometer. The peptides are then isolated based on their mass/charge 
properties, fragmented using low energy collision with inert gas (or with 
resonance excitation), and the fragments are analyzed using a second round of 
5 mass spectrometry. 

The relative abundance of daughter product ions in peptide tandem mass spectra 
varies considerably, and some are not observed. This variation reflects subtle 
differences between favored and disfavored fragmentation sites, the nature of the 

10 amino acid side chains, and their position on the peptide backbone. CID of 

protonated peptides also leads to other fragmentation reaction products that can 
complicate spectral interpretation. Molecular losses of water or ammonia for 
instance, are commonly observed in the product ion scans of tryptic peptide ions. 
Spectra often also contain non-peptide noise peaks. Because of this, de novo 

15 interpretation of spectra is extremely difficult to automate and most MS-based 
identification techniques rely on reducing the computational scale of the problem 
by searching protein sequence databases using a relatively simple correlation 
algorithm. 

20 The fragmentation patterns of the peptides can be used to obtain amino acid 
sequence information by comparison with predicted patterns obtained from 
translated protein databases. In addition, advances in tandem mass spectrometry 
mean that polypeptides can now be identified at a low picomolar to femtomolar 
ievei in a rapid, sensitive, and versatile manner. By revealing the composition of 

25 biologically relevant, low abundance protein complexes, the technology can 
provide fundamental insight into the circuitry of interacting proteins. 

Tryptic peptides are particularly amenable to MS/MS analysis since mobile 
protons localize to the N-terminal amine and the side chains of the carboxy- 
30 terminal arginine or lysine residues at which proteolysis occurs. These protons 
cause peptides to fragment in a somewhat predictable manner following 
activation in a tandem MS, leading to production of two broad classes of fragment 
ions - the so-called amino-terminal b-type ions and carboxy-terminal y-type ions 
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(a typical MS/MS peptide spectra showing prominent b- and y-ions is shown 
below). 



The fragmentation pattern reflects the dissociation of the peptides along the 
5 peptide bond backbone, and therefore correlates with the sequence of amino 
acids for those peptides. Recognition of the members of the b- and y-ion series 
is a fundamental process of MS-based protein sequence interpretation. Since de 
novo interpretation of spectra is difficult to automate, most MS-based 
identification techniques rely on reducing the computational scale of the problem 
10 by searching protein sequence databases using a relatively simple correlation 
algorithm. The SEQUEST program (US Patent 5,538,897), for instance, uses 
uninterpreted product ion spectra to search databases of theoretical spectra ( 
derived from protein and translated gene sequence databases. 

15 Recent developments in tandem mass spectrometry (MS/MS) now allow for the 
identification of hundreds of proteins per sample in a single run using available 
technology. This represents a major breakthrough compared to traditional 
methods, for example, 2D gel electrophoresis, and permits, for the first time, 
protein analysis on a truly proteomic scale. 

20 

Accurate mass measurement of peptides derived from proteins provides 
information not available from DNA sequence, such as pbst-tfanslational 
modifications and correction to errors in the DNA databank. Database searching ( 
with masses of peptides obtained from proteolytic digests is a well-established 
25 technique in many laboratories around the world. The searching of databases 
with partial sequence information obtained from MS/MS sequencing experiments 
is even more reliable because it imposes statistical constraints on the 
identification. 

30 The ability of mass spectrometry techniques to quantify the levels of individual 
peptides in a sample has been limiting. Recent approaches, such as ICAT 
(isotope-coded affinity tags; Gygi et al, 2000), have begun to address this issue. 
Using ICAT and similar strategies, the proteins of two samples are differentially 
modified with a reagent that quantitatively adds a molecular tag of defined 
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molecular mass to one of the protein samples. By combining the samples after 
this treatment, the relative abundance of different protein species in each sample 
can be estimated by comparing the signal intensities of the corresponding 
peptides in the mass spectrometer. 

Another quantitative approach, limited to culturable organisms, is to label growth 
media with stable isotopes such as N15. The isotope becomes incorporated into 
the peptide or protein and the isotope-treated peptide is offset in the mass 
spectrum by multiples of 1 amu (the difference in mass between the naturally 
abundant isotope N14 and the heavy isotope derivative N15) depending on the 
number of N atoms in the peptide. These spectra can be deconvolved to 
determine the relative abundance of the labeled and unlabeled peptide species. 
Alternatively, non-isotopic mass tags, whereby the 'labeled' or tagged species is 
offset by the mass of the tag, can be used. Thus methods suitable for high- 
throughput and efficient identification and quantitation of large numbers of 
proteins from complex mixtures are now available. 

HPLC 

High-resolution separation techniques are required to separate the peptide 
components of complex biological mixtures prior to mass spectrometry. A 
particularly powerful approach to identifying the components of complex protein 
mixtures is direct analysis of the protease-digested proteins using high- 
performance, high-resolution multi-dimensional liquid separation techniques 
coupled online to mass spectrometry/database searching (HPLC-MS/MS)(Link et 
al., 1999). This strategy enables the separation of very complex peptide mixtures, 
such as the whole cell extracts or nuclear extracts (Washburn, 2000). One aspect 
of the method separates complex peptide mixtures by strong cation exchange in 
the first dimension and by reverse phase in the second. However, many 
combinations of separation media and more than two dimensions could be used. 
One advantage of the strategy is that it eliminates the need to separate proteins 
on gels or to identify them using antibody- or affinity-based techniques that are 
both time-consuming and difficult to standardize. Therefore this technique 
circumvents the technical and analytical limitations associated with traditional 
proteomics technologies. 

8 
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Bioinformatics 

The interpretation of peptide mass spectra for the purposes of generating protein 
identifications can be carried out manually but requires experience and skill and 
is prohibitively time-consuming. For this reason, computer algorithms have been 
developed that, while not capable of interpreting all spectra they encounter, can 
easily outperform human identifications for even minimally complex peptide 
mixtures. Any of several generally available algorithms may be used for this 
purpose. For instance, the SEQUEST program (Eng et al., 1994) uses 
uninterpreted product ion spectra to search databases of theoretical spectra 
derived from protein and translated gene sequence databases. SEQUEST first 
generates a list of theoretical peptide masses for each entry in the database that 
match the experimentally determined peptide mass, producing a list of candidate 
peptides. The program then calculates the fragment ion masses expected for 
each of the candidate peptides, generating a predicted MS/MS spectrum. Finally, 
the experimentally determined MS/MS spectrum is compared with the predicted 
spectra using a correlation function. Each comparison receives a score, and the 
highest-scoring peptide(s) are reported. When high scoring matches are 
detected, one effectively jumps from spectral data directly to a peptide identity, 
which in turn can be linked to the entire amino acid and DNA sequence of the 
corresponding gene. Ideally, a protein is positively identified when the spectra of 
one or more peptides in a tryptic digest can be matched unambiguously. 

Mass spectral reference libraries representing stored tandem mass spectra, or 
validated chemical signatures, are routinely used for the identification of small 
chemical compounds by MS (eg. Wiley Registry, NIST database). Unknown 
compounds can then be both identified by searching experimental spectra 
against a comprehensive database of these reference mass spectra, which are in 
turn derived from pure compounds, so that only hits of strong similarity or identity 
are produced. A similar reference spectral database approach would likewise 
facilitate MS-based identification of proteins. 

Compared to mRNA expression analysis the development of corresponding 
'proteomics' technologies has lagged, with only a few laboratories addressing 
complex phenotypes on a global scale. Nonetheless, protein expression profiling 
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holds great promise for rapid genome functional analysis. It is plausible that the 
protein expression profile could serve as a universal and rich cellular phenotype: 
provided that the cellular response to disruption of different steps of a given 
biochemical process or pathway is similar, and that there are sufficiently unique 
5 cellular responses to the perturbation of most cellular pathways, systematic 
characterization of novel genetic mutants could be carried out with a single 
genome-wide protein expression measurement 

To date the only studies focusing on peptides or proteins that includes a 
10 quantitative component has been the separation of bacterial and yeast cell 

lysates on 2-dimensional electrophoretic gels (refe). These approaches do not 
directly identify the resolved proteins, are relatively insensitive, and are unlikely to 
scale up to the study of larger proteomes (e.g. that of vertebrates). Furthermore, 
no attempt was made to use the data to identify or characterize unknown 
is samples. 

SUMMARY OF THE INVENTION 

The protein profiling approach proposed has both a qualitative and a quantitative 
component such that each profile generated can be directly compared to other 
20 profiles present in a reference database. 

This invention describes the use of peptide profiling to identify, characterize, and 
classify biological samples. In complex samples, many thousands of different 
peptides will be present at varying concentrations. The invention uses liquid 

25 chromatography and similar methods to separate peptides, which are then 

identified and quantified using mass spectrometry. By identification it is meant 
that the correct sequence of the peptide is established through comparisons with 
genome sequence databases, since the majority of peptides and proteins are 
unannotated and have no ascribed name or function. Quantification means an 

30 estimate of the absolute or relative abundance of the peptide species using mass 
spectrometry and related techniques including, but not limited to, pre- or post- 
experimental stable or unstable isotope incorporation, molecular mass tagging, 
differential mass tagging, and amino acid analysis. 

10 
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The principle experimental strategy of the present invention is centered on rapid 
high-throughput protein identification using coupled tandem mass spectrometry 
(MS/MS) and sequence database searching. Quantitation is based on either 
metabolic labeling with stable isotopes or with chemical derivation. Below, an 
example of a non-isotopic tag based on the lysine^specific guanidylation reagent 
O-methylisourea is described in detail. Significant patterns of peptide expression 
are identified with software and data mining algorithms. Below, a method is 
described for identifying, classifying and characterizing functions of known and 
unknown gene products, peptides and proteins , for characterizing metabolic and 
other functional pathways in cells, and for identifying the proteins and pathways 
targeted by drugs and other reagents. The method is based on the comparison 
of protein profiles obtained following global proteomics or other comprehensive 
protein studies from cells, cell fractions, tissues, organisms or other defined 
sources. 

The invention further contemplates the use of high-throughput robotic screening 
of diverse chemical compound libraries to systematically identify small molecules 
that perturb cellular pathways associated with disease. The protein targets of the 
lead compounds will be isolated and identified by the tandem mass spectrometry 
profiling techniques described herein. Protein profiling acts as an optimal assay 
since the profile of a healthy cell or tissue is the goal. 

Theinvention = relates tO"a-m"eth6iaT6riae7itifyinf*the constituent proteins for a cell 
type, tissue or pathological sample using a database comprising peptide profile 
libraries wherein the libraries have multiple peptide sequences, comprising: 

1. deriving a plurality of peptides from the cell type, tissue or pathological 
sample; 

2. identifying the peptide species by liquid phase tandem mass spectroscopy 
sequencing; 

3. compiling a data set or peptide profile containing the collection of peptide 
sequences obtained thereby; and 

4. cross-tabulating with a collection of peptide sequences in the database. 

The step of deriving a plurality of peptides from the cell type, tissue or 
pathological sample preferably further comprises the step of: 
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a) obtaining a peptide-containing extract of the cell type, tissue or pathological 
sample; 

b) digesting the extract producing peptides with an enzyme, the enzyme 
capable of localizing mobile protons to the N-terminal amine and the side 

5 chains of the carboxy-terminal arginine or lysine residues; 

c) separating the peptides by high pressure liquid chromatography apparatus; 

The enzyme preferably comprises one selected from the group consisting of 
trypsin and endoproteinase LysC. The step of digesting the extract producing 
10 peptides preferably further comprises the steps of: 

a) dividing the extract into two equal portions; 

b) derivatizing completely one of the two equal portions with a reagent, the 

reagent comprising one selected from the group consisting of o- 
methylisourea, homoarginine, canavanine, hydrazine, phenylhydrazine, and 
15 butyric acid derivatives. 

c) combining the two portions. 

The methods of the invention may be used in toxicology analysis. The methods 
optionally comprise administering a candidate compound to a cell. As described 
20 above, samples suitable for MS anaylsis are generated and a peptide profile is 
produced. Relative abundance of peptides in samples is also preferably 
determined. This candidate compound peptide profile is compared to peptide 
profiles in a database or library (for example, profiles showing the cell in a normal 
state and in varied states of toxicity). If the candidate compound sample profile is 

25 highly similar to (for example, greater than 90%, 95%, or 99% similarity), or 
identical to a profile in the database or library, then that similarity shows the 
amount of toxicity of the candidate compound to the cell. If the candidate 
compound sample profile is highly similar to a normal cell profile, then the 
candidate compound is less likely to be toxic than if the candidate compound 

30 sample profile is similar to the peptide profile of the cell in state of toxicity. The 
relative abundance of the test sample peptides is also preferably compared to 
other profiles to determine the amount of toxicity of a candidate compound. In a 
similar manner, candidate drugs compounds may be screened against cells, such 
as diseased cells. If the candidate drug shifts the profile from a disease profile 

12 
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and relative abundance towards a normal, healthy profile and relative abundance 
with substantial similarity (eg. Over 90%, 95%, 95% similarity), or identical to the 
healthy profile and relative abundance, the drug compound is likely to be useful 
as a therapeutic. 



Another embodiment relates to a method for identifying a peptide sequence for a 
cell type, tissue or pathological sample using a database comprising peptide 
profile libraries wherein the libraries have multiple peptide sequences, 
comprising: 

a) obtalning a peptide-containing extract of the cell type, tissue or pathological 

sample; 

b) digesting the extract producing peptides with an enzyme capable of localizing ( ~ 

mobile protons to the N-terminal amine and the side chains of the carboxy- 
terminal arginine or lysine residues; 

c) separating the peptides by high pressure liquid chromatography apparatus; 

d) identifying the peptide species by tandem mass spectroscopy sequencing; ' 

and 

e) compiling a data set or peptide profile containing the collection of peptide 

sequences obtained thereby. 
The enzyme is preferably selected from the group consisting of trypsin and 
endoproteinase LysC. The step of digesting the extract producing peptides 
preferably-further comprises" th^ steps bf: 

a) dividing the extract into two equal portions; C 

b) derivatizing completely one of the two equal portions with a reagent, the 

reagent comprising one selected from the group consisting of o- 
methylisourea, homoarginine, canavanine, hydrazine, phenylhydrazine, and 
butyric acid derivatives. 

c) combining the two portions. 

Another aspect of the invention includes a method for quantitating the relative 
abundance of proteins in two samples of a cell type, tissue or pathological 
sample using a database comprising peptide profile libraries wherein the libraries 
have multiple peptide sequences, comprising: 
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a) deriving a plurality of peptides from each sample of the cell type, tissue or 
pathological sample; 

b) identifying the peptide species by tandem mass spectroscopy sequencing; 

c) compiling a data set or peptide profile containing the collection of peptide 
5 sequences obtained thereby; 

d) cross-tabulating with a collection of peptide sequences in the database of 
peptide sequences; and 

e) determining the relative abundance of the proteins. 

In the methods of the invention, a pathological sample may have been contacted 
10 with a candidate drug compound and the peptide profile and/or relative 

abundance of the peptides and/or proteins is compared to a database comprising 
peptide profile libraries of the cell in varied states of toxicity (ie. exposed to known 
toxic compounds which injure and/or kill the cell). The toxicity of the candidate 
drug compound may be determined by comparison of the profile and relative 
15 abundance for the cell type, tissue or pathological sample exposed to the 
candidate drug compound with the profile and relative abundance for the cell 
type, tissue or pathological sample in varied states of toxicity and a normal state. 
A similar method may be used to determine whether a compound is likely to be 
useful as a therapeutic, for example by comparison of the profile and relative 
20 abundance for a pathological (diseased) cell type, tissue or sample exposed to 
the candidate drug compound with the profile and relative abundance for the cell 
type, tissue or sample in a normal, healthy state. 

The invention includes a method for quantitating the relative abundance of 
25 proteins in two samples of a cell type, tissue or pathological sample using a 
database comprising peptide profile libraries wherein the libraries have multiple 
peptide sequences, comprising: 

a) deriving a plurality of peptides from each sample of the cell type, tissue or 
pathological sample; 

30 b) identifying the peptide species by tandem mass spectroscopy sequencing; 
c) compiling a data set or peptide profile containing the collection of peptide 
sequences obtained thereby; 
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d) determining the degree of relatedness of a collection of peptide sequences 
in the database of peptide sequences using clustering and related statistical 
methods 

The step of deriving a plurality of peptides in two samples preferably further 
5 comprises the step of: ----- 

a) obtaining a peptide-containing extract of each sample; 

b) digesting separately the extracts producing peptides with an enzyme, the 

enzyme capable of localizing mobile protons to the N-terminal amine and 
the side chains of the carboxy-terminal arginine or lysine residues; 
10 c) combining the two extracts; and 

d) separating the peptides by high pressure liquid chromatography. 

The enzyme preferably comprises one selected from the group consisting of 

trypsin and endoproteinase LysC. 



15 
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The step of digesting the extracts preferably further comprises the step of 
derivatizing completely one of the two extracts with a reagent, the reagent 
comprising one selected from the group consisting of o-methylisourea, 
homoarginine, canavanine, hydrazine, phenylhydrazine, and butyric acid 
derivatives. 

The invention also includes a method for identifying a peptide sequence for a cell 
type, tissue or pathological sample, comprising: 

a) obtaining a peptide-containing extract of a cell type, tissue or pathological 

sample; 

b) digesting the extract producing peptides with an enzyme capable of localizing 

mobile protons to the N-terminal amine and the side chains of the carboxy- 
terminal arginine or lysine residues; 

c) separating the peptides by high pressure liquid chromatography apparatus- 

d) identifying the peptide species by tandem mass spectroscopy sequencing- 

and 

e) compiling a data set or peptide profile containing the collection of peptide 

sequences obtained thereby. 
The enzyme preferably comprises one selected from the group consisting of 
trypsin and endoproteinase LysC. 
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The step of digesting the extract producing peptides preferably further comprises 
the steps of: 

a) dividing the extract into two equal portions; 

b) derivatizing completely one of the two equal portions with a reagent, the 
reagent comprising one selected from the group consisting of o- 
methylisourea, homoarginine, canavanine, hydrazine, phenylhydrazine, and 
butyric acid derivatives. 

c) combining the two portions. 

Another embodiment of the invention is a computer system for identifying 
quantitative peptide profiles, comprising: 

(a) a database including peptide profile libraries for a plurality of types of 
organisms wherein the libraries have multiple peptide profiles each profile 
comprising an array of at least 50 peptide species each having a unique 
identifier cross-tabulated with quantitative data indicating relative and/or 
absolute abundance of each peptide species in a sample; and 

(b) a user interface capable of receiving a selection of one or more queries to 
the database for use in determining a rank-ordered similarity of peptide 
profiles in the database. 

The invention includes a method of producing a computer database comprising a 
computer and software for storing in computer-retrievable form a collection of 
peptide profiles for cross-tabulating with data specifying the source of the 
peptide-containing sample from which each peptide profile was obtained. 
Optionally, at least one of the sources is from a sample known to be free of 
pathological disorders. Optionally, at least one of the sources is a known 
pathological specimen. 

The invention also includes a method of comparing quantitative peptide profiles 
using a database of a plurality of peptide profile libraries, the method comprising: 

a) receiving a selection of two or more of the peptide profile libraries; 

b) determining the peptide profiles common to the selected peptide profile 
libraries and identifying profiles unique to each of selected peptide profile 
library; and 

c) displaying the results of the determination. 

16 
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The correlation of a peptide profile against selected peptide profile libraries may 
be determined by 

P X .y = [1/n (M t0 n) I (Xj - Px) (Yj - Py)]/[d X . dy ] 

where peptides common to two profiles score '1' and peptides not shared 

between profiles score '0'. 
The peptides profiles are preferably of cell fractions, the cell fractions comprising 
high molecular weight proteins, soluble proteins, membrane proteins, modified 
proteins, phosphoproteins, peptides terminating in lysine or arginine or the 
specific products of proteolytic enzymes or chemical derivatives of those 
products, peptides containing rare amino acids, and proteins isolated by binding 
to disease-specific affinity reagents. 

The specific products of proteolytic enzymes may be comprise chemical 
derivatives of these products wherein de novo sequencing or relative abundance 
measurements of the peptides is facilitated. 

The chemical derivatives may be obtained by guanidinylation and related 
modifications. The rare amino acids may comprise tryptophan and cysteine and 
amino acids comprising 5% or less of the amino acid representation. 
The disease-specific affinity reagents may comprise polyclonal antibodies, toxin 
or drugs. The peptide profiles may be of peptide sequences, the peptide 
sequences comprising mammalian peptide sequences. Thee peptide profiles 
may be of peptide sequences, the peptide sequences comprising microbial 

peptide sequences. ™ - — 

The step of receiving a selection of two or more of the peptide profile libraries for 
comparison may include receiving a user selection from two or more pull-down 
menus using a graphical user interface. The step of receiving a selection of two 
or more of the peptide profile libraries for comparison may comprise command 
line entry using a computer. The step of receiving a selection of two or more of 
the peptide profile libraries for comparison may comprise receiving an 
electronically transmitted file containing sequence and quantitative data. The 
results of the determination may comprise a unique identifier for related peptide 
profiles. The results of the determination may comprise annotated information 
relating to the related peptide profiles obtained from a public database. The 
results of the determination may comprise quantitative or relative abundance 
information relating to the related peptide profiles obtained from a public 
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database. The method may further comprise the step of displaying the peptide 
profiles common to the selected peptide profile libraries. The method may further 
comprise the step of displaying the peptide profiles unique to the selected peptide 
profile libraries. 

5 

The invention also includes a method of identifying peptide profiles common to a 
set of environments, organisms, organs, tissues, cells, cellular fractions or 
isolated molecular complexes using a database comprising peptide profile 
libraries for a plurality of types of organisms wherein the libraries have multiple 
10 peptide sequences, the method comprising: 

(a) displaying at least one list of peptide profile libraries; 

(b) receiving a selection of one or more peptide profile libraries from at least 
one list of peptide profile libraries; 

(c) determining peptide profiles common to the selected peptide profile 
15 libraries; and 

(d) displaying the results of said determination. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Embodiments of the invention will be described by way of example and 
20 with reference to the drawings in which: 

FIG- 1 is a diagram of the MCAT approach for peptide sequencing and relative 
protein abundance determination. 

Fig. 2 is diagram showing how MCAT enables identification and quantitation of 
complex protein mixtures. 
25 Figs. 3A and 3B are diagrams showing de novo sequencing of a yeast peptide 
and a human peptide using MCAT approach. 

Figs. 4A and 4B are diagrams showing relative abundance ratios of positively- 
identified peptides. 

Fig. 5 is a peptide profile generated by a one-dimensional LCMS from diverse 
3p human tissues. 

Fig. 6 shows proteins identified using MCAT based peptide profiling of seven 
human tissues. 
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Fig. 7 shows the differences between protein expression of the seven human 
tissues highlighted by applying agglomerative clustering algorithms. 
Fig. 8 is a similarity dendrogram for different human tissue constructed using 
peptide profiling. 

5 Fig. 9 is a comparison of peptide profiles of different cell compartments. 

Fig. 10 is a comparison of peptide profiles for untreated and leptin-treated human 
muscle cells. 

Fig. 11 shows peptide profiling to distinguish species. 

Fig. 12 is a representation of a reference database of protein profiles. 

10 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

A Quantitative Peptide Profile serves as a precise fingerprint of peptides that can 
be successfully isolated, identified and quantified from the myriad of proteins 
expressed in cells under any given condition. This profile, in turn, can serve as a 
15 unique identifier of cell state. This document describes a method to use 

quantitative peptide profiles to compare biological samples, from any tissue or 
cell, among different types of cell (e.g. nervous tissue cells), or even in samples 
where little or no mRNA is made (e.g. blood platelet cells). 

20 The present invention is distinct from the established method of mRNA 
expression profiling in three important respects. 



First, as mentioned above, the relative abundance of an mRNA is not predictive 
of the abundance of the corresponding protein or cognate peptides. This is 
25 because many factors affect protein expression subsequent to the event of 
mRNA production, including splicing, protein terminal processing, protein 
localization, protein degradation, protein modification, codon usage, the levels of 
available amino acids and the subcellular localization of the protein . mRNA 
expression profiling is unable to account for or predict these events. 

30 

Second, the technology used to acquire mRNA and peptide expression data is 
fundamentally different, the former using nucleic acid hybridization and 
fluorometric quantitation, with the latter, in this embodiment of the invention, 
using mass spectrometry and related ionization techniques. The invention 
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includes a method for detecting and quantitatively analyzing peptides in a 
biological sample, comprising: 

a) obtaining a biological sample in a form suitable for coded abundance tagging; 

b) identifying and quantitating the peptides in the sample by mass coded 
5 abundance tagging. 

In one aspect, the method involves: 

obtaining an extract of the biological sample, such as a cell extract, 
digesting the sample, preferably with an enzyme, such as trypsin, to generate 
10 peptides with a terminal amine group, such as a terminal lysine, 

contacting the peptides with mass differential reagent, such as a guanidination 
compound (eg. Lysine guanidination compound, such as o-methylisourea, which 
modifies the epsilon-amine of the C-terminal lysine\), 

separating the peptides, preferably with liquid chromatography, such as high 
15 throughput capillary liquid chromatography, and 

generating mass spectra for the peptides, preferably with electrospray tandem 
mass spectrometry. 

The method is preferably carried out in both orientations, with a sample divided in 
two and either modified or unmodified. Peptides are alternatively unmodified and 
modified with o-methylisourea differ by the mass differential encoded by the mass 
differential reagent (e.g. 42 amu for O-methylisourea). The method preferably 
further involves sequencing the peptides and/or determining relative abundance 
of the peptides. Methods of sequencing and determining relative abundance are 
described below. Sequencing preferably involves comparing pair-wise sets of 
spectra (MS/MS spectra) to identify identities of y-ion peaks. One can use a 
short sequence of contiguous amino acid sequence from a peptide (e.g. 5-10 
amino acids or greater than 10 amino acids) to identify a corresponding protein. 

30 For identified peptides, the single ion intensity profile is reconstructed from the full 
scan data and the relative abundance of modified peptides is determined by 
integrating the area under the curve. 
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The invention includes a method of identifying a test sample by obtaining a 
peptide profile for the test sample, preferably by MS. This peptide profile is then 
compared to peptide profiles in a database or library to determine if the test 
sample profile is highly similar to (for example grater than 90%, 95% or 99% 
5 similarity) to a profile in the database or library. Relative abundance information 
may similarly be used to identify the test sample. 

The methods of the invention may be used in toxicology analysis. The methods 
optionally comprise administering a candidate compound to a cell. As described 

10 above, samples suitable for MS anaylsis are generated and a peptide profile is 
produced. Relative abundance of peptides in samples is also preferably 
determined. This candidate compound peptide profile is compared to peptide 
profiles in a database or library (for example, profiles showing the cell in a normal 
state and in varied states of toxicity). If the candidate compound sample profile is 

15 highly similar to (for example, greater than 90%, 95%, or 99% similarity), or 
identical to a profile in the database or library, then that similarity shows the 
amount of toxicity of the candidate compound to the cell. If the candidate 
compound sample profile is highly similar to a normal cell profile, then the 
candidate compound is less likely to be toxic than if the candidate compound 

20 sample profile is similar to the peptide profile of the cell in state of toxicity. The 
relative abundance of the test sample peptides is also preferably compared to 
other profiles to determine the amount of toxicity of a candidate compound. In a 
similar manner, candidate drugs compounds may be screened against cells, such 
as diseased cells. If the candidate drug shifts the profile from a disease profile 

25 and relative abundance towards a normal, healthy profile and relative abundance 
with substantial similarity (eg. Over 90%, 95%, 95% similarity), or identical to the 
healthy profile and relative abundance, the drug compound is very likely to be 
useful as a therapeutic. 



30 



Although mRNA expression profiles from cells treated with different drugs have 
been compared to each other in order to determine which existing profile most 
closely matches a 'novel' profile (Hughes et al., 2000), this approach has been to 
date confined to one type of organism, the yeast Saccharomyces cerevisiae. 
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* • 

Using a comprehensive database of reference peptide expression profiles, the 
pathway(s) perturbed as a consequence of an uncharacterized mutation, 
pharmaceutical treatment, or developmental or disease state would be 
ascertained by simply asking which expression patterns in the database the 

5 resulting profile most strongly resembles. The database or library will include one 
or more profiles and/or relative abundance determination and may be electronic 
or in a hard copy form. A sufficiently large and diverse set of profiles obtained 
from different mutants, chemical treatments, and environmental conditions would 
also result in a relatively comprehensive identification of coordinate protein 

10 expression sub patterns, allowing hypotheses to be drawn regarding the 

functions of gene products based on their relationship to other proteins (Eisen et 
al., 1998). 

There are several advantages to this profiling approach compared to the analysis 
15 of single peptides or proteins. First, there is no requirement for prior knowledge 
about the functions of the responsive peptides or parental proteins. Second, 
protein functions deduced from comparisons of profiles in a database can be 
derived from very subtle physiological responses. For instance, even though 
peptide levels may change only slightly in response to an experimental treatment, 
20 coordinate changes among many measured peptide abundances can be 
sufficient to characterize that phenotype. The large numbers of peptides 
measured make it unlikely that an unrelated physiological state will have an 
identical profile, even though this may not be apparent when using conventional 
experiments that measure the levels of one or a few proteins. Third, closely 
25 related profiles can be classed together, thus improving our understanding of the 
underlying biological basis of the classifications. 

The invention includes proteins, including drugs, and other compounds identified 
using methods of the invention. 
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Examples 

Example 1: Measurement of protein relative abundance in complex mixtures 

the method relies on modification of peptides at e-amine of lysine residues with 
O-methylisourea. Peptides so modified can be readily detected by mass 
spectrometry because their mass is increased by 42Da (per lysine residue in the 
sequence). Therefore, the relative abundance of a single peptide from two 
different samples can be determined following differential modification with O- 

methylisourea by comparing the signal intensities for the pair in a mass 

spectrometer. 

The steps of the MCAT procedure are as follows (Fig.1): 

(1) Two protein mixtures, obtained following different experimental treatments 
of a sample, are digested enzymatically with trypsin. 

(2) One digest is treated with O-methylisourea and the other with control 
buffer. 

(3) The digests are desalted using ZipTip reverse phase extraction. 

(4) The two mixtures are combined and analyzed by automated electrospray 
LC -MS/ MS. Using eithe r one-dimensional ( reverse phase) or t wo- 
dimensional (cation exchange and reverse phase) liquid chromatography, 
the peptides are separated as they are introduced to the mass 
spectrometer. The instrument is run in automated multistage mode, 
whereby the following cycle is implemented. First, a full MS scan (400- 
1600 m/z) is used to record the relative intensities of peptide ions 
emerging from the column. Next, MS/MS scans of selected ions are used 
to collect spectra suitable for peptide identification. The instrument then 
reverts back to full scan mode, but is programmed to exclude MS/MS 
analysis of ions that have been identified in the previous cycle(s). 

(5) The MS/MS spectra are used to identify the peptides using protein 
database searching algorithms. 
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(6) For identified peptides, the single ion intensity profile is reconstructed from 
the full scan data and the relative abundance of modified and unmodified 
peptides calculated by integrating the area under the curve. 

In order to correct for systemic errors, for instance preferential labeling by O- 
methylisourea of one sample, the experiment is carried out in both orientations, 
that is both samples are divided in two and either modified or unmodified. The 
fractions are then combined with the corresponding modified or unmodified 
fracton from the other sample. 



Table 1 shows some top scoring peptides from this analysis and their relative 
abundance as estimated by the area-under-curve of their respective selected ion 
tracings. For nearly all peptides, the ratio of unmodified to modified signal is 
slightly less than the expected 1:1. The variation from ideal 1:1 ratio is not the 
result of reduced ionization efficiency or MS signal of the modified peptides 
relative to their unmodified forms because the effect was consistently observed in 
subsequent experiments independently of which sample was chosen for 
modification. More likely, it results from preferential recovery of unmodified 
peptides during the Zip Tip desalting step. 

For this reason, when comparing two samples A and B using the MCAT 
procedure, four mass spectrometry analyses are routinely carried out: I) A versus 
A mod , II) A versus B mod , III) B versus B mod , and IV) B versus A mod . The ratios of 
unmodified to modified peptide signais obtained in i and Hi were used to 
normalize II and IV respectively, and the combination of III and IV served to 
independently confirm the quantitative observations. 
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Table 1. Identification and quantitation of peptides from a yeast whole cell digest 



Protein 


Peptide 


Z a 


Score" 


Observed 
ratio 


Expected 
ratio 


-YLR044C 


AQYNEIQGWDHLSLLPTF 
GAK (SEQ ID NO:1) 


2 


2.3993 


1:0.29 


1:1 


YLR044C 


TTYVTQRPVYLGLPANLV 
DLNVPAK (SEQ. ID. NO:2) 


2 


2.6639 


1:0.2 


1:1 


YLR044C 


KLIDLTQFPAFVTPMGK 
(SEQ ID NO:3) 


2 


3.3881 


1:0.67 


1:1 


YHR174W 


WLTGVELADMYHSLMK 
(SEQ ID NO:4) 


2 


4.0552 


1:0.73 


1:1 


YHR174W 


GVMNAVNNVNNVIAAAFV 
K (SEQ ID NO:5) 


2 


3.2283 


1:0.48 


1:1 


YBR118W 


TLLEAI DAI EQ PSRPTDKP 
LRLPLQDVYK (SEQ ID 
NO:6) 


3 


3.3888 


1:0.63 


1:1 


YBR118W 


VETG VI KPGM WTFAPAG 
VTTEVK (SEQ ID NO:7) 


2 


2.5458 


1:0.23 


1:1 


YEL034W 


VHLVAIDIFTGK (SEQ ID 


1 


3.0798 


1:0.15 


1:1 




NO:8) • — ' 






YKL060C 


SPIILQTSNGGAAYFAGK 
(SEQ ID NO:9) 


2 


3.6709 


1:0.73 


1:1 


YCR012W 


ALENPTRPFLAILGGAK 
(SEQ ID NO:10) 


2 


2.7650 


1:0.33 


1:1 


YDR441C 


GFVPIRRVGKLPGEC* 
(SEQ ID NO:11) 


2 


1.1770 


1:1.07* 


1:1 


YGR192C 


VINDAFGIEEGLMTTVHSL 
TATQK (SEQ ID NO: 12) 


2 


3.1456 


1:0.31 


1:1 



a. Peptide charge 

b. SEQUEST Cross-correlation score 
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Next, mixtures derived from yeast whole cell extracts containing varying 
proportions of MCAT-treated and MCAT-untreated sample were analyzed (Fig. 
2). 

5 Relative abundance signal from five peptides with high SEQUEST scores 

showed linearity across two orders of magnitude (Fig. 2). Beyond this range, the 
weaker signal of the two abundances is indistinguishable from background noise. 

Table 2 shows variation in the measured relative abundance for two peptides 
10 from the same parent protein (and therefore are present in equimolar 
concentrations) in three repJicate experiments. Experiment-to-experiment 
variation for these peptides is within 25% and variation within a single experiment 
for peptides derived from the same protein is within 20% (Table 2). 

15 Table 2. Identification and quantitation of two peptides derived from YLR044C in 
three replicate experiments (A, B, C). 



Protein 


Peptide 


Ratio 
A: A 


Ratio 
A:B 


Ratio 
A:C 


YLR044C 


KLIDLTQFPAFVTPM 
GK (SEQ ID NO:13) 


1.00:1.00 


1.00:0.78 


1 .00:0.87 


YLR044C 


AQYNEIQGWDHLSL 
LPTFGAK (SEQ ID 
NO: 14) 


1.00:1.00 


1.00:0.79 


1.00:1.03 



Ratio of unmodified to modified peptides (normalized to A:A) 

20 

This invention also includes computer systems including software and hardware 
to implement the above methods. Such systems include a database with the 
peptide profiles. 
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Example 2: De Novo Peptide Sequencing and Quantitative Profiling of Complex 
Protein Mixtures Using Mass Coded Abundance Tagging 

Introduction 

5 There is growing recognition that qualitative and quantitative analysis of proteins 
on a genome-wide scale will accelerate the development of powerful new 
diagnostic tools and therapeutics, and lead to a better understanding of the 
molecular logic that governs cell behavior. This is because regulation of protein 
abundance holds the key to the proper function of most biological processes 
io (Pandey & Mann, 2000). Proteomics studies depend on scalable, robust, and 
automated methods for protein identification and quantitation that can routinely 
characterize the numerous diverse proteins typically found in biological samples. 

Mass spectrometry (MS) is currently the technology of choice for identifying 
is proteins present in biological mixtures. The primary advantages of MS are its 
h.gh sensitivity, accuracy and capacity. Tandem mass spectrometry (MS/MS) 
provides a means for fragmenting mass-selected precursor peptide ions and 
measuring the mass-to-charge ratio (m/z) of any product daughter ions produced 
(Andersen et al., 1996). The process usually produces two principle classes of 
20 fragment ions, the so-called N-terminal b-type ions and C-terminal y-type ions 
Informative high quality MS/MS spectra of tryptic peptides typically show 
. PromjnenLb^and^^ 

MS/MS analysis since mobile protons that stimulate the fragmentation process 
readily associate with the side chains of the C-terminal arginine or lysine residues 
25 at which proteolysis occurred ' 

If accurate sequence information is available, computer database search 
algorithms can rapidly and accurately identify proteins analyzed by MS/MS (Eng 
et al., 1994; Mann & Wilm, 1994; Taylor & Johnson, 1997, Qin et al., 1997) in 
30 effect linking the spectra to a corresponding cognate protein or DNA sequence 
When combined with recent developments in tandem mass spectrometry this 
approach allows for routine identification of dozens to hundreds of proteins in a 
single analysis. However, because the possibility of alternative splicing, mutation 
and/or post-translational modification is likely to be a significant feature of the 
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proteomes of higher organisms, a facile peptide sequencing method that is 
independent of sequence databases is desirable. 



Manual interpretation of peptide MS/MS spectra for the purposes of protein 
5 identification (a process usually referred to as de novo sequencing) is often 
prohibitively challenging. Factors such as variation in favored fragmentation sites, 
the effects of the chemical nature of the amino acid side chains and their relative 
order in a peptide backbone, and the presence of side-products such as neutral 
loss ions and non-peptide noise peaks. To address this issue, Mann and 
10 coworkers pioneered a post-experiment stable isotope labeling strategy whereby 
the C-termini of tryptic peptides are labeled with deuterated water in order to 
reduce spectral complexity. Comparison of the modified and unmodified peptide 
MS/MS product ion spectra allows the C-terminal y-ions to be readily 
distinguished and, hence, the peptide sequence discerned. The impact of this 
15 approach has been restricted, however, by the prohibitive cost of the stable 
isotope and the high mass resolution required to distinguish the labeled products. 

Functional genomics studies using DNA microarray technologies have been used 
successfully to compare the abundance of thousands of mRNA species from 
distinct cell states. In contrast, only limited analogous quantitative data has been 
obtained for protein abundance. As the scope of protein analysis has shifted from 
a molecule-by-molecule approach to a genomic scale, the ability to generate 
quantitative protein data has lagged considerably. Chait and coworkers reported 
the potential of stable N 15 isotope labeling of proteins as a means to determine 
the relative abundance of select subsets of proteins isolated from cultured yeast 
cells (Oda et al., 1999). As the isotope becomes incorporated, the mass of the 
protein becomes offset in a mass spectrum by multiples of 1 amu (the difference 
in mass between the naturally abundant N 14 isotope and the heavy N 15 isotope 
derivative) depending on the number of labeled N atoms. Although powerful, this 
approach is restricted to organisms that can be grown in defined media. 

Aebersold and coworkers recently introduced an alternative protein quantitation 
strategy based on post-experiment stable isotope labeling (Gygi et al, 1999). The 
ICAT (isotope-coded affinity tag) chemistry uses isotopic variants of a biotin- 
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containing moiety to differentially label cysteine-containing peptides as a means 
to obtain relative abundance data for proteins found in two distinct samples in a 
single analysis. Other approaches based on differential stable isotope labeling 
have been devised (Munchbach et al., 2000). The ICAT method is unique in that 
it specifically enriches for peptides containing the relatively rare amino acid 
cysteine, thereby simplifying complex protein mixtures for subsequent MS 
analysis. The relative abundance of proteins can then be determined by 
monitoring the ratios of pairwise sets of selected peptide species which are offset 
by 8 amu. While representing a major advance, the ICAT approach is based on a 
sophisticated proprietary chemistry that analyzes relatively rare cysteine- 
containing peptides. 

Here, a complementary protein identification and quantitation strategy is 
described, which is termed Mass Coded Abundance lagging (MCAT), based on 
the differential post-experiment labeling of tryptic peptides with the lysine 
guanidation agent O-methylisourea followed by high throughput capillary liquid 
chromatography electrospray tandem mass spectrometry (LC-MS/MS). MCAT 
permits facile de novo sequencing of proteins present at pico- to femtomole 
levels in complex biological mixtures and provides for robust determination of the 
relative abundance of proteins in various cell states in a systematic, reproducible 
and straightforward manner. The development and applications of a systematic 
protein-expression-profilin^ 

should serve as a powerful means for characterizing the physiological, 
development or disease state of cells or organisms at the proteome level. 

Results 

De novo Peptide Sequencing using MCAT 

The MCAT sequencing method relies on the selective and quantitative (ie. 
complete) modification of the e-amine of C-terminal lysine residues of tryptic 
peptides with O-methylisourea (Fig. 1A). This reagent specifically and efficiently 
transforms lysine into homoarginine but does not react with the peptide amino 
terminus or other side groups (Kimmel, 1967). Peptide derivatization with O- 
methylisourea has previously been shown to facilitate peptide sequencing by 
MALDI post-source decay (Hale et al., 2000; Beardsley et al., 2000). Here, it is 
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shown that it can be used to sequence multiple individual peptides from, complex 
mixtures in a single high-throughput electrospray LC-MS/MS analysis. 

The MCAT de novo sequencing approach is based on two principles. First, a 
5 short sequence of contiguous amino acid sequence from a peptide (5-10 
residues) usually contains sufficient information to identify a corresponding 
unique protein. Second, peptides alternatively unmodified and modified with O- 
methylisourea differ by the mass differential encoded by the MCAT reagent (42 
amu). This allows the identities of the informative y-ion peaks to be readily 
10 delineated by comparing pair-wise sets of MS/MS spectra, allowing for 
systematic sequence determination. The MCAT labeling procedure is simple, 
economic and easy to perform with complex protein mixtures. 

The steps of the MCAT peptide sequencing procedure are as follows: (1) A 
15 protein mixture, which can be a purified polypeptide or protein complex, a cell 
fraction, or a crude cell extract, is first digested enzymatically with trypsin; (2) Half 
of the digest is derivatized to completion following incubation with an excess O- 
methylisourea; (3) The digests are desalted by C18 solid phase extraction and 
combined; (4) The pooled peptide mixture is fractionated by reverse phase HPLC 
20 and analyzed by automated ESI MS/MS. The mass spectrometer is operated in 
an automated dual mode whereby successive scans alternatively record a) the 
m/z of modified/unmodified peptide pairs as they elute from the column and b) 
the MS/MS fragmentation pattern of each peptide that has undergone collision- 
induced dissociation (CiD); (5) Following MS analysis, the data are processed to 
25 obtain the amino acid sequence identities of the components of the protein 
mixture. The process is illustrated schematically in Figure 1B. 

Inspection of pair-wise peptide spectra indicates that most ion peaks, notably the 
b-ion and y-ion series, are retained upon modification (Table 1). Since the C- 
30 terminal lysines of completely-processed tryptic digests are specifically labeled, 
the C-terminal y-ions produced during the MS/MS fragmentation reaction are 
mass shifted by the addition of the MCAT moiety. The y-ion peaks of the MCAT- 
modified peptides are offset by 42 amu (Fig. 2), or by factors of 42 resulting from 
the addition of a second or a third charge (ie. 21, 14 amu). In contrast, the 
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recorded m/z values for b-ions and chemical noise remain unchanged. Therefore 
comparison of MS/MS spectra for each unmodified/modified peptide pair allows 
ready determination of the y-ion peaks. With high quality spectra, discrimination 
of a well-defined and continuous y-ions series allows the amino acid sequence of 
a peptide to be readily deduced. This simplifies the spectral interpretation 
process, allowing for systematic sequence determination by assigning amino acid 
masses that correspond to y-ion peak distances using a reference table of 
monoisotopic amino acid masses. If required, a delta mass corresponding to a 
possible post-translational modification (e.g. +80.0 amu for phosphorylation on 
serine, threonine or tyrosine residues) or neutral loss (eg. water or ammonia) can . 
be incorporated into this table. 



In a systematic series of studies using a crude yeast cell extract (Table 1) it is 
established that MCAT provides an effective method for sequencing multiple 
peptides analyzed by LC-MS/MS. First, the ionization, charge and fragmentation 
properties of peptides were not greatly affected by the chemical derivatization 
procedure. Peptides generally have one of three different charge states (+1, +2f 
or +3), each of which results in a unique spectrum for the same peptide The 
spectra of numerous unmodified and modified peptide forms showed similar 
information content and could be correctly interpreted using database search 
algorithms with similar efficiency. Second, the modification of lysine-containing 
peptides occurred in- a robust, unbiased and reproducible manner. Third the 
mass tag (42 amu) added to the treated peptides was easily resolvable by MS 
regardless of charge state and did not overlap with other common adducts or 
peptide modifications. Even for a charge state of + 3. the delta mass is 14 units 
well within the resolution of a mass spectrometer. Fifth, the process simplified the 
spectral interpretation process so that the area of combinatorial sequence space 
to be searched was easily within the limits of modern computing technology. 

High confidence amino acid sequence was readily obtained for ten peptide 
spectra using the MCAT approach (Table 1). Good quality spectra were chosen 
from MS runs analyzing complex protein mixtures from various sources (a 
..bacterial cell lysate, a yeast cell lysate, and a human nuclear extract) Two 
representative analyses are shown in Fig. 2. The identifications were confirmed 
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using a computer database search algorithm. The SEQUEST algorithm (and 
similar algorithms) can detect MCAT modified lysine residues unequivocally 
because modification of a C-terminal lysine following trypsin digestion alters the 
m/z of y-series ions but not b-series ions relative to the unmodified peptide. 

5 

Although carried out manually here, the MCAT sequencing process may be 
formalized to facilitate automation. First, the mass of the tag (or a factor of it 
resulting from multiple charges) is added to each peak observed in the 
unmodified spectrum (above some threshold). The spectrum of the modified 

10 peptide is searched for peaks corresponding to these 'mass-tagged' peaks, any 
such peaks being candidate y-ions. Peaks appearing in both spectra are likely to 
represent b-ions or other ion products and are excluded from the initial analysis. 
Next, the mass differences between all candidate y-ions are calculated. Mass 
differences matching the known masses of single or double amino acids are 

is noted and attempts are made to extend the sequence from this starting point in 
both directions (i.e. higher and lower m/z) using known single or double amino 
acid masses. The putative sequences can be ranked using a score incorporating 
factors such as unbroken peak series and correlation of observed peaks with 
theoretical peaks. Moreover, for each putative y-ion series, the remaining peaks 

20 (L those conserved in the unmodified and modified spectra) are candidate b-ions 
and therefore can be used to impose further statistical limits on the y-ion 
designations. In other words, for any identified y-ion sequence ACDEFG, the 
corresponding sequence GFEDCA should be observed, and the extent of the 
presence or absence of the corresponding peaks can be factored into the overall 

25 score. 

Our results are typical of peptide MS/MS experiments in that incomplete y-ion 
series were generally observed. For high mass y-ions (yn, yn-1), this may occur 
because of charge repulsion; for low mass y-ions (y2, y3), because ion trap 
30 instruments generally fail to resolve ions lower than -1/3 the m/z of the precursor 
ion. Nonetheless, for most peptides examined, up to 8 to 15 continuous y-ions 
were detected, covering the bulk of the predicted amino acid sequence (Table 1). 
A properly ordered stretch of 6-7 amino acids is usually sufficiently informative to 
identify a corresponding protein using the BLAST algorithm. 
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Table 2 shows that MCAT reagent selectively modifies all lysine-terminated 
tryptic peptides present in the mixture in a quantitative and robust manner. In 
order to show that modification by the MCAT reagent is specific and that peptides 
so modified are recognizable by spectral identification algorithms, LC-MS/MS on 
a control yeast extract and a yeast lysate that had been treated with O- 
methylisourea was performed. The acquired MS/MS spectra were typically of 
high quality, with distinct b-series ion patterns the same for modified and 
unmodified spectra and the y-series offset by 42 Da, confirming that a C-terminal 
lysine had been modified (Fig . 2). Moreover, the SEQUEST scores for both 
modified and unmodified peptides were comparable and typical of high fidelity 
identifications. Importantly, in no case was an unmodified peptide detected in the 
treated sample (i.e. yielding high SEQUEST scores). The corollary was also true, 
with no peptides being significantly scored as being modified in an untreated 
15 sample (Table 2). 

Comprehensive LC-MS/MS analysis of an untreated and an O-methylisourea 
modified yeast cell lysate yielded significant SEQUEST scores for 291 peptides. 
For peptides treated with O-methylisourea, the rate of modification of non-lysine 

20 residues, such as arginine or alanine, by O-methylisourea was negligible (data 
not shown), as reported by others (Kimmel, 1967; Hale et al., 2000; Beardsley et 
al., 2000)7 Greater than 95% of SEQUEST-validated peptides containing lysine 
residues were classified as modified at lysine. In contrast, less than 3% of 
untreated peptides were scored as modified by SEQUEST, the same rate of 

25 false-positive scoring observed for arginine-containing peptides. These false- 
positives may result from poor quality spectra, or from acetylation or 
trimethylation of amino acids that generate a gain in mass (monoisotopic) of 
42.0106 Da or 42.0471 Da respectively. Such false positives can be easily 
eliminated upon inspection of MS/MS spectra because the y-ions series do not 

30 show the characteristic 42 amu shift 

Limitations to the MCAT sequencing method include the need for good quality 
spectra exhibiting a near continuous y-ion series. Furthermore, as with all de 
novo sequence efforts, some ambiguity remains due to the isobaric or near- 
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isobaric nature of certain amino acids (e.g. leucine and isoluecine). The MCAT 
approach is limited to peptides that terminate with a lysine residue. Tryptic 
fragments ending with arginine resdues are not modified and, therefore, cannot 
be sequenced by this approach. If necessary, endoproteinase LysC can be used 

5 instead of trypsin to generate peptides ending exclusively in lysine residues 
(apart from peptides derived from the C-terminus). Finally, it should be noted that 
incomplete trypsin or LysC digestion can potentially complicate the MCAT 
sequencing process by causing a mass shift in a subset of b-ions. However, the 
presence of modified internal lysine residues can be readily detected a priori by 

10 searching for parent ion mass shifts of multiples of 42 amu (adjusted for the 
charge on the ion). 

Relative Protein Abundance Determination Usin g MCAT 

The MCAT approach allows the relative abundance of proteins to be compared in 
15 two different samples following differential modification of peptides from one of 
the samples with O-methylisourea. By combining the peptides after treatment, the 
relative abundance of different protein species present in each sample can be 
estimated by measuring the signal intensities of the peptide pairs in a full scan 
MS analysis. The basic MCAT approach for measuring protein abundance is 
20 outlined in Figure 1C. 

In general, a first test sample and a second test sample may be an experimental 
sample (e.g. a sample exposed to a test compound of interest) and a control 
sample (not exposed to the test compound), respectively. Both samples are 
preferably enzymaticaiiy digested, for example in trypsin, and then one of the 
25 samples is treated (derivated) with a reagent to create a mass differential. This 
reagent may be called a mass differential reagent and is preferably a lysine 
guanidination compound. It may be, for example, o-methylisourea or any 
compound suitable for MCAT, that creates amino acids terminating in lysine or a 
homoarginine ending group or variant (memetic) thereof. The peptide of each 
30 test sample are then separated, for example ligand chromatography such HPLC, 
and subjected to MS. The MS spectra is obtained and the peptides in the first 
and second samples are identified, for example, by protein database searching. 
Optionally, the relative abundance of the peptides in the first sample and the 
second sample are determined, for example, by integrating the area under the 
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curve in a single ion intensity profile. Preferably, the peptide profile and relative 
abundance in the first and second sample is carried out in both orientations. 

MCAT protein quantitation is based on two principles: First, pairs of peptides 
alternatively unmodified and modified with O-methylisourea can be discriminated 
during a single MS run, thereby serving as mutual internal references for 
accurate relative quantitation. In MS, the ratios between the recorded signal 
intensities of the lower and upper mass components of these ion pairs provide a 
direct measure of the relative abundance of the two forms of a peptide and, by 
inference, the corresponding proteins in the. original cell pools. Second, the 
identity of the peptides can be obtained by performing MS/MS during the same 
analysis. 

The steps of the MCAT peptide quantitation procedure are as follows: (1) Two 
protein mixtures to be compared are obtained following different experimental 
treatment of a cell or tissue and are digested enzymatically with trypsin; (2) One 
digest is derivatized with O-methylisourea; (3) The peptides are desalted by C18 
solid phase extraction, combined, and the isolated peptides are separated and 
analyzed by automated multistage LC-MS/MS. The mass spectrometer is 
operated in a dual mode where two alternative scans cycle repeatedly. First, a full 
MS scan monitors the signal intensity of peptides eluting from the capillary 
column. Second, peptide sequence information is generated by selecting peptide 
ions for CID fragmentation in MS/MS mode. Sequence identification can be done 
usin 9 the de novo approach described above or using a protein database search 
algorithm. (4) Peptides are quantified by comparing the relative signal intensities 
of pairs of peptide ions with identical sequence that differ in mass due to lysine 
guanidination. In practice, an ion intensity profile is reconstructed for each 
sequenced peptide using the MS data and the relative abundance of modified 
and unmodified peptides calculated by integrating the area under the curve. The 
combination of MS and MS/MS data therefore determines the relative quantities 
and identities of the components of protein mixtures in a single analysis. The 
approach is illustrated schematically in Figure 1C. 
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The MCAT approach serves as an effective method for determining relative 
abundance of proteins by LC-MS/MS since: (1) O-methylisourea derivatizes all 
lysine-containing peptides present in the mixture in a quantitative manner; (2) the 
agent adds a mass tag to the treated peptide that is easily resolvable by the 
5 mass spectrometer and that does not overlap with common adducts or peptide 
modifications; (3) the modification preserves the charge and ionization properties 
of peptides such that the efficiency of ionization and signal intensity are 
equivalent; and (4) the modified peptides generally co-elute during standard 
reverse phase chromatographic separation. 

10 To illustrate the process, the relative abundance determination of the peptide 
LPWFDGMLEADEAYFK (SEQ ID NO:15) from two replicate yeast whole cell 
extract experiments is shown in Figure 3. Base peak chromatograms show many 
peptides eluting over a 60min run, while selected ion tracings for the predicted 
doubly-charged unmodified and modified forms of the peptide show both eluting 

is at 35-36min (Fig. 3A). A single full scan of an ion trap mass spectrometer 
operated in MS mode is shown in Figure 3B. Two prominent ion species are 
discemable and indicated with respective m/z values 21 m/z units apart (Fig. 3B). 
The fact that the ions co-elute, have a detected mass difference of 21 m/z units, 
and have identical sequences (data not shown) identifies them as a pair of doubly 

20 charged sister peptides. Over the course of the 60 minute elution gradient, more 
than 2,000 MS scans were automatically acquired. Figure 3C shows 
reconstructed ion chromatograms for each of the peptide species. The relative 
quantities were determined by integrating the curves contouring the respective 
eluting peaks. The ratio (unmodified:modified) was determined as 0.88 (Table 2). 

25 The peaks in the reconstructed ion chromatograms appear serrated because the 
MS system alternates between MS and MS/MS modes in order to both measure 
ion intensity as well as generate a mass spectrum of selected peptide ions for the 
purpose of protein identification. 

30 Table 2 shows some representative high-scoring peptides from a representative 
MCAT LC-MS/MS analysis of a yeast cell extract. In these experiments a 1:1 
mixture of unmodified:modified peptides was analyzed, and single ion tracings for 
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select peptides throughout an entire chromatographic run typically showed 
isolated peaks with the unmodified form co-eluting, or eluting slightly earlier, than 
the modified form (Fig. 3A and C). For nearly all peptides examined, the ratio of 
unmodified to modified signal was close to the expected 1:1. The range of signal 
intensities were generally within two-fold of the unmodified form and the 
percentage error (the difference between the observed and expected 
abundances) ranged from 1 to 62% (Table 2). Some exceptions were evident and 
excluded from the analysis. These included peptides that could be positively 
identified but whose signal is very weak, and peptides containing arginines that 
were modified in addition to lysine at low frequency. Another category of ion. 
found unsuitable for quantitation were singly-charged ions. It is unclear why this 
is the case but the signal from singly-charged ions is typically lower than that for 
doubly- or triply-charged ions, possibly rendering them less likely surpass the 
intensity threshold required for accurate quantitation. 

Figure 4 shows variation in the measured relative abundance for two peptides 
from the same parent protein (and therefore are present in equimolar 
concentrations) in three replicate experiments. Importantly, multiple peptides 
independently analyzed for several proteins gave similar linear responses. 
Experiment-to-experiment variation for these peptides is within 25% and variation 
within a single experiment for peptides derived from the same protein is within 
20%. The variation from ideal 1:1 ratio is not the result of reduced ionization 
efficiency or MS signal of the modified peptides relative to their unmodified forms 
because the effect was consistently observed in subsequent experiments 
independently of which sample was chosen for modification. More likely, it results 
from modest variations in peptide recovery during sample workup. 

In order to correct for any possible systemic labeling errors, for instance 
preferential labeling by O-methylisourea of one sample, MCAT quantitation can 
be carried out in reciprocal orientations. For this reason, when comparing two 
independent protein samples (A and B), derived for instance from two distinct cell 
states, the basic MCAT procedure can be carried out in four complementary and 
reciprocal mass spectrometry analyses: I) unmodified sample A versus modified 
sample B; II) unmodified sample B versus modified sample A; III) unmodified 
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sample A versus modified sample A; IV) unmodified sample B versus modified 
sample B. The ratios of unmodified to modified peptide signals obtained in 
experiments III and IV can be used to systematically normalize and control for 
variations in the data obtained in experiments I and II, respectively. In practice, 
5 the MCAT analysis can be simplified into a two-tiered reciprocal experiment set, 1 
and II, which should independently confirm any significant quantitative 
observations obtained in a sample comparison. 

To confirm the quantitative nature of the MCAT approach, mixtures of modified 
10 and unmodified peptides derived from a common crude yeast cell extract were 
prepared at various ratios and analyzed by a 30 minute LC-MS/MS analysis. The 
MS/MS spectra acquired were used to search a non-redundant genome 
database using the SEQUEST algorithm (Eng et al., 1994) to identify the proteins 
present in mixtures. The relative ratios of 5 peptide sister pairs was quantified as 
is described above (Fig. 4B). This analysis shows the relative abundance of 
proteins can be accurately determined (i.e. exhibits a linear response) over a >30 
fold dilution series. Beyond this range, the weaker signal of the two abundances 
was indistinguishable from background noise in these experiments. 

20 It should be emphasized that the data were acquired for polypeptides present at 
a pico- to femtomole level in a highly complex protein mixture. The loading 
capacity of capillary reverse phase columns for complex peptide mixtures 
imposes a strict limit on the detection of low abundance proteins by LC-MS/MS. 
With a purified protein, most current MS systems generally exhibit a practical 

25 dynamic range of roughly three orders of magnitude based on maximal signal to 
noise ratios that can be acquired (using a purified or low complexity protein 
preparation). However, sophisticated chromatographic separation techniques can 
be coupled to fractionate complex peptide mixtures prior to MS in order to 
substantially improve the detection limits of MS protein analysis (Link et al., 1999; 

30 Washburn et al., 2001). Hence, when combined with the MCAT approach, 
determination of the relative abundance of moderate to low abundance proteins 
should be achievable even in the absence of enrichment. 
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An experimental approach for systematically sequencing and quantifying proteins 
isolated from complex biological mixtures using basic chemistry and mass 
spectrometry techniques is described and validated. De novo sequencing 
expands the range of organisms that can be analyzed and removes the reliance 
on DNA sequence databases that may be incomplete, erroneous- or that fail to 
account for complexities introduced by alternative splicing, protein modifications, 
or protein polymorphism. The quantitative capabilities of the method also 
overcome a significant limitation of current proteomics technologies, whereby the 
determination of protein abundance on a large-scale is generally low throughput, 
expensive, and tedious, for instance, radiolabelling of proteins before analysis by 
two-dimensional gel electrophoresis and quantitation following isolation of 
individual spots (that may contain one or more polypeptides). 

The ICAT method reported by Aebersold and coworkers (Gygi et al., 1999) may 
significantly improve throughput and reduce sample complexity by enriching for 
proteins containing the underrepresented amino acid cysteine. These features 
are useful for sampling a mixture whose proteome complexity could overwhelm 
the ability of current LC-MS technology to resolve it. The MCAT strategy 
described here is not limited to any particular affinity chemistry and in principle 
can be coupled to analogous affinity-based enrichment steps. For this reason, 
MCAT can potentially be used to identify and quantify all the proteins present in a 
biological sample. In combination with powerful multi-dimensional LC protein 
separation techniques, such as that described by Yates and coworkers (Link et 
al., 1999; Washburn et al., 2001), considerable depth in proteome coverage may 
be achieved. Quantitative data describing patterns of peptide or protein 
expression for many hundreds or thousands of proteins can be used to identify or 
classify protein 'profiles' in a similar manner to that routinely used for gene 
expression data. The combined MCAT approach can therefore be used for 
identifying, classifying and characterizing functions of known and unknown gene 
products, for characterizing metabolic and other functional protein pathways in 
cells, and for identifying proteins and pathways targeted by drugs and other 
reagents. 



( 



The MCAT method offers key experimental gdvantaaes . 
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First, the approach is simple and effective. It builds on established MS techniques 
and principles that are flexible and can easily be adjusted for large-scale projects, 
including efforts to generate peptide or protein profiles describing the effects of 
5 environment, mutation, disease or experimental interventions such as drug 
treatment. Significant patterns of expression can be identified with appropriate 
software and data mining algorithms. 

Variations of the MCAT approach can easily be devised, including strategies to 
10 address other quantitative aspects of protein expression, those searching for 
post-translational modifications, or those screening for mutant proteins. It is likely 
that the number of unique peptide species per organism will be multiplied 
significantly by the presence of post-translational modifications compared to 
genome predictions. Because the mass of many common important modifying 
15 groups are known, and because their preferences for particular amino acids are 
often known, the database can be searched for ions predicted to result from 
peptides with specific modifications. 

Finally, the addition of a dynamic component to the molecular descriptions of 
20 protein activities is likely to prove critical to our understanding of the biochemical 
circuitry within cells. Consequently, the development of robust analytical 
methods, such as the MCAT approach described here, that allow for efficient 
identification and quantitation of large numbers of proteins from complex mixtures 
can be expected to have a major impact. 
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Experimental protocols 

Materials. Media, standard-grade and HPLC-grade laboratory chemicals were 
5 obtained from Fischer Scientific (Fair Lawn, NJ). O-methylisourea (S- 
methylisothiourea hemisulfate salt) was from Sigma-Alderich (St. Louis, MO). 
Poroszyme immobilized trypsin was from Applied Biosystems (Framingham, MA). 

Preparation of protein extracts. The protease-deficient S. cerevisiae yeast 
10 strain BJ5460 was grown to late-log phase (OD ~3) at 30°C and protein whole 

cell extracts prepared as follows: Cells were harvested, frozen, and mechanically ( 
lyzed by grinding in the presence of dry ice. The cells were thawed in lysis buffer 
(8M urea, 1 mM CaCI 2 , 100 mM Tris-HCL, pH8.5). Insoluble debris was pelleted 
by a high-speed (20 K x g) spin and the supernatant diluted to 2M urea using 
15 digestion buffer (100 mM Ammmonium bicarbonate, pH8.5, 1 mM CaCI2. A 
bacterial whole cell extract was similarly prepared using the E. coli DH5a strain. 
Human nuclear extracts were prepared using a commercial kit (Pierce), and 
diluted into digestion buffer. 



20 Tryptic Digestion and Peptide Derivatization. Porozyme immobilized trypsin 

beads were added to an aliquot of each protein extract at a 1:500 protein ratio ( 
and the digests incubated at 30°C for two days with tumbling. The extracts were 
aliquoted into two microtubes. Solid O-methylisourea was added to one of the 
tubes to achieve a final concentration of 1M. Base (NaOH) was added to 0.5N to 

25 adjust the pH to >10. The reaction was incubated at 37°C overnight. The peptide 
mixtures were extracted by solid-phase extraction using SPEC-PLUS PTC 18 
cartridges (Ansys Diagnostics, Lake Forest, CA) according to the manufacturer's 
instructions and buffer exchanged into a 5% ACN, 0.1% formic acid solution. 
Samples not immediately analyzed were stored at -80°C. 
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MCAT peptide sequencing. Each sample was subjected to microcapillary LC- 
MS/MS analysis with modifications to the general method described by Link and 
coworkers (1999). A quaternary Surveyor HPLC pump (ThermoFinnigan Canada) 
was directly coupled to a Finnigan LCQ-DECA ion trap mass spectrometer 

5 equipped with a custom microLC electrospray ionization source. A fused-silica 
microcapillary column (100 pm i.d. x 365 pm i.d.) was pulled with a Model P-2000 
laser puller (Sutter Instrument Co., Novato, CA) as described. The microcolumn 
was packed with 10 cm of 5 pm Ci 8 reverse-phase material (Zorbax XDB-C18, 
Hewlett-Packard). Approximately 100 pg of the unmodified fraction and 100 pg of 

10 the derivatized peptide fraction were combined and loaded onto a single 

microcolumn for sequence analysis. After loading, the column was placed in-line 
with the ion source system setup as described (Link et ai, 1999). A fully 
automated 30 min 100% buffer A (5% ACN, 0.1% formic acid) to 80% solvent B 
(95% ACN, 0.1 % formic acid) binary gradient was run at a flow rate of -0.3 

15 ul/min. Eluted peptides were analyzed by automated MS/MS as described by 
Link and coworkers (1999) except that a full scan range of 400-1600 m/z was 
used. 



SEQUEST analysis. The SEQUEST algorithm (Eng et al., 1994) was run on 
20 each dat set against sequence databases obtained from the National Center for 
Biotechnology Information (Bethesda, MD). Positive sequence identification was 
based on several criteria (XCorr and DCn score, and the presence of tryptic 
termini) described at http, and all identifications were confirmed manually. 

MCAT protein quantitation. Pairs of samples to be compared were subjected to 
automated uLC-MS/MS analysis with modifications to the general method 
described above. Approximately 200 pg of the unmodified fraction and 200 pg of 
the derivatized peptide fraction were combined and loaded onto a microcolumn. 
After loading, a fully automated 30 or 60 min 0-80% A:B gradient chromatography 
run was carried out on each sample. The buffer solutions used for the 
chromatography were 5% ACN/0.1% Formic acid (buffer A), 80% ACN/0.1% 
Formic acid (buffer B). Eluting peptides were analyzed by coupled automated 
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uLC-MS-MS/MS techniques as described above. There was a consistent slight temporal 
difference in the elution of unmodified/modified peptide pairs, with the unmodified light 
analog eluting slightly before the heavy form. Selected ion traces for each peptide pair 
were quantified using the ADDXPRESS program by which the peak area of each eluting 
peptide was reconstructed and used in the ratio calculation. 



Table 1. De novo peptide sequencing from complex mixtures using MCAT 



Identified 
peptide 


b-ion s 


eries* 




b*-io 


i series" 


V-ion series" 


V*-ion series" 




E Predicted 
AA° 


tu 

1 

(A 

A 


5 

£ S 

717.8 ' 


•o 

I s 

O 


JC 
o 
<s 


tu 

7173 


0) 

Si 

■c 

o 


D 

JZ 

5 




■a 

u 

•c 

S.E 

X 
IU 

7483 


9 

o 

7483 


1 
S 
A 


is 

! ^ 

7903 


"a 
S> 

si 

jo fe 
o 

7913 


a 
O 

To 


V 

42.2 


+ 
>> 

137.0 


Yeast 
YGR912C 

VINDAFGIE 
cGLMTTVHS 
LTATQK 
(SEQ. ID. 
NO: 16) 
m = 2575.9 
2=2 


831.0 
960.1 


831.6 


& 


831.0 
960.1 


831.6 


A 


0.0 


886.0 
985.1 


686.3 
985.4 


«A 

A 


928.0 
1027.1 


928.0 
1027.7 


A 
A 


41.7 
423 


99.7 
101.1 


V 
T 


A 


1089.2 






1089. 
2 


1089.2 


A 




1086.2 


1086 
4 


A 


1128.2 


11283 


A 


42.4 


100.5 


T 


A 


1146.2 






1146. 
2 


1146.2 


A 




11873 


1 1 BTt 
J.J.O/. 

6 


■A 


1229.3 


12293 


A 


41.7 


1313 


M 


A 


1259.4 






1259. 
4 








13183 


1318. 
3 


«* 


13603 


1360.6 


.A 


423 


1133 




A 


1390.6 






1390. 
6 








1431.7 


1431. 
7 




1473.7 


1473.9 


A 


42.2 


57.2 


G 


A 


1491.7 


1491.9 


A 


1491. 
7 


14913 


A 


0.1 


1488.7 


1489. 
0 




1530.7 


1531.1 


A 


42.1 


129.0 


E 


A 


1592.8 






1592. 
8 








16173 


1617. 
9 




1659.8 


1660.1 


A 


42.2 


129.2 


E 




1691.9 






1691. 
9 








1747.0 


1747. 
4 


.A 


1789.0 


17893 


A 


41.9 








1829.1 






1829. 
1 


1829.1 


A 




1860.1 






1902.1 














1916.1 
340.5 
453.6 
567.7 


19163 
340.5 
453.6 
5673 


A 
A 
A 


1916. 
1 

3403 
453.6 
567.7 


^19163 
340.5 
453.5 
5673 


A 

—So* 

A 


0.0 
0 0 
0.1 
a n 


1917.2 

317.4 

4313 


1917. 
3 


A 


1959.2 

359.4 

4733 


1959.4 


A 


42.1 








RBSB 


664.9 
766.0 
881.1 
968.1 


766.2 


A 
A 


664.9 
766.0 
881.1 
968.1 


665.4 
880.7 


-A 
A 
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658.8 
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658.2 
773.6 


tA 
■A 

A 
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629.7 
7003 
8153 
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629.4 


-A 
A 


40.9 
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99.1 


V 




IUJNPTD5D 

AVGNAVK 

(btvj. ID. 

NO:17) 

m = 1740 0 

2=2 


1083 .2 






1083. 
2 








860.9 
976.0 


975.4 


A 


902.9 
1018 


903.5 
1018.4 


A 
A 


43 


114.9 


D 


A 


11543 


11543 


A 


3 








1077.1 


^lO??^ 
5 




1119.1 


1119.6 




42.1 


101.2 


T 


A 


1253.4 


12 c 533 


A 


1253. 
4 


12533 


A 


0.2 


1174.2 


1174. 
5 


.A 


1216.2 


12163 


A 


42.0 


96.9 


P 


-A 


1310.5 






1310. 
5 








12883 


1288. 
5 


A 


13303 


13303 


A 


42.0 ' 


114.0 


N 


A 


1424.6- 


1424.6 


A 


1424. 
6 


1424.0 




0.6 


1401.5 


1401. 
6 




14433 


14433 


A 


41.9 


113.0 


I i 


A 




1495.7 






1495. 
7 








1514.7 


1514. 
1 




1556.7 
















1594.8 

526.6 
663.7 


1594.6 
663.4 


A 
A 


1594. 
8 

526.6 
663.7 


1594.6 
663.4 


A 


0.0 


1627.8 

568.7 
639.7 


5683 
639.4 


A 
..A 


1669.8 

610.7 
681.7 


610.7 
681.7 


A 
A 


42.4 
423 


71.0 


P 
A 




HumanACT 
B 


7603 
859.9 

973.1 


760.8 


A 


7603 
859.9 

973.1 


859.6 
972.5 


A 
A 




768.9 
870.0 

983.1 i 
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869.4 

983.4 


A 
A 
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912.0 

1025.1 
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911.5 

1025.1 


A 
A 

-A 


41.9 
42.1 

41.7 


1283 
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E 
T 


A 

A 
A 

A 
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A 
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3 
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A 


0.2 
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5 


.A 
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-A 


43.1 
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I 




LTEAPLNPK 


1187.4 
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4 
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A 
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A 
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NO:18) 
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5 ■ 
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6 


A 
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A 


41.9 
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1387.4 
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6 
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A 


0.1 
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7 


.A 
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-A 
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A 
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A 
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7 
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A 
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a. b and b* refer to unmodified and modified b-ion series respectively 

b. y and y* refer to unmodified and modified y-ion series respectively 

c. g£ indicates a match between expected and observed m/z values (tolerance of 2.0 m/z units) 

d. Ab, Difference between observed b and b* m/z values 

e. Ay, Difference between observed y and y* m/z values 

f. A(y,y+1), Difference in observed m/z between successive y series ions, adjusted for charge state of ion 

g. Predicted AA, Amino acid residue predicted using A(y,y+1) 

h. s& indicates a match between MCAT-predicted and SEQUEST-predicted amino acid. 



Table 2. Identification and quantitation of peptides from a yeast whole cell digest 















Identification* 


Quantitation 0 


Protein 


Peptide 




z 


m/z 6 


Soore c 


-MCAT 


+MCAT 


Measured 
abundance 


% 
error 














P 


P* 


P 


p 

* 


P 


p* 




YBR118W 


SVEMHHEQLEQGVPGDNV 
GFNVKfSEQ. ID. NO:19) 


2550.8/ 
2592.8 


2 


1276.4/ 
1297.4 


2.2433/ 
2.5321 




X 


X 




1.00 


0.76 


24±4 




Tl 1 FATDATFOP^RPTDKPL 
RLPLQDVYK# (SEQ. ID. 
NO"20) 


3320.8/ 
3404.8 


3 


1107.9/ 
1135.9 


33888/ 
3.3370 




X 


X 




1.00 


0.63 


37±5 




VETGVIKPGMWTFAPAGV 
TTEVX* (SEQ. ID. 
NO:21) 


2430.9/ 
2472.9 


2 


1216.4/ 
1237.4 


2.5458/ 
2.1831 


Eft 


X 


X 


en 


1.00 


0.38 


62 ±12 


YCR012W 


ALENPTRPFLAILGGAK 
(SEQ ID NO*22) 


1768.1/ 
1810.1 


2 


885.0/ 
906.0 


1.7773/ 
1.4083 




x 


x 


e£ 


1.00 


0.57 


43 ±5 


YDFU55C 


HWFGEWDGYDIVK j 
(SEQ. ID. NO:23) ! 


1675.9/ 
1717.9 


2 


838.9/ 
859.9 


3.7988/ 
3.6211 






x 


s£ 


1.00 


0.71 


29 ± 5 


YDR487C 


HGIPUSIEEIAQYLK 
fCi=n TD 


1824.2/ 
1866.2 


2 


913.1/ 
934.1 


2.1238/ 
1.6387 






x 




1.00 


0.86 


14 ± 1 


YGR063C 


LPAEWELLPHYKPR 
(SEQ. ID. NO: 25) 


1761.1/ 
1803.1 


2 


881.5/ 
902.5 


2.0444/ 
1.9739 




x 


x 




1.00 


0.66 


34 ± 6 


YGR192C 


INDAFGIEEGLMTTVHSLT 
ATQK 

( SEO. ID. NO:26) 


2476.8/ 
2518.8 


2 


1239.4/ 
1260.4 


2.9164/ 
4.1100 




X 


X 




1.00 


0.52 


48 ±28 




V1NDAFGIEEGLMTTVHSL 
TATQK 

(SEO. ID. NO:27) 


2575.9/ 
2617.9 


2 


1288.9/ 
1309.9 


3.1456/ 
3.3717 




X 


X 




1.00 


0.44 


56 ±17 




VPTVDVSWDLTVK 

(SEQ. ID. NO:28) 


1512.7/ 
1554.7 


2 


757.3/ 
778.3 


3.2279/ 
3.1548 




X 


X 


A 


1.00 


1.29 


29±11 


YGR214W 


NVQVHQEPYVFNARPDGV 
HVINVGK (SEQ. ID. NO: 
29} 


2817.2/ 
2859.2 


3 


940.0/ 
954.0 


1.8494/ 
2.2204 




X 


X 




1.00 


0.61 


39 ±10 


YGR254W 


AQYNHQGWDHLSLLPTF 
GAK (SEQ. ID. NO:30) 


2388.7/ 
2430.7 


2 


1195.3/ 
1216.3 


2.4748/ 
3.0844 




X 


X 


eA 


1.00 


0.81 


19 ±2 




YPIV5IEDPFAEDDWEAW 
SHFFK 

(SEO. ID. NO:3n 


2829.1/ 
2871.1 


3 


944.0/ 
958.0 


3.1108/ 
3.2183 




X 


X 




1.00 


0.61 


39 ±9 


YHR174W 


WLTGVELADMYHSLMK 
(SEQ. ID. NO:32 


1894.2/ 
1936.2 


2 


948.1/ 
969.1 


4.0552/ 
3.8246 


i* 


X 


X 


E* 


1.00 


0.77 


23 ± 3 


Y3R105C 


TVIFTHGVEPTWVSSK 
(SEQ. ID. NO:33) 


1800.1/ 
1842.1 


2 


901.0/ 
922.0 


1.5600/ 
1.8810 




X 


X 




1.00 


0.75 


25±4 


YKL060C 


SPIILQTSNGGAAYFAGK 
(SEQ. ID. NO:34) 


1795.0/ 
1837.0 


2 


898.5/ 
919.5 


3.6709/ 
4.2032 




X 


X 


E^ 


1.00 


0.73 


27±5 




TGVIVGEDVH NLFTYAK 
(SEQ. ID. NO:35) 


1863.1/ 
1905.1 


2 


932.5/ 
953.5 


3.2735/ 
2.6813 




X 


X 




1.00 


0.75 


25±4 


YLR044C 


KLIDLTQFPAFVTPMGK# 
(SEQ. ID. NO:36) 


19063/ 
1948.3 


2 


954.1/ 
975.1 


3.5845/ 
3.9361 




X 


X 




1:00 


0.83 


17±2 


YLR058C 


EVIYDLENPINFSVFPGHQ 
GGPHNHTIAALATALK 
(SEQ. ID. NO:37) 


3772.2/ 
3814.2 


3 


1258.4/ 
1272.4 


1.8356/ 
2.5693 




X 


X 


E* 


1.00 


0.73 


27 ±6 



a. Molecular mass of unmodified/modified peptides ions. 

b. Mass-to-charge ratio of unmodified/modified peptides. 

c. SEQUEST cross-correlation score for unmodified/modified peptide. 

d. Identifications were determined in untreated samples (-MCAT) or samples modified using MCAT (+MCAT). ^ or x 
indicates that the unmodified (P) or modified (P*) peptides were observed (b£) or not observed (x) in the respective 
sample. 

e. Relative abundance measurements are for 1:1 mixtures of unmodified and modified samples. Percentage error refers to 
deviation from ideal (1:1) ratio ± standard deviation for multiple measurements. 

# These peptides were modified at more than one lysine residue. 
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Further discussion of the figures related to MOAT 

(1) The MCAT approach for peptide sequencing and relative protein 
abundance determination. 



See Figure 1. (A) The guanidination reaction is specific for the side chains of lysine, 
which is selectively converted to homoarginine. (B) For sequencing using MCAT, 
protein mixtures are first digested with trypsin, which generates peptides suitable for MS 
analysis that terminate with lysine or arginine residues. Half of the sample is treated 
with the MCAT reagent O-methylisourea. Peptides ending in lysine are modified, which 
adds 42 amu to the mass of the peptide but does not alter the properties of the peptide 
during LC-MS analysis. The peptides mixtures are combined at a 1:1 ratio, separated 
by reverse phase LC and introduced online into a MS instrument using electrospray 
ionization. Following tandem MS analysis, peptide sequence is determined by 
comparing MS/MS spectra of unmodified and modified peptides. The fragmentation 
pattern of both sister peptide pairs are similar except for the shifted y-ion series, which 
can be deconvoluted to reveal the amino acid sequence of the peptide. (C) For relative 
abundance measurements, samples representing different cell states are alternatively 
modifie d or unmodified with MCAT. Full MS spectra are recorded for sister peptide 
species and their relative abundance determined by measuring the respective trace 
intensities on reconstructed single ion chromatograms. 
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(2) MCAT enables identification and quantitation of complex protein 
mixtures. 

5 See Figure 2. (A) Ion chromatograms recorded for the base peak (top), an 
unmodified peptide ion [LPWFDGMLEADEAYFK+2H] +2 (middle) and its 
corresponding 0-methylisourea(MCAT)-modified form (bottom). When mixtures 
of untreated and MCAT-treated protein digests are resolved by reverse phase 
LC, the modified peptides elute with a minor delay compared to the respective 

10 unmodified forms (35.9 vs. 35.7 min respectively in this example). (B) Depending 
on charge and the number lysine residues, the m/z signals observed for pairs of 
unmodified or modified peptide ions during MS are offset by 42, 21 or 14 m/z 
units (for plus 1 , 2 or 3 ions respectively). In this example, the peak signals 
recorded for the unmodified (967.07 m/z) and modified (988.08 m/z) forms of the 

15 peptide are offset by 21 m/z units, indicating a +2 charge. The peptide ions are 
then independently selected and automatically fragmented by MS/MS. 
Comparison of the y-ion series allows the amino acid sequence to be determined. 
(C) The relative abundance of individual peptides can be determined by 
reconstructing the chromatograms for the unmodified and modified forms of the 

20 peptide ions and calculating the ratio of signal intensities using area under curve 
integration. 

(3) be novo sequencing of a yeast peptide and a human peptide using 
MCAT approach. 

25 

See Figures 3A and 3B. (A) The peptide WDLVEHVAK (SEQ ID 
NO:38)analyzed by MCAT LC-MS/MS in a digest of yeast whole cell extract. A 
representative MS/MS spectrum of the unmodified peptide (top) and the 
corresponding spectrum for the modified form (below) are shown. Because the 
30 MCAT reagent reacts specifically with lysine residues, the carboxy-terminal lysine 
of a tryptic peptide is uniquely modified. Therefore, the signals for the y-series of 
ions (where charge localizes to the carboxy-terminal lysine) are shifted +42 m/z 
units and can be immediately identified, whereas the b-series of ions (where 
charge is retained at the amino terminus) are unaltered. The expected m/z values 
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for b- and y-series ions of the unmodified and modified peptides are given (right), 
with those observed in the experiment underlined. The amino acid order is 
resolved by measuring the mass difference between successive y-ion peaks. (B) 
The peptide VAPEEH PVLLTEAPLN PK (SEQ ID NO:39)was identified in a digest 
of nuclear extract from HeLa cells. In this peptide a stretch of ten amino acids (A- 
E-T-L/l-L/l-V-P-H-E-E) can be identified by mapping y-ions to the bands shifted 
by 42 m/z units in the modified spectrum (bottom) relative to the unmodified 
spectrum (top). The dominant peak at 892.9 in the unmodified spectrum is 
approximately 21 m/z units from an dominant unassigned peak at 914,4 in the 
modified spectrum. These peaks probably represent doubly-charged y16 ions 
that terminate in with proline, an amino acid commonly observed to form 
dominant peaks during CID. The other major peak in both spectra (1292.6 and 
1334.5 in the upper and lower panels respectively) is a singly-charged y12 ion 
that also terminates wtih proline. Therefore, an additional advantage of the MCAT 
technique is the resolution of such ambiguous peaks through charge 
determination. In the case of both yeast and human peptides, the identical 
molecular masses of leucine and isoleucine prevent their resolution by MS. 

(4) The MCAT method is reproducible and quantitative. 



See Figures 4 A and 4 B. (A) A yeast wh ole cell was digested with trypsin in three 
replicate experiments (A, B, C). Each digest was divided into two equal portions, 
one of which was treated with O-methylisourea. Each pair of mixtures was then 
recombined at a 1:1 ratio and protein quantitation determined by the MCAT LC- 
25 MS/MS. The relative abundance ratios (expressed at the ratio of modified to 

unmodified peptide signal) of a subset of positively-identified peptides is given for 
each analysis. (B) Untreated and MCAT-labeled yeast protein tryptic digests 
were combined in varying proportions ranging from from 16:1 (modified to 
unmodified) to 1:16 effective concentrations. The measured relative abundance 
ratios for five representative peptides are plotted versus the log(10) of the dilution 
ratio. 
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Peptide Profiling 

Below examples are shown of the utility of peptide profiling as a means to 
characterize and classify diverse human tissues, to characterize subcellular 
fractions of individual tissues, and to illustrate how a database of such peptide 
5 profiles can serve as a depository of protein expression information that can be 
mined rapidly and accurately for knowledge about the status of an unknown 
sample. This process is robust, sensitive and reproducible. Although the method 
is generally applicable, the following serve to illustrate select uses of the 
approach. 

10 

Example 3: Use of peptide profiles to characterize human tissue 
The invention includes methods of characterizing human tissue. The method 
comprises generating samples suitable for MS analysis and producing a peptide 
profile. The relative abundance of peptides in samples is also preferably 
15 determined. The peptide profile that is generated is compared to peptide profiles 
in a database or library using common algorithms in order to identify cognate 
proteins, preferably those that are considered important therapeutic targets, as 
well as metabolic enzymes and structural proteins. 

20 Table 1 shows 40 peptides sequenced and quantified from a human lung tissue 
lysate sample in a single LC-MS analysis that are then used to construct a unique 
peptide profile. The peptides in turn allowed for the identification of cognate 
corresponding proteins present in the sample (a total of 867 proteins were 
unambiguously identified in this analysis). Note that the peptides sequences 

25 obtained by a generic database search algorithm were both preceded by, and 
terminated with, a K or R residue as a result of cleavage of the input proteins by 
trypsin. The sequence of a total of 1896 peptides were determined in this one 
analysis with high accuracy and sensitivity, demonstrating the ability of the 
approach to generate a detailed profile or fingerprint of protein expression of a 

30 complex tissue. 

Table. 1. Partial List of Peptides observed in human lung tissue used for peptide 
profiling, 

K.AAIANLCIGDLITAIDGEDTSSMTHLEAQNK.I 

(SEQ. ID. NO:40) K.AGNNMLLVGVHGPR.T (SEQ. ID. NO:59) 
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I^tn ^™' DHWPE PGTSLLAAFDQWR.E ICAHGPGLEGGLVGKPAEFTIDTKG 
KAAPLSLCALTAVDQSVLLLKPEAK.L 

!?f ?J D " NO:42> KAHSPQGEGEIPLHR.G (SEQ ID NO"61) 

^SISI^^ 0 ' 43 ' ^"VSFKPTVAOOR, (S ( EQ . ,0. NO^ 

if f ?J D ' N ° :44) KAIEVIRPAHILQEK.E (SEQ. ID N0 63) 

KAASWLTHQGSFHGAFR.S (SEQ. ID. N0.45) KAIQDAGCQVLK.C (SEQ ID N0 64) 

K.ACISPKPQKPWDK.D (SEQ. ID. NO:49) KALESPERPFLAILGGAK.V (SEQ ID NO 68) 

ICADIIYPGHGPVIHNAEAKI (SEQ. ID. NO:50) ^%^^^^L^E 
ICAEEVAFWTELLAK. N (SEQ. ID. NO:51) KALHASGAK.V (SEQ. ID NO70V 

^p?m^^L L ^ (SEQ - ,D - NO:52) K.ALHASGAKWAVTR.T (SEQ. ID NQ 71) 

ICAFAMIIDKLEEDISSSMTNSTAASRPPVTLR.L J a 

Jf E °: IPrJUSSS 1 ICALLNNSHYYHMAHGK.D (SEQ. ID N072) 

K.AFAQAQSHIFIEK.T (SEQ. ID. NO:54) KALN RPPTYPTK Y (SEQ ID NO-73) 

K.AFISNVKTALAATNPAVR.T (SEQ. ID. NQ:55) UoTm^^^ 

SSS^"^* 8 "* J^DH^GTLLKPNMVTPGHACTQ^F 
K.AGAPPGLFNWQGGAATGQFLCHHR E * 

!?A E GHP^ KALTGGIAHLFKQ (SEQ. ID. NO:76) 

(SEQ. ID. NO:58) KALVKPQAIKPK.M (SEQ. ID. NO:77) 

A further embodiment of the invention includes using profiles such as this to 
compare different tissues or experimental samples. For instance, a comparison 
5 of the peptide profiles for human pancreatic and heart tissues can be made with a 

simple 2-dimensional plot that can be extended to 'n' different planes as required C 
(for 'n' types of tissue, samples, or patients). Comparison of the peptide profiles 
of these samples can be done using standard computational methods (e.g. 
agglomerative clustering). In the case of human pancreatic tissue, the analysis 
10 showed that although several proteins are shared between the tissues, many are 
not. Therefore, a further embodiment of the invention is the use of peptide 
profiles to characterize tissues and thereby categorize samples. 

Although this patent describes primarily approaches involving peptide profiling, 
is the approach can be extended to whole protein profiling (and to other 

applications where separation techniques compatible with mass spectrometry 
may be used to elicit a profile, for instance lipid profiling, phosphoproteins 
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profiling, small molecule metabolite profiling; these methods preferably involve 
tagging the compounds of interest and performing LC-MCAT to generate a lipid 
profile, phosphoprotein profile, small molecule metabolite profile. The methods 
can provide identity and relative abundance information by readily adapting the 
5 methods described herein with peptides.). 

Table 2 shows some of the corresponding proteins (of the 867 unique proteins 
identified in this analysis) identified by searching the SwissProt Protein database 
using the identified peptide sequences (http://www.expasy.ch/sprot/). 

10 

Table 2. Proteins identified using peptides isolated from human lung tissue. 

P47915 60s ribosomal protein 129. 5/2000 [MASS=17456] 
P48025 tyrosine-protein kinase syk (ec 2,7.1.1 12) (spleen tyrosine kinase). 11/1997 
P48147 prolyl endopeptidase (ec 3.4.21 .26) (post-proiine cleaving enzyme) (pe). 10/1 
P48444 coatomer delta subunit (delta-coat protein) (delta-cop) (archain). 11/1997 [M 
P48634 large proline-rich protein bat2 (hla-b-associated transcript 2). 2/1996 [MASS 
P48735 isocitrate dehydrogenase [nadp], mitochondrial precursor (ec 1.1.1 .42) (oxaio 
P49023 paxillin. 7/1998 [MASS=60937] 

P49137 map kinase-activated protein kinase 2 (ec 2.7.1.-) (mapk-activated protein ki 
P49182 heparin cofactor ii precursor (hc-ii) (protease inhibitor leuserpin 2). 11/19 
P49321 nuclear autoantigenic sperm protein (nasp). 7/1998 [MASS=85191] 
P49327 fatty acid synthase (ec 2.3.1 .85) [includes: ec 2.3.1 .38; ec 2.3.1 .39; ec 2.3 
P49407 beta-arrestin 1. 7/1999 [MASS=46969] 

P49411 elongation factor tu, mitochondrial precursor (p43). 12/1998 [MASS=49542] 

P49773 hint protein (protein kinase c inhibitor 1) (pkci-1). 7/1998 [MASS=1 3671] 

P50096 inosine-S'-monophosphate dehydrogenase 1 (ec 1.1.1.205) (imp dehydrogenase 1) 

P50552 vasodilator-stimulated phosphoprotein (vasp). 11/1997 [MASS^39830] 

P50748 hypothetical protein kiaa0166. 1 1/1997 [MASS=250749] 

P50851 cdc4-like protein (fragment). 7/1998 [MASS=213599] 

P51174 aeyl-coa dehydrogenase, long-chain specific precursor (ec 1.3.99.13) (lead). 

P51660 estradiol 17beta-dehydrogenase4 (ec 1.1.1.62) (17-beta-hsd 4) (17-beta-hydr 

P51790 chloride channel protein 3 (clc-3). 7/1998 [MASS=84793] 

P51812 ribosomal protein s6 kinase ii alpha 3 (ec 2.7.1.-) (s6kii-alpha 3) (p90-rsk 

P51885 lumican precursor (lum) (keratan sulfate proteoglycan). 7/1998 [MASS=38351] 

P51991 heterogeneous nuclear ribonucleoprotein a3 (hnrnp a3) (fbrnp) (d10s102). 7/19 

P52272 heterogeneous nuclear ribonucleoprotein m (hnrnp m). 10/1996 [MASS=77469] 

P52480 pyruvate kinase, m2 isozyme (ec 2.7.1.40). 7/1999 [MASS=57756] 

Cursory examination of this list shows that many interesting and therapeutically 
important proteins are identified by this process, including low abundance 
15 regulatory proteins such as signaling proteins, transport channels, and nuclear 
proteins. 
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A common criticism of current proteomics technologies based on two- 
dimensional polyacrylamide gels is that they are insensitive and only identify high 
abundance metabolic proteins, ie. proteins that are not normally critical 
determents of disease (although these can be important effectors of disease) 
especially since drug development strategies nearly always target low abundance 
proteins important for counteracting a disease phenotype. 

It is clear from the above table that peptide profiling can successfully describe 
many proteins that are considered important therapeutic targets, and not just 
metabolic enzymes and structural proteins. 

Table 3 shows how proteins from various therapeutically important categories 
were readily identified and quantified in a single analysis. This list was made 
using keywords present in the sequence annotation databases and therefore 
represents the minimum representation of such classes - the vast majority of 
sequenced mammalian proteins await functional annotation. 

By contrast, a recently published study (Proteomics 1 ,1303-1 9 A database of 
protein expression in lung cancer.Oh JM, Brichory F, Puravs E, Kuick R, Wood C, 
Rouillard JM, Tra J, Kardia S, Beer D, Hanash S. 2001) where over 1300 2D gels 
were analyzed from a variety of different lung cell lines and tumors, identified less 
than 200 proteins, the majority of which were metabolic and structural proteins of 
high abundance, and provided no quantitative information. 



Table 3. Peptide profiling identifies therapeutically important proteins. 
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Conventional approach (Oh et 
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Kinases j 
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Phosphatases 




12 


1 
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Channel proteins 




12 
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Apoptosis proteins 


1 
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Proteins contributing to cancer : 


10 


0 








Proteins with homoloqv to viral proteins 


27 


0 
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Example 4. Peptide profiling to characterize diverse human tissues 



One-dimensional LCMS was used to obtain peptide profiles from diverse human 
tissues (Fig.5). The one-dimensional approach has 2- to 10-fold lower resolution 
5 compared to two-dimensional approaches but was used in this case to example a 
large number of samples to illustrate the principle. Table 4 shows the number of 
peptides and proteins identified for different human tissues. 

Table 4. The peptide profiling approach can be applied to diverse tissues. 

10 





Proteins 


Peptides 


Brain 


359 


734 


Heart 


114 


231 


Testes 


78 


136 


Liver 


56 


83 


Muscle 


72 


66 


Plasma 


288 


846 


Pancreas 


202 


283 



J I 



It is assumed that diverse tissues may express many similar proteins (for 
instance ribosome associated proteins), yet express a subset of unique proteins 

15 that functionally distinguishes one tissue from another. Similarly, the proteome of 
diseased tissue may be different to healthy tissue. Although this may seem self- 
evident, very few studies have addressed these issues by directly comparing the 
proteomes from different samples. This is largely because of the technical 
impediments mentioned above - conventional techniques generally characterize 

20 only the most abundant proteins and peptides, and these peptides are least likely 
to differ from tissue to tissue. Figure 6 shows how many proteins were identified 
using MCAT based peptide profiling for a preliminary study of seven human 
tissues. Notably, the peptide and protein profiles of each tissue is distinct. Even 
with this preliminary low resolution analysis, each tissue evokes a different 

25 signature when subjected to peptide profiling. 

When the proteins identified for different tissues are compared, it is clear that 
some proteins are common to several tissues, while some are tissue-specific 
(Fig. 6). These differences can be highlighted by applying agglomerative 
30 clustering algorithms to the data (Fig. 7). In this figure as an example, common 
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proteins are highlighted in the large rectangular box, while heart- and brain- 
specific proteins are highlighted in the smaller rectangular boxes. 
Furthermore, the degree of relationship between these tissues can be 
established by comparing such peptide profiles (Fig. 8). Although the principle 
was illustrated here Using different human tissues, such analysis can be used to 
detect other proteomic changes, for instance human heart tissue following 
exercise or myocardial infarction, or following administration of drugs. 

Example 4. Peptide profiling to chamrt*ri 7 e subcellular factions of a Wn^ fa 



io tissue 



In another embodiment of the invention, peptide profiling can be used to analyze 
the subtractions of a cell, preferably into nuclear, cytoplasmic and membrane 
fractions. This discriminatory power of peptide profiling is illustrated here, where 
the method is used to examine the subtractions of a single clonal cell line. 
Cultured human myoblast cells were processed into nuclear, cytoplasmic and 
membrane fractions and analyzed using the peptide profiling technique (Fig. 9). 
Significantly, over 400 membrane-localized proteins were identified. This class is 
normally very difficult to analyze using conventional proteomics methods yet is of 
particular pharmacologic/therapeutic interest, being the site of receptors and 
channels with critical signaling and transport functions. 

Tables 7 and 8 show how peptide profiling can be applied to different cellular 
subtractions and used to identify compartment-specific proteins. 

Table 7. Peptide profiling applied to different cell compartments. 





Peptides 


Proteins 


Cytoplasmic 


2220 


994 


Nuclear 




804 


428 


Membrane 




727 


403 



( 
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Table 8. Peptide profiling identifies compartment-specific proteins 





CytoplasmicMembrane 


Nuclear 


Unique 


805 


249 


262 


Total 


994 


428 


403 


Percent unique 


80 


58 


65 



5 Example 5: Use of peptide profiles to characterize human cell lines 

In another embodiment of the invention, this invention includes methods of 
characterizing human cell lines. The method comprises generating samples 
suitable for MS analysis and producing a peptide profile. The relative abundance 
of peptides in samples is also preferably determined. The peptide profile that is 

10 generated is compared to peptide profiles in a database or library using common 
algorithms in order to identify cognate proteins, preferably those that are 
considered important therapeutic targets, as well as metabolic enzymes and 
structural proteins. In a further embodiment, these profiles can comprise a small 
prototype database or library, against which novel samples may be screened. 

15 

A number of peptides from four human cell lines of distinct cellular origin are 
identified by mass spectrometry and linked to their parent proteins. This profile is 
one-dimensional because no addition information about the peptides (e.g. 
quantitative information) is included. Table 6 shows the number of peptides and 
20 proteins identified for the different human cell lines. 

Table 6. Peptide profiling of different cultured cells 





Proteins 


Peptides 


Myoblasts 


576 


1373 


HeLa 


974 


2067 


NYP17 


192 


290 


Raji-Jurka 


233 


376 



25 

Here, an independent extract of one of the four cell lines is screened and 
demonstrates how this extract can be conclusively shown to be highly similar or 
identical to a profile in the database. 

30 Method 

Cell extracts derived from four human cell lines (MCF7, TPA, Jurkat, K566) were 
digested with trypsin (Porozyme, Perceptive Biosystems, USA) and analyzed 

54 



ISDOCID: <WO_02097703A3JA> 



10 



15 



WO 2002/097703 

PCT/CA2002/000801 

using an ion trap mass spectrometer (Deca, Thermoquest, USA) following 
separation of digested peptides using online HPLC. The mass spectrometer was 
programmed to collect primary MS spectra from parent ions, as well as tandem 
mass spectra of daughter ions generated from the first, second and third most 

abundant ions observed in the program window. These spectra were then used 
to search nonredundant genome databases using the SEQUEST algorithm 
(Yates et a.., 1995) to identify the peptides and proteins present in the samples. 

The following table shows the top-scoring peptides identified in the analysis of 
one of these cell lines, Jurkat: experiment, After statistical filtering 74 91 96 
123 peptides were used to identify 55, 62, 49, 59 different proteins in the ' 
respective cell. The peptides for all four cell lines were deposited into a database 
m th«s case a Microsoft Access file. The protein profiles are graphically 
represented below (5922, 4091, 5644 and 4166 tryptic peptides were observed 
from MCF7, TPA, Jurkat and K566 cells respectively, visual representation)- 
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If these profiles are considered as a small Index or database, novel profiles can 
be searched against them using any common correlation test. For instance here 
the correlation is calculated by: 

Px,y = [1/0 fl =i to n) I (Xj - |Jx) (Yj - |j y )]/[a x . d y ] 



where peptides common to two profiles score T and peptides not shared 
between profiles score '0\ 
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Correlation scores, P x , y , for one-dimensional peptide profiles obtained from four 
human cell lines: 





MCF7 


TPA 


Jurkat 


K556 


? 


MCF7 


1 


0.0105 


0.33596 


0.09 


0.07 


TPA 


0.0105 


1 


0.33596 


0.31714 


0.26733 


Jurkat 


0.33596 


0.33595 


1 


0.09 


.8644 


K556 


0.09 


0.31714 


0.09 


1 


.0.09 



5 „ Tnis P re,imina, y analysis suggests that the.peptide.profiles obtainedfrom Jurkat 
and MCF7, and Jurkat and TPA nuclear extracts are more similar than those 
obtained for other combinations. More importantly, when the peptide profile 
obtained from an independent preparation of Jurkat nuclear extract (labeled '?' in 
the above Table), it received a high score and could be identified as being most 

10 closely related to the Jurkat cells. 

Applications of Protein Expression DatasPte 



( 



15 



20 



25 



Relevance to Disease 

As an example of the approach, its potential use in the diagnosis and study of 
human disease is described, for example in infectious disease or a genetic 
disease=such=as-caneer^TheHnvent^ 

compare, classify, and characterize and investigate biological or clinical samples 
from normal and virus- or bacterially-infected cells and tissues, similar cells 
obtained over a course of infection, or similar cells obtained over the course of a 
therapeutic treatment. Similarly, the invention may be used to systematically 
identify, compare, classify, and characterize and investigate biological or clinical 
samples from normal and cancerous cells and tissues, cancerous cells and 
tissues obtained from a variety of related or unrelated liquid or solid tumors, cells 
obtained over time that follow the development of a progressive cancer, or cells 
similarly obtained over time that follow the progression of a therapeutic 
intervention. 



c 
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The resulting datasets or profiles may therefore (i) identify robust signatures of 
disease states that can be used to facilitate diagnostic and prognostic medical 
procedures, (ii) refine current models of disease and highlight productive areas 
for focusing further basic and applied investigative approaches. 

5 

Uses in Toxicology Studies 

As another example of the use of the invention, quantitative peptide profiles may 
be used for investigation of toxic effects in human or other tissues or cells, for 
instance the side-effects of candidate drug compounds. This is because the 

10 toxicity may be represented by changes in the expression patterns of peptides 
and proteins in the cells. Currently, such toxic effects are investigated using 
general marker enzymes such as cytochrome oxidase. In many ways, this is a 
'blunt tool', failing to differentiate between different types of toxicity, and/or the 
severity of the toxic effect. Quantitative peptide profiles are likely to be discrete 

15 for individual compounds while profiles generated in response to related 
compounds would be expected to be also related to each other. 

A database of profiles can be assembled that describes the protein complements 
of tissues treated with known toxic agents. Large numbers of drug candidates 

20 can then be screened and their profiles compared to those in the reference 
database. Accordingly, the invention includes methods of determining the toxicity 
of a candidate drug compound. The method comprises administering the 
candidate compound to a cell. As described above, samples suitable for MS 
anaylsis are generated and a peptide profile is produced. Relative abundance of 

25 peptides in samples is also preferably determined. This candidate compound 
peptide profile is compared to peptide profiles in a database or library (for 
example, profiles showing the cell in a normal state and in varied states of 
toxicity). If the candidate compound sample profile is highly similar to (for 
example, greater than 90%, 95%, or 99% similarity), or identical to a profile in the 

30 database or library, then that similarity shows the amount of toxicity of the 
candidate compound to the cell. If the candidate compound sample profile is 
highly similar to a normal cell profile, then the candidate compound is less likely 
to be toxic than if the candidate compound sample profile is similar to the peptide 
profile of the cell in state of toxicity. The relative abundance of the test sample 
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peptides is also preferably compared to other profiles to determine the amount of 
toxicity of a candidate compound. 

Profiles obtained from drug candidates that are similar to those obtained from 
5 damaged tissue alert the investigators to potential toxicity problems associated 
with that compound. Because each single profile comprises a large dataset 
(many individual proteins and their relative abundances), comparison of the 
profiles is statistically powerful. This reduces dependence on animal toxicity 
trials, where large numbers of animals may be necessary to obtain statistically 
10 relevant data. 

Healthy cells, and cells treated with toxic agents, will be analyzed by liquid 
chromatography-tandem mass spectrometry (LC-MS/MS) using a novel semi- 
quantitative approach, resulting in a protein profile for each treatment that serves 
as a signature of the cell state. The profile comprises data relating tens to 
hundreds of individual proteins and therefore represents a highly specific and 
sensitive description of the protein complement of the cell or tissue in that 
particular state. 



15 



20 



25 



Even without knowledge of protein function, the profiles from cells treated with 
novel compounds can be compared to those from healthy cells or cells treated 

-with-texfe^p^^ 

an early stage of drug development. Further, where the test profile matches the 
profile produced by treatment with a characterized compound or family of 
compounds, the mechanism of toxicity may be similar to that produced by the 
reference class. This application of the invention can be applied to any primary or 
transformed cell line, or to tissues obtained from animal models, preferably 
mammalian and more preferably human, or to experimental or clinical samples. 



30 




A further embodiment of the peptide profiling invention is to characterize and 
identify the effect of drugs and other experimental treatments on the proteome 
In this example, cultured human muscle cells were treated with the hormone drug 
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leptin. For both treated and untreated samples, over 400 proteins and 900 
peptides were identified. Of these, 170 were uniquely observed in one or other 
sample. In Figure 10, a screenshot of this analysis shows peptides present in 
one or other sample (green or red) and peptides unique to either sample (blue). 
5 This experiment demonstrates that the invention can be used to examine the 
effect of drugs and other treatments on proteome mixtures. 

Example 7. Peptide profiling to characterize tissue from different organisms 

10 As further proof of principle, the peptide profiling approach was applied to 
different organisms - two microbes (Escherchia coli and Saccharomyces 
cerevisiae) and two mammals (Homo sapiens - humans and Mus musculis - 
common lab mouse). A standard MCAT LC-MS peptide profiling analysis was 
used to follow expression of hundreds of proteins for each species (Tables 9, 10). 

15 

Table 9. Peptide profiling of microbial species. 





Proteins 


Peptides \ 


Yeast 


233 


519 


Bacteria 


542 


1647 



20 

When the peptide profiles of the highly divergent microbial species were 
compared, 516 of the 519 yeast proteins were unique. In contrast, when a similar 
analysis was done for peptide profiles of the two mammalian species, 44 of 1 97 
mouse peptides were similarly observed in the human profile (representing 
25 homologous protein/peptide species). Thus, these preliminary analyses indicate 
that peptide profiling can both distinguish species, and that the peptide profile 
may reflect the degree of relatedness of organisms (Figure 1 1). 

Table 10. Peptide profiling of mammalian species. 

30 





Proteins 


Peptides \ 


Mouse 


142 


197 


Human 


256 


445 



Example 8. Peptide profiling is reproducible 
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Because peptide profiling relies on the use of many data points to assess the 
degree of relatedness of many different samples, it is critical that the method be 
reproducible. This is confirmed on the samples described here. One such 
example, involving the peptide profile of yeast whole cell lysate, is shown here 
(Table 11, Table 12). 



Table 1 1. Peptides observed for two repeat samples. 





? Total" ! Shared 


Samplel 


! 776j 686 


Sample2 


723! 686 



Table 12. Proteins observed for two repeat samples. 



! Total Shared 


Samplel 


304. 259 


Sample2 


1 288! 259 



This analysis establishes the reproducibility of the process. 

Figure 12 is a representation of a reference database of protein profiles, 
incorporating both the identity, relative quantities, and overlap of peptides or 
proteins in various samples. 

It will be appreciated that the description above relates to the preferred 
embodiments by way of example only. Many variations on the computer system 
and methods for delivering the invention will be obvious to those knowledgeable 
In the field, and such obvious variations are within the scope of the invention as 
described and claimed, whether or not expressly described. 

All references, including journal articles, patents and patent applications, in this 
application are incorporated by reference herein in their entirety. 
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We claim: 

1 . A method for identifying the constituent proteins for a cell type, tissue or 
pathological sample using a database comprising peptide profile libraries 
wherein the libraries have multiple peptide sequences, comprising: 

5 a) deriving a plurality of peptides from the cell type, tissue or pathological 
sample; 

b) identifying the peptide species by liquid phase tandem mass spectroscopy 
sequencing; 

c) compiling a data set or peptide profile containing the collection of peptide 
10 sequences obtained thereby; and 

d) cross-tabulating with a collection of peptide sequences in the database. 

2. The method of claim 1 , wherein the step of deriving a plurality of peptides 
from the cell type, tissue or pathological sample further comprises the step 

is of: 

a) obtaining a peptide-containing extract of the cell type, tissue or 
pathological sample; 

b) digesting the extract producing peptides with an enzyme, the enzyme 
capable of localizing mobile protons to the N-terminal amine and the side 

20 chains of the carboxy-terminal arginine or lysine residues; 

c) separating the peptides by high pressure liquid chromatography 
apparatus; ~~~ " ' 

3 - The method of claim 2, wherein the enzyme comprises one selected from 
25 the group consisting of trypsin and endoproteinase LysC. 

4. The method of any of claims 2 to 3, wherein the step of digesting the 
extract producing peptides further comprises the steps of: 
a) dividing the extract into two equal portions; 
30 b) derivatizing one of the two equal portions with a reagent, the reagent 
comprising one selected from the group consisting of o-methylisourea, 
homoarginine, canavanine, hydrazine, phenylhydrazine, and butyric acid 
derivatives. 

c) combining the two portions. 
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5. A method for identifying a peptide sequence for a cell type, tissue or 
pathological sample using a database comprising peptide profile libraries 
wherein the libraries have multiple peptide sequences, comprising: 

5 a) obtaining a peptide-containing extract of the cell type, tissue or 
pathological sample; 

b) digesting the extract producing peptides with an enzyme capable of 
localizing mobile protons to the N-terminal amine and the side chains of the 
carboxy-terminal arginine or lysine residues; 
10 c) separating the peptides by high pressure liquid chromatography 
apparatus; 

d) identifying the peptide species by tandem mass spectroscopy sequencing; 
and 

e) compiling a data set or peptide profile containing the collection of peptide 
15 sequences obtained thereby. 

6. The method of claim 5, wherein the enzyme comprises one selected from 
the group consisting of trypsin and endoproteinase LysC. 

20 7. The method of any of claims 5 to 6, wherein the step of digesting the 
extract producing peptides further comprises the steps of: 

a) dividing the extract into two equal portions; 

b) derivatizing completely one of the two equal portions with a reagent, the 
reagent comprising one selected from the group consisting of o- 

25 methylisourea, homoarginine, canavanine, hydrazine, phenylhydrazine, and 

butyric acid derivatives. 

c) combining the two portions. 

8. A method for quantitating the relative abundance of proteins in two 
30 samples of a cell type, tissue or pathological sample using a database comprising 
peptide profile libraries wherein the libraries have multiple peptide sequences, 
comprising: 

a) deriving a plurality of peptides from each sample of the cell type, tissue or 
pathological sample; 
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b) identifying the peptide species by tandem mass spectroscopy sequencing; 
compiling a data set or peptide profile containing the collection of peptide 
sequences obtained thereby; 

c) cross-tabulating with a collection of peptide sequences in the database of 
5 peptide sequences; and 

d) determining the relative abundance of the peptides and/or proteins. 

9. A method for quantitating the relative abundance of proteins in two samples 
of a cell type, tissue or pathological sample using a database comprising 

10 peptide profile libraries wherein the libraries have multiple 

sequences, comprising: 

a) deriving a plurality of peptides from each sample of the cell type, tissue or 
pathological sample; 

b) identifying the peptide species by tandem mass spectroscopy sequencing; 
15 c) compiling a data set or peptide profile containing the collection of peptide 

sequences obtained thereby; 

d) determining the degree of relatedness of a collection of peptide sequences 
in the database of peptide sequences using clustering and related statistical 
methods 

20 

10. The method of any of claims 8 to 9, wherein the step of deriving a plurality of 
peptides in Itwo samples further comprises the step of: 

a) obtaining a peptide-containing extract of each sample; 

b) digesting separately the extracts producing peptides with an enzyme, the 
25 enzyme capable of localizing mobile protons to the N-terminal amine and the side 

chains of the carboxy-terminal arginine or lysine residues; 

c) combining the two extracts; and 

d) separating the peptides by high pressure liquid chromatography. 

30 11. The method of claim 10, wherein the enzyme comprises one selected from 
the group consisting of trypsin and endoproteinase LysC. 

12. The method of any of claims 9 to 1 1 , wherein the step of digesting the 
extracts further comprises the step of derivatizing completely one of the two 
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extracts with a reagent, the reagent comprising one selected from the group 
consisting of o-methylisourea, homoarginine, canavanine, hydrazine, 
phenylhydrazine, and butyric acid derivatives. 



5 13. A method for identifying a peptide sequence for a cell type, tissue or 
pathological sample, comprising: 

a) obtaining a peptide-containing extract of a cell type, tissue or pathological 
sample; 

b) digesting the extract producing peptides with an enzyme capable of 
10 localizing mobile protons to the N-terminal amine and the side chains of the 

carboxy-terminal arginine or lysine residues; 

c) separating the peptides by high pressure liquid chromatography 
apparatus; 

identifying the peptide species by tandem mass spectroscopy sequencing; and 
15 d) compiling a data set or peptide profile containing the collection of peptide 
sequences obtained thereby. 

14. The method of claim 13, wherein the enzyme comprises one selected from 
the group consisting of trypsin and endoproteinase LysC. 

20 

15. The method of any of claims 12 to 14, wherein the step of digesting the 
extract producing peptides further comprises the steps of: 

a) dividing the extract into two equal portions; 

b) derivatizing completely one of the two equal portions with a reagent, the 
25 reagent comprising one selected from the group consisting of o- 

methylisourea, homoarginine, canavanine, hydrazine, phenylhydrazine, and 
butyric acid derivatives: 

c) combining the two portions. 



30 16. A computer system for identifying quantitative peptide profiles, comprising: 
a) a database including peptide profile libraries for a plurality of types of 
organisms wherein the libraries have multiple peptide profiles each profile 
comprising an array of at least 50 peptide species each having a unique identifier 
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cross-tabulated with quantitative data indicating relative and/or absolute 
abundance of each peptide species in a sample; and 

b) a user interface capable of receiving a selection of one or more queries to 
the database for use in determining a rank-ordered similarity of peptide profiles in 
the database. 



17. A method of producing a computer database comprising a computer and 
software for storing in computer-retrievable form a collection of peptide profiles 
for cross-tabulating with data specifying the source of the peptide-containing 

10 sample from which each peptide profile was obtained . 

18. The method of claim 17, wherein at least one of the sources is from a 
sample known to be free of pathological disorders. 

15 1 9. The method of claim 1 8, wherein at least one of the sources is a known 
pathological specimen. 



. 20. A method of comparing quantitative peptide profiles using a database of a 
plurality of peptide profile libraries, the method comprising: 
20 a) receiving a selection of two or more of the peptide profile libraries; 

b) determining the peptide profiles common to the selected peptide profile 
libraries and identifying profiles unique to each of selected peptide profile library; 
and 

c) displaying the results of the determination. 

25 

21 . The method of claim 20, wherein the correlation of a peptide profile 
against selected peptide profile libraries is determined by 

Px,y = [1/n 0=! to n) I (Xj - Ux) (Yj - Uy)]/[d x . dy ] 

where peptides common to two profiles score '1' and peptides not shared 
30 between profiles score '0'. 

22. The method of claim 21 , wherein the peptides profiles are of cell fractions, 
the cell fractions comprising high molecular weight proteins, soluble proteins, 
membrane proteins, modified proteins, phosphoproteins, peptides terminating in 
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lysine or arginine or the specific products of proteolytic enzymes or chemical 
derivatives of those products, peptides containing rare amino acids, and proteins 
isolated by binding to disease-specific affinity reagents. 

5 23. The method of claim 22, wherein the specific products of proteolytic 
enzymes comprise chemical derivatives of these products wherein de novo 
sequencing or relative abundance measurements of the peptides is facilitated. 

24. The method of claim 23, wherein the chemical derivatives are obtained by 
10 guanidinylation and related modifications. 

25. The method of any of claims 21 to 24, wherein the rare amino acids 
comprise tryptophan and cysteine and amino acids comprising 5% or less of the 
amino acid representation. 

15 

26. The method of any of claims 21 to 25, wherein the disease-specific affinity 
reagents comprise polyclonal antibodies, toxin or drugs. 

27. The method of any of claims 21 to 26, wherein the peptide profiles are of 
20 peptide sequences, the peptide sequences comprising mammalian peptide 

sequences. 

28. The method of any of claims 21 to 26, wherein the peptide profiles are of 
peptide sequences, the peptide sequences comprising microbial peptide 

25 sequences. 

29. The method of any of claims 21 to 28, wherein the step of receiving a 
selection of two or more of the peptide profile libraries for comparison includes 
receiving a user selection from two or more pull-down menus using a graphical 

30 user interface. 

30. The method of any of claims 21 to 28, wherein the step of receiving a 
selection of two or more of the peptide profile libraries for comparison comprises 
command line entry using a computer. 
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31 . The method of any of claims 21 to 28, wherein the step of receiving a 
selection of two or more of the peptide profile libraries for comparison comprises 
receiving an electronically transmitted file containing sequence and quantitative 

5 data. 

32. The method of any of claims 21 to 31 , wherein the results of the 
determination comprise a unique identifier for related peptide profiles. 

10 33 - The method of any of claims 21 to 31 , wherein the results of the 

determination comprise annotated information relating to the related peptide 
profiles obtained from a public database. 

34. The method of any of claims 21 to 31 , wherein the results of the 

is determination comprise quantitative or relative abundance information relating to 
the related peptide profiles obtained from a public database. 

35. The method of any of claims 21 to 34, further comprising the step of 
displaying the peptide profiles common to the selected peptide profile libraries. 

20 

36. The method of any of claims 21 to 34, further comprising the step of 
displaying the peptide profiles unique to the selected peptide profile libraries. 

37 - A method of identifying peptide profiles common to a set of environments, 
25 organisms, organs, tissues, cells, cellular fractions or isolated molecular 

complexes using a database comprising peptide profile libraries for a plurality of 
types of organisms wherein the libraries have multiple peptide sequences, the 
method comprising: 

a) displaying at least one list of peptide profile libraries; 
30 b) receiving a selection of one or more peptide profile libraries from at least 
one list of peptide profile libraries; 

c) determining peptide profiles common to the selected peptide profile 
libraries; and 

d) displaying the results of said determination. 
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The peptide or protein profile of different tissues are unique. 
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Bioinformatics techniques reveal general and tissue-specific proteins. 
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Similarity dendrogram for different human tissue constructed 
using peptide profiling. 
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Note: Fract refers to a chromatographic subtraction. 
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Comparison of peptide profiles for untreated- and leptin-treated human 
muscle cells. 
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